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Ever hear of 
irty-five header? 


~ 
Tandem locomotives, teamed to pull a train — 
up a grade, are often called a double header. 
It would take a “35 header’? to equal the 
power of the RMI rocket engine shown here. 
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RMI propulsion systems are now being per- RMI since 1941 . . . reliability in application. 
fected for use in ballistic and guided missiles Another reason why RMI’s program of 
as well as in the manned aircraft of the advanced research and engineering contin- 
future. These advanced propulsion systems ues to be a good investment for taxpayers, 
incorporate the principle which has guided government and industry! s ] \ 
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As the technology of missile machining progressed missilemen looked forward to oO 
the many advantages of combining hydrospinning and contour finish machining. - 
But the technical obstacles were great and no one dared the expensive risk of i 
research. 
Diversey tackled this technical research and now after many months of experi- 
mentation Diversey has solved this problem. Above you see the successful produc- 
tion line results. You are looking at a Diversey hydrospun nose section for an IRBM 
being contour finish machined to extremely close tolerances. SEND 
os This advance of Diversey is important to the missile industry because of its FOR 
flexibility, because it saves valuable metal, and because it sharply reduces cost. 
At Diversey you have the LARGEST FACILITIES for your missile metal machining FREE 
and hydrospinning problems. You have a team of alert, experienced, BOOKLET 
and progressive technical people. Bring your tough jobs to Diversey. © 
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Artificial Satellites—A Bibliography of Recent Literature 


_— annotated listing of about 340 refer- 
ences,! alphabetical by author, brings up 
to date the bibliography on *‘A History of the 
Artificial Satellite’ by Alan R. Krull, pub- 
lished in JeT Proputsion, May 1956, pages 
369-383. 

Since, with the exception of one or two 
1956 entries, the Krull bibliography covered 
the literature through 1955, this listing pur- 
posely includes references selected from 
books, unclassified research reports and tech- 
nical periodical articles published during 
1956, 1957 and January-February 1958. 

Much has been written since 1955, for 
much has transpired to make the artificial 
satellite one of the chief topics of the times. 
Instead of theorizing, man has actually 
launched artifical satellites, so this bibliog- 
raphy is really a recorded reference story of 
Vanguard, Sputnik and Explorer, their 
progress, development and scientific uses. 

The same guide lines mentioned by Mr. 
Krull in his earlier bibliography have been 
followed in this one; that is, fiction, juvenile 
literature and news items are omitted. 
There are a few entries connected with space 
law, but no attempt has been made to in- 
clude everything; also, there is nothing on 
space medicine, manned flight, interplane- 
tary flight, trips to the moon or to Mars or 
space stations. The following bibliographies 
are suggested for further information on 
these phases: 

Benton, Mildred. “Earth Satellites, 
Guided Missiles, Rockets and Space Flight,” 
Wilson Lib. Bull. 32: 412-419, Feb. 1958. 

Department of the Army Library. ‘“‘Guided 
Missiles, Rockets and Artificial Satellites 
(Including Project Vanguard),”’ 153  pp., 
Jan. 1957 (Spec. Bib. 11). 

Hogan, J.C. ‘‘Space Law Bibliography,” 
J. Air Law & Commerce 23: 317-325, 
Summer 1956. 

Krzywoblocki, M. Z. ‘‘The Forces Acting 
on an Air Vehicle. A Review of the Litera- 
ture,”’ 112 pp., Chicago University, June 
1955 (Aerodynamic Studies, Suppl. 10). 

Ley, Willy. ‘‘Rockets, Missiles and Space 
Travel,” 528 pp., New York, Viking, 1957 
(Bibliography, pp. 489-520). 

‘Missiles and Myths,” 
Rockets 2: 74, Aug. 1957. 

National Science Foundation. ‘A Bib- 
liography for the International Geophysical 
Year,’ 51 pp., Washington, D. C., July 1, 
1957 (NSF-57-25). 

Rosen, M. W. “Twenty-Five Years of 
Progress Toward Space Flight’’ (includes 
literature review). JET Propu.sion 25: 
623-626, Nov. 1955. 


Missiles and 
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' Part Two of this Bibliography covering 
the period from January 1957 through 
February 1958 will appear in a subsequent 
issue of JET PROPULSION. 
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‘Part One—1956 


MILDRED BENTON? 


~ 
U. S. Naval Research Laboratory, Washington, D. C. 


Anderton, D. A. 


Gas Dynamics, Part I-—-Key Tool for 
ICBM, Satellite Studies. Part Il1—Hyper- 
sonic Quest Yields New Techniques. Missile 


Eng. 1: 4-11, illus., Sept. 1956. A discussion 
of the two extremes of aerodynamic flow 
problems—hypersonics (ICBM) and super- 
aerodynamics (IGY satellite). 

Vanguard Depends on Tools and Gages. 
Aviat. Wk. 65: 50-52, Nov. 5, 1956. Tools 
and gages will determine whether the Van- 
guard will be a satellite or a meteor. 


Army Signal Corps Engineering Labora- 
tories, Fort Monmouth, N. J 

The Environment of an Earth Satellite, by 
R. Griffith, W. Norberg and W. G. Stroud. 
41 pp., diags., 1956. (Tech. Memo. 1747). 
Graphs, tables and other data summarizing 
what is known of the composition and proper- 
ties of the upper atmosphere and the pres- 
ence of radiation, gravity, geomagnetism and 
meteorities at high altitudes. 

Baum, W. A. 

A Fundamental Cosmological Experiment 
for the Artificial Satellite. Astron. Soc. Pac. 
Publ. 68: 118-120, Apr. 1956. Information to 
be sought by means of instruments is bright- 
ness and color of the extragalactic sky. 
Bennett, W. H. 

Proposed Measurement of Solar Stream Pro- 
tons. In Van Allen, J. A., ed., Scientific 
Uses of Earth Satellites, pp. 194-197, Ann 
Arbor, University of Michigan Press, 1956. 
The solar protons in their initial approach in 
free orbits and any protons traveling later in 
captive orbits can be detected and measured 
only in the vicinity of the auroral zones of 
magnetic latitude, and an artificial satellite 
which passes approximately over the two 
poles of the earth will be needed for this pur- 
pose. 

Bergaust, Erik 

Satellite Science Not So Simple. U. S. 
Satellite Designs Copied? Missiles and 
Rockets 1: 18-19, Nov. 1956. Includes 
illustration from a Russian magazine which 
indicates what the Russian public think their 
IGY satellite will look like. It is an exact 
copy of 8. F. Singer’s early MOUSE (1954). 
Bergaust, Erik, and Beller, William 

Satellite. 287 pp., illus., New York, Han- 
over House, 1956. A general nontechnical 
discussion of the projected program for the 
launching of earth satellites—their operation, 
their possible uses and their likely effects on 
earth. 

Berning, W. W. 

Tonospheric Structure as Determined by a 
Minimal Artificial Satellite. In Van Allen, 
J. A., ed. Scientific Uses of Earth Satellites, 
pp. 253-262, Ann Arbor, University of 
Michigan Press, 1956. Discusses the prob- 
lems and develops the analytical methods 
necessary for determining ionospheric elec- 
tron densities from trajectory measurements. 


Blitzer, Leon, Weisfeld, M., 
A. D. 

Perturbations of a Satellite’s Orbit Due to— 
Earth’s Oblateness. J. Appl. Phys. 27: 1141- 
1149, Oct. 1956. Points up possible use of © 
the US-IGY satellites for determining the — 
earth’s oblateness to a new precision. (For — 
extension of this paper, see Blitzer entry in 
Part Two, 1957.) 


Chubb, T. A., Friedman, H., and Kupperian, — 


and Wheelon, © 


= 


A Lyman Alpha Experiment for the Van- 
guard Satellite. In Van Allen, J. A., ed. 
Scientific Uses of Earth Satellites, pp. 147- 
151, Ann Arbor, University of Michigan 
Press, 1956. <A satellite experiment 
the possibility of monitoring the sun 
Lyman alpha almost continuously for many 
days, and would be worth hundreds of 
individual rockets for the same purpose. 

A Satellite Experiment to Determine the 
Distribution of Hydrogen in Space. In Van 
Allen, J. A., ed. Scientific Uses of Earth 
Satellites, pp. 152-156, Ann Arbor, ine 
sity of Michigan Press, 1956. Proposed use 
of a photon counter to measure simultane- 
ously the intensity of hydrogen Lyman alpha 
radiation (1215.7 A) received directly from 
the sun and the resonance radiation of the 
same wave length produced by sunlit hy -dro- 
gen atoms in space. 


Cornig, R. 

Economics of Rocket Propelled Aeroplanes. 7 
II. Aero. Eng. Rev. 15: 49-58, Oct. 1956. 
An appendix is included which indicates costs — 
for placing material in a satellite orbit. 


Cushman, Robert 

Next Satellite Problem: Data Descent. — 
Aviat. Wk. 64: 53, 56-57, May 14, 1956. g 
Discusses satellite data recovery and men- | 
tions proposal by R. W. Porter to recover 
data from a future earth satellite by a stain- 
less steel balloon. Article with same title 
appears in Missile Eng. 1: 48-49, illus., Sept. — 
1956. 

Daniels, F. B. 

Electromagnetic Propagation Studies with a 
Satellite Vehicle. In Van Allen, J. A., e 
Scientific Uses of Earth Satellites, pp. 276-282, 
Ann Arbor, University of Michigan Press, 
1956. The proposed experiment is a deter- 
mination of the integrated value of the ioniza- 
tion from the earth’s surface to the satellite, | 
by measuring the rotation of the plane of | 
polarization of a radio wave transmitted — is 
from the satellite. ; 
Davis, R. J., Whipple, F. L., and Zirker, J.B. 

The Orbit of a Small Earth Satellite. In 
Van Allen, J. A.,ed. Scientific Uses of Earth — 
Satellites, pp. 1-22, Ann Arbor, University of | 
Michigan Press, 1956. Considers motion of | 
a satellite launched into an inclined, moder- 
ately eccentric orbit from a point approxi- 
mately 500 km above sea level. Discusses — 
effects on desired initial orbit produced by — 
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deviations from desired launching conditions 
and atmospheric drag. Additional perturba- 
tions are discussed, such as action of the sun 
and moon, and variations of the density of 
the upper atmosphere with time and latitude, 
Bernoulli forces acting on a spinning object, 
and effect of upper-atmospheric winds. 
Davis, R. J. 

Ultraviolet Stellar Magnitudes. In Van 
Allen, J. A., ed. Scientific Uses of Earth 
Satellites, pp. 157-165, Ann Arbor, Univer- 
sity of Michigan Press, 1956. Computations 
are performed from which it is possible to 
construct a tentative map of the sky at 1249 
A. Results of this survey should prove 
valuable in the selection of the most interest- 
ing regions for study when extraterrestrial 
photometry of high resolution from artificial 
satellites becomes possible. 
deBey, L. G. 

Systems Design Considerations for Satellite 
Instrumentation. In Van Allen, J. A., ed. 
Scientific Uses of Earth Satellites, pp. 49-54, 
Ann Arbor, University of Michigan Press, 
1956. Discusses some of the factors which 
must be considered in establishing sound 
objectives and adequate observing methods. 


de Nike, John 

The Effect of the Earth’s Oblateness and 
Atmosphere on a Satellite Orbit. In Franklin 
Institute. Earth Satellites as Research Ve- 
hicles. Proceedings of a Symposium, Phila- 
delphia, Franklin Institute, April 18, 1956, 
pp. 79-86, Lancaster, Pa., Journal of the 
Franklin Institute, June 1956. (Mono- 
graph 2.) A member of The Martin Com- 
pany, designers and builders of the rocket 
vehicle for the Vanguard satellite, relates 
some of the problems. 

Summary of article with same title, Brit. 
Interplan. Soc. J. 16: 29-30, Jan.-Mar. 1957. 


Dubin, Maurice 

Meteoric Bombardment. In Van Allen, 
J. A., ed. Scientific Uses of Earth Satellites, 
pp. 292-300, Ann Arbor, University of Michi- 
gan Press, 1956. An experiment for the 
detection of meteoric particles entering the 
earth’s atmosphere is proposed for an earth 
satellite. 


Easton, R. L. 

Radio Tracking of the Earth Satellite. An 
Opportunity for Amateur Collaboration. QST 
40: 38-41, 134, illus., July 1956. Describes 
simplified tracking system for amateur use. 


Ehricke, K. A. 

Aero-Thermodynamics of Descending Orbital 
Vehicles. Astronautica Acta 2: 1-19, 1956. 
The theory of laminar and turbulent friction 
and heat transfer in the compressible boun- 
dary layer with and without dissociation is 
discussed with respect to the descent of 
orbital vehicles through the atmosphere. 

Ascent of Orbital Vehicles. Astronautica 
Acta 2: 175-190, 1956. Ascent into a satel- 
lite orbit along a minimum energy ellipse 
(elliptic ascent) is compared with ascent 
along an arc of an elliptic trajectory (ballistic 
ascent) where the satellite stage is trans- 
ferred into the oribit at or near the summit. 
Discussed are basic energy requirements 
(other than losses during powered ascent), 
thrust program, effect on the size of stages, 
visibility during ascent, and flight time. 

Astronautical and Space-Medical Research 
with Automatic Satellites. In Franklin Insti- 
tute. Earth Satellites as Research Vehicles. 
Proceedings of a Symposium, Philadelphia, 
Franklin Institute, April 18, 1956, pp. 25-62, 
Lancaster, Pa., Journal of the Franklin Insti- 
tute, June 1956. (Monograph 2.) Tech- 
nical and scientific research aspects in con- 
junction with technological satellites are out- 
lined and their correlation with various phases 
of manned astronautics is shown. The sur- 
vey of the space-medical research by means 
of bio-satellites also considers the relevant 
technical and scientific aspects. 

Instrumental Satellites and Instrumental 
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Comets. {nteravia 11: 960-962, Dec. 1956. 
Discussion of future developments in satel- 
lites and their scientific applications concerns 
flight performance and accuracy, tracking 
and data transmission. 

Solar Power for Space Craft. Missiles and 
Rockets 1: 44-46, illus., Nov. 1956. Hydro- 
gen-manufacturing satellites could rely on 
concentration of solar power in power collec- 
tors. The manufacturing process could be 
completely automatic, requiring only oc- 
casional human supervision and maintenance. 


Findlay, D. A. 

Electronics in the IGY Program. Electron. 
29: 138-142, Dec. 1956. Projects planned 
for the IGY will use electronic instrumenta- 
tion to sense, record and telemeter data. 
Earth satellite program will use transistorized 
oscillators for tracking and telemetering of 
data to ground from subminiature measuring 
equipment. 


Franklin Institute 

Earth Satellites as Research Vehicles. 
Proceedings of a Symposium held April 18, 
1956. 115 pp., illus., The Institute, June 
1956. (Monograph 2.) In a conference on 
rockets for peace, the discussion was de- 
voted to different kinds of research that we 
may hope to accomplish with transistors and 
rocketry. See also entries under de Nike, 
Ehricke, Lawrence, Newell, Porter, Spitzer. 


Furth, F. R. 

Project Vanguard—The International Geo- 
physical Year Earth Satellite. Aero. Eng. 
Rev. 15: 55-59, Mar. 1956. Survey of prob- 
lems, trends and practical progress in tech- 
niques of very-high-altitude research under 
the USAF-Army-Navy sponsorship in terms 
of the design and construction requirements 
of the launching rocket and the satellite, 
selection of the launching site and tracking 
stations, and the determination of the satel- 
lite orbit. 

Also in Eng. J. (Canada) 39: 1531-1536, 
Nov. 1956, and Shell Aviation News 215: 6-9, 
May 1956. 


Gast, P. R. 

Insolation of the Upper Atmosphere and of a 
Satellite. In Van Allen, J. A., ed. Scientific 
Uses of Earth Satellites, pp. 73-84, Ann 
Arbor, University of Michigan Press, 1956. 
A consideration of temperatures. 


Gatland, K. W. 

The Vanguard Project. Spaceflight 1: 15- 
28, Oct. 1956. Vehicle design, launching and 
tracking problems, with a review of rocket 
development. 


Ginzburg, V. L. 

Eksperimental’ Naya Proverka Obshchei 
Teorit Otnositel’nosti i Iskusstvennye Sputniki 
Zemli (Experimental Verification of the 
General Theory of Relativity and Artificial 
Earth Satellites). Priroda 45: 30-39, Sept. 
1956. In Russian. In a suggestion for 
measuring the gravitational displacement of 
frequency it is proposed that one measures 
not the frequency but the difference between 
the readings on clocks on earth and the read- 
ings on clocks on the satellite. Also in Usp. 
Fiz. Nauk 63: 119-122, Sept. 1957. 

Possibility of Using Artificial Earth Satel- 
lites for the Experimental Verification of the 
Theory of General Relativity. Soviet Phys. 
JETP 3: 136-138, Aug. 1956. The reception 
of the radio signals on earth from the satellite 
can be used to determine the gravitational 
shift of frequency in the earth's field. The 
second effect that may be verified with the 
help of the satellite is the precession of the 
perihelion of the orbit. Translation of 
article appearing in Zh. Eksper. Teor. Fiz. 
30: 213-214, 1956. : 


Great Britain, Royal Aircraft Establishment, 
Farnborough 

The Descent of an Earth-Satellite Through 
the Upper Atmosphere, by D. G. King-Hele. 
9 pp., figs., Aug. 1956. (Tech. Memo GW- 


277) <A theoretical study is made of the 
flight path of an uncontrolled but aero- 
dynamically stable satellite as it spirals 
down through the earth’s atmosphere from 
an initially circular orbit under the action of 
air drag. The earth and its atmosphere are 
taken as spherically symmetrical. For alti- 
tudes above about 125 n. miles a simple 
solution is found; the velocity of the satel- 
lite is independent of its size, shape, weight 
and initial altitude and equal to the oribital 
velocity appropriate to its current altitude; 
while its angle of descent, in radians, is twice 
its drag/weight ratio. Estimates of lifetirie 
are also made. 

Haber, Heinz 

Space Satellites, Tools of Earth Research. 
Natl. Geog. Mag. 109: 487-509, illus., Ay 
1956. A general and lucid explanation of 
plans for a true space vehicle. 

Hagen, J. P. 

The Exploration of Outer Space with «n 
Earth Satellite. Inst. Radio Engrs. Proc. 4}: 
744-747, June 1956. A brief account of tiie 
beginnings of Project Vanguard: Why it is 
important to launch an earth satellite, sonie 
of the problems involved, and what is to ‘ve 
gained, research-wise. 

Haley, A. G. 

Basic Concepts of Space Law. JET PrRi- 
PULSION 26: 951-957, 968, Nov. 1956. The 
unmanned earth satellite, pp. 955-957. 
Happner, J. P. 

Satellite Geomagnetic Measurements. In 
Van Allen, J. A.,ed. Scientific Uses of Earth 
Satellites., pp. 234-246, Ann Arbor Univer- 
sity of Michigan Press, 1956. Practical 
requirements and techniques for obtaining 
useful data concerning magnetic storms and 
other disturbances. 

Hartman, L. M., and Haviland, R. P. 

A Satellite Propagation Experiment. In 
Van Allen, J. A., ed. Scientific Uses of 
Earth Satellites, pp. 268-275, Ann Arbor, 
University of Michigan Press, 1956. A 
satellite ionosphere and propagation experi- 
ment is proposed which will enlist the aid of 
operators in the Radio Amateur Service. 
Particular measurements to be made at 
specially equipped stations are reviewed. 
Haviland, R. P. 
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= Testing Ai 


Basie information concerning the free jet method of 
supersonic development testing is presented. Parameters 
necessary to specify test hardware, define test facility 
capabilities, and optimize facilities for this type of testing 
are discussed. The substantial gains obtainable in a 
facility’s testable range by separately exhausting engine 
gas flow and free jet air spilled around the engine are de- 
fined, as are performance levels of various free jet diffuser 
4 nomograph is included which permits rapid 


types. 
cseeneumes of facility requirements and/or capabilities. 
Introduction 


development of supersonic of ever- 
- increasing power to higher operating speeds has created 
demand for more powerful ground test facilities. The ex- 
pense of such facilities makes necessary the optimization of 
the various components which indirectly or directly affect 
the power requirements. 

Three principal methods of engine development currently 


in use are wind tunnel, free jet and connected pipe (or bell A 


mouth) inlet testing. The free jet method of engine devel- 
opment testing is discussed herein. Performance levels 
obtainable with the components which comprise a free jet 
test arrangement are presented and parameters necessary to 
define test hardware, testability of an engine in a facility and 
optimization of a facility for this type of test are discussed. 


General Considerations 


It is considered essential in the development testing of air- 
breathing powerplants to simulate as closely as possible those 
conditions which would exist in flight. Full-scale develop- 
ment of a complete engine is desirable, if not imperative, be- 
cause of inadequate scaling criteria for the combustion proc- 
esses, structural components and possible coupling effects 
between diffuser and compressor or diffuser and combustor. 

It is desirable, of course, to provide exact flight simulation 
of all internal and external flow conditions. However, this 
can be accomplished only by flight testing or by use of a wind 
tunnel. Development by flight testing can be dangerous 
as well as expensive. Although wind tunnel testing probably 
costs less than flight testing, it is still prohibitively expensive 
for testing all but small powerplants. This can be attributed 
to the large tunnel size and power requirements to provide 
complete flight simulation. Generally, a nonreturn type wind 
tunnel is necessary because the engine exhaust products con- 
taminate the tunnel air, thereby necessitating large-capacity 
driers and heaters for the continual inflow of atmospheric air, 
or an elaborate scavenging system. 

To date, the simplest and most economical engine devel- 
opment testing method available is considered to be the di- 
rect-connect or bell mouth inlet system. This system re- 
quires the least facility power because only that air flow re- 
quired by the engine in flight must be supplied; however, 
no external simulation is provided and internal simulation is 
limited. Such factors as boundary layer shock wave inter- 
actions and their resultant effects upon diffuser velocity 
profiles, angle-of-attack effects on diffuser performance and 


Presented at the ARS Semi-Annual Meeting, San Francisco, 
Calif., June 10-13, 1957. 

1 Project Engineer. 

2 Group Engineer. 


May 1958 


A, (ANNULUS) 
489 


diffuser-combustor (or diffuser-compressor) coupling effects 
‘annot be evaluated. 

A suitable compromise between the aforementioned 
methods can be obtained by use of the free jet test arrange- 
ment where exact flight conditions of Mach number, pres- 
sure, temperature and angle-of-attack are simulated at the 
diffuser inlet. This provides flight simulation of all flow 
parameters inside the engine, provided the engine exit flow 
is sonic. Consequently, all internal components and the 
aerothermodynamic system can be evaluated with much 
less power than a closed test section wind tunnel and a much 
larger degree of flight simulation than that offered by the 
direct-connect type test. 

The expression ‘“‘free jet’’ as used in this article refers to a 
nozzle with no constant area test section downstream as com- 
pared to a wind tunnel. 


Flight Simulation Requirements 
Inlet Flow Field 


"Free jet nozzle configurations which have either rectangu- 
lar or circular exit cross sections are most common. Gener- 
ally, the rectangular end section is used when optical instru- 
mentation is desired since it is readily adaptable for such a 
purpose. In addition, variable Mach number free jets 
are more easily developed in rectangular dimensions. How- 
ever, nacelle type engines are usually circular and, for a given 
inlet diameter, a circular free jet of smaller size than a square 
one can be used, thus reducing facility air flow requirements. 
Also, as noted in (1),* better free jet diffuser performance is 
obtained if the diffusing passage is geometrically similar to 
the entering air flow. Generally, the free jet and diffuser 


parameters discussed herein are for the circular type, al- 
though the same general considerations prevail regardless of 
shape. 

Minimum free jet sizes are desirable in order to minimize 
conversely, to 
The 


facility air supply power requirements or, 
permit tests of larger engines in a particular facility. 


Ay 


Me 


— 


DATA 
< 
fo) 
° 
134 
< 
z 
‘ A,*FREE JET EXIT AREA 
124 
A, "MINIMUM DIFFUSER AREA 
z 
14 


1.0 Ls 20 25 3.0 3.5 4.0 
FREE JET MACH NUMBER-Mo 
Fig. 1 Maximum contraction ratio that permits starting of 
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SUPERCRITICAL SUBCRITICAL minimum size is defined as that which will permit starting 
of supersonic flow into the engine and free jet diffusers, and 
which will contain the engine inlet flow field within the test 
rhombus. 

The requirement for starting supersonic flow into the free 
Ete jet and engine diffusers is shown in Fig. 1. As shown, this 
theoretical curve can be approached very closely in prac- 
tice. The flow area which must start (Ao) is the free jct 
xit area and the contracted area (A;) is the sum of the mini- 
mum areas of the engine and free jet diffusers. 

The requirement for simulation of the inlet flow field after 
starting is accomplished is somewhat more complex since it 
is a function of engine inlet type, engine operating condition 
and free jet nozzle operation. For example, it is postulate: | 
that for proper inlet flow field simulation, the Mach or shoc' 
wave emanating from the free jet nozzle must intersect the 
engine inlet somewhat downstream of the cowl lip and en- 
close the inlet subsonic flow region (if one exists). Howeve 
shock systems for external compression inlets and partici 
larly isentropic spikes, operating below design Mach numbe: 
are sometimes impossible to predict theoretically. Fig. 2 d 
picts typical shock systems which would be encountered wit! 
two different inlet types. These may be obtained experi- 
mentally by model scale tests utilizing conventional optica! 
methods. The definition of minimum size free jet require- 
ments is in many instances made only by analysis of shadow 
graph or schlieren pictures. 

Fig. 3 presents the experimentally determined free jet 


FREE JET 


ISENTROPIC INLET diameter requirements for an isentropic spike inlet designed 
Fig. 2. Flow patterns required for inlet simulation for Mach 2.50. For comparison, the analytically deter- 
+h 5 mined requirements for an isentropic spike inlet or design 
at each Mach number are also included. The calculations 
EXPERIMENTAL ve were based upon design information presented by J. Conners 
REQUIRED TO ENCLOSE SONIC LINE - and R. Meyer (2) with the assumption that the minimum 
—— °° —_ annulus area required for starting occurred at the inlet plane 
4 ee Consequently, it can be seen that the governing requirements 
\ basibeeinca. for above design Mach number are those which can be at- 
SUBCRITICAL CALCULATED tributed to starting supersonic flow and those below design 
\ jp Mach number to the enclosure of the inlet subsonic flow re- 
. S rd : gion in the test rhombus. The latter requirements require 
ws 24 Se quite large free jet nozzles, as shown. 
7 is 20 25 30 35 40 45 5.0 55 
Proper simulation of inlet flow conditions described above 
Fig. 3 Typical minimum ratios of free jet nozzle exit diameter would be of small value in free jet testing unless provisions 
to engine inlet diameter were also made for obtaining exhaust flow simulation. This 
) can be accomplished by choking the engine exit nozzle throat. 


rani sf Hence, with inlet flow established and a sonic condition exist- 
. . . . 
Pexnaust ing at the exit nozzle throat, complete internal flight simu- 
ai lation (with the exception of the divergent portion of the 
Pa exit nozzle, if such exists) is obtained. Therefore, the maxi- 
2 mum altitude that can be simulated (as far as the exit is con- 
. cerned) is a function of the efficiency with which the sonic 
& exit flow is diffused to external cell pressure according to the 
relationship 
= (Pen) (Pa/Pem) 
exh = minimum facility exhaust pressure available 
= P. facility exhaust p labl 
(atmospheric pressure for pure blowdown 
2 operation) 
= engine exit nozzle minimum pressure ratio re- 
z oa quired for choking 
org a ae = minimum total pressure which will maintain a 
choked exit; defines maximum test altitude at 
. 
ad 18 20 a given Mach number 
NOZZLE PRESSURE RATIO Ptg/Pexnaust 
Fig. 4 Exhaust diffuser performance obtainable with engine This relationship reveals the importance of efficient exhaust 
divergent exit flow diffusion (minimum which directly affects 
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OW ] . FREE JET MACH NUMBER, My 
8 | Fig. 6 Minimum operating pressure ratios of circular free jet 
J VA diffusers and nozzles; angle-of-attack = 0 deg; free jet exit 
ne’ to enginefinlet diameter = approx. 1.3 
ter- 
igh 1.0 | 4 
ons 10 Ls 2.0 2.5 3.0 35 4.0 48 
ers MACH NUMBER M, 
uni Fig. 5 Minimum supersonic diffuser operating pressure ratios 
ne. 
nts 
at- maximum test altitude and/or exhauster power require- 
ign ments. 
re- A simple exhaust diffuser such as that formed by a diver- 3 - 
ire gent exit nozzle will yield a good pressure recovery. Fig. 4 = °F P 
shows that a value = 1.2 is obtainable. For simpli- wiTh MODEL 
fication, this is obtained by plotting the engine cold flow (no —_ 
burning) pressure recovery vs. exit nozzle pressure ratio. At & 
values of P,;/P.x, above which the exit is choked, the in- w sob 
ve ternal engine pressure recovery P;5/P1. is fixed as shown by 2 — 
+a the horizontal portion of the curve. The break point in the : ao} 
MS curve defines the minimum value of P,;/P.., which will main- g 
at. tain choked exit conditions. Although these data are for z : 
st- air, data for burning conditions show only slight variation : sof 
from this minimum pressure ratio. 
he 2 
4 Free Jet Diffusers 
| The free jet nozzle which has no physically bounded test 7 TYPE 0 
* section (hence the name free jet) may be operated over a | a acca 
wide range of operating pressure ratios depending upon the 
degree of efficiency with which the free jet kinetic energy is 
converted into pressure energy. : 
. . > 
In a simple stored air blowdown facility, where P..;, equals : 
atmospheric pressure, high efficiencies are mandatory since 
the maximum test altitude and test duration are directly a : 
n hi i } d fo oki diffusers and nozzles; angle-of-attack = 0 deg; free jet width 
much higher t tan the pressure ratios required for cho ing t 1e to engine inlet diameter = approx. 1.3 
: engine exit. This importance is further evidenced in an 
altitude facility where exhausters are required to evacuate 
the test cell. Reference (1) notes that starting pressure jet arrangement will vary between the limit of free jet nozzle 
" ratios for free jet nozzle diffusers may be slightly higher than overexpansion and that of wind tunnel diffusers. The band 
‘ the minimum operating pressure ratio, particularly at higher of optimum free jet diffusers serves as a guide to show the © 
Mach numbers. However, this is relatively unimportant in relative performance level that can be currently achieved. 
all but extremely marginal facilities. Generally, if any range Due to the numerous arrangements that can be evolved with — 
t of testable altitudes is available, sufficient pressure ratio will free jet systems, a complete survey of free jet diffuser mini-_ 
3 exist for starting at the lower altitudes. mum pressure ratio performance would be difficult to pre- 
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Fig. 8 Typical altitude free jet test facility and its test capa- 
bilities; Mach 2.50, ramjet at minimum power 


Since the engine itself forms part of the free jet diffuser 
passage, some compromise has to be made between good 
diffuser performance and engine accessibility and viewability. 
The extent of the compromise is dictated by the free jet 
nozzle size required and the ability of the facility to operate 
these nozzles. 

A comparison of free jet diffuser operating pressure ratios 
is shown in Figs. 6 and 7 for circular and square geometrical 
arrangements. The optimum pressure ratio can be obtained 
from a configuration such as type A, where the diffuser pass- 
age is geometrically similar to the pattern of inflowing air and 
a long constant area section is followed by a gradual sub- 
sonic diffuser. This makes the engine quite inaccessible, 


however, for repair, modification and instrumentation. = 5 
with this type are reported in (1). : 
Two types of free jet diffusers that are used by Marquardt 


Aircraft Co. for full-scale ramjet engine testing are the “‘tai- 
lored diffuser” and the “alternate diffuser.” The performance 
of both are shown in Fig. 6 as types B and C, respectively. 
The tailored diffuser is comprised of a relatively long circular 
diffuser which, together with the engine cowl, forms the 
annular duct passage. Some flexibility and accessibility is 
sacrificed, but good minimum operating pressure ratios are 
obtained. 

The alternate diffuser is a design compromised in favor of 
flexibility. This short segment which is added to the free 
jet nozzle markedly improves overexpanded nozzle perform- 
ance and makes the engine more accessible for inspection, 
repair and modification. In addition, remotely controlled 
transient angle of attack may be accomplished without alter- 
ing the test setup. 

A square free jet diffuser, similar to the circular version 
(type A), performs as shown in Fig. 7, type D. Slightly 
poorer performance is shown (1), probably due to the circular 
engine and square diffuser forming a less efficient diffusing 
passage. 

An arrangement commonly called a “pressure box” system 
which permits shadowgraph or schlieren coverage of the inlet 
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Fig. 9 (a) Supersonic engine total pressure ratios. (b) Maxi- 
mum blowdown test altitudes for supersonic ramjet, turbojet 
and typical free jet diffuser 


flow field is shown in Fig. 7 as type E. Good performane 
without an engine installed has been shown by Ohio State 
University (3), while the University of Minnesota (4) shows 
somewhat poorer performance with a model installed, the 
inlet area of which was 50 per cent that of the free jet nozzle. 
This is to be expected considering the less favorable turning 
angles and diffusion passages associated with a model in- 
stalled in the pressure box. 

It can be seen that various arrangements of free jets and 
diffusers which fall into the optimum free jet diffuser band of 
Fig. 5 5 are possible. The selection is dictated by engine inlet 
size, desired test altitudes and facility cap: abilities. 


Testability 


Testing capabilities of an engine test facility may be most 
conveniently illustrated by presenting a facility operation 
curve. The capacity is set by the operating limit of the 
equipment which makes up the basic testing facility. For 
convenience, the AEDC engine test facility (Tullahoma, 
Tenn.) was chosen as a representative test facility which is 
capable of conducting development testing on a supersonic 
air-breathing powerplant (5). Fig. 8 presents an arrange- 
ment of an altitude free jet test cell and its operational capa- 
bilities. The inlet test air is supplied continuously through 
the supersonic free jet nozzle by a system of compressors, 
where part of this air passes through the engine and part is 
spilled around the engine. In the test cell downstream of 
the engine, the spilled air mixes with the hot gases being 
ejected by the engine, the mixture then being evacuated by a 
system of exhausters. 

The effect of operating pressure ratios upon the range of 
test altitudes available from this particular facility for types 
B and C free jet diffusers at Mach 2.5 is also presented in 
Fig. 8. The free jet size was arbitrarily chosen to be limited 
to a 70,000 ft test altitude with a type B diffuser. The simpler 
type C diffuser does not permit free jet nozzle operation 
at any altitude, whereas the type B diffuser permits opera- 
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tion from 48,000 to 70,000 ft (with a small untestable region 
around 55,000 ft). Required operating pressure ratios are 
shown in Fig. 6. Only zero degrees angle of attack testing 
was considered and flight inlet temperature simulation was 
assumed to vary according to an NACA Standard Day. 
The high altitude limit can be attributed to the maximum 
capability of exhauster equipment to evacuate the total mass 
flow handled to the required cell pressure, while the inlet 
compressor limitations are responsible for the lower altitude 
limit. The determination of these limits is obtained from 
the intersections of the free jet nozzle characteristic lines 
with the facility operational curves. It is evident that 
choice of proper test hardware plays an important role in 
determining the testability of a given supersonic air-breathing 
powerplant in a particular facility. With improved effi- 
cencies in free jet diffusers, exhaust diffusers and optimiza- 
tion of exhauster systems, significant increase in testable 
\\titude can be made. 


Facility Exhauster Optimization 


Altitude free jet test cells have generally been designed to 
operate in the manner described above. Namely, the free 
jet spilled air is mixed with engine exhaust gases and both 
are evacuated to a pressure which is required to provide inlet 
flow simulation (operate the free jet diffuser). However, 
« much higher pressure ratio is generally required to permit 
inlet flow simulation than to provide engine exit simulation 
(choke the exit). Consequently, it follows that a more effi- 
cient utilization of test facility capabilities can be obtained 
hy handling the diffusion of free jet spilled air and engine ex- 
haust products separately. This type of “split cell” arrange- 
ment can be accomplished by having either the engine gases 
exhausting to atmosphere or to a system of exhausters sep- 
arate from those exhausting free jet spilled air. 

For the case where the engine exhausts to atmosphere (and 
assuming that sufficient facility exhaust capacity is available 
for the purpose of exhausting the free jet spilled air), engine 
operating pressure ratios can then be converted directly into 
lines of maximum altitude on a Mach number vs. altitude 
plot. This can be illustrated by considering two typical 
cases of supersonic air-breathing powerplants, the ramjet 
and the turbojet. Fig. 9(a) presents representative pressure 
ratios for each engine as a function of Mach number. The 
conversion of these pressure ratios into altitude limitations is 
also presented in Fig. 9(b) where they are compared with the 
pure blowdown limits for a type B free jet diffuser. As 
shown, much higher altitude simulation at high Mach num- 
bers is possible by arranging the engine exhaust products and 
free jet diffuser spilled air to discharge to separate pressure 
levels. 

For both cases a supersonic pressure recovery of an opti- 
mum axisymmetric isentropic spike inlet was obtained from 
(2), and a pressure ratio Pi5/P.., = 1.2 was assumed for the 
exhaust diffuser. 

For testing at altitudes above those shown in Fig. 9(b), ex- 
hausters must be used to reduce engine exit pressure as well 
as free jet diffuser discharge pressure. Significant gains can 
be made over the upper testable altitude shown in Fig. 8 by 
dividing the exhauster capacity such that a portion exhausts 
the engine gases and the remainder exhausts the free jet 
spilled air. For example, the maximum altitude testable at 
Mach 2.5 in Fig. 8 is 70,000 ft. This is determined by the 
intersection of the combined engine gases plus spilled free 
jet air characteristic with the facility exhauster operational 
curve. However, in Fig. 10 the free jet spilled air and engine 
gas flow characteristics are plotted separately. Conse- 
quently, if the total facility exhaust capacity is split such that 
half exhausts engine gases and the other half spilled air, then 
the maximum testable altitude becomes 90,000 ft. Con- 
versely, for the same altitude limit of 70,000 ft a 40 per cent 
larger frontal area engine could be tested. 
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Fig. 10 Increase in maximum test altitude with split facility 
exhauster arrangement; Mach 2.50, ramjet at minimum power 
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with and without split facility exhauster arrangements; angle-of- 
attack = 0 deg 


Complete optimization of an altitude exhaust facility 
would require very complex ducting to a large number of 
individual exhausters. The exhauster capacity would have 
to be divided between engine and free jet spilled air at differ- 
ent ratios, depending upon Mach number, free jet to inlet 
diameter ratio, and free jet diffuser type. The testable re- 
gion gained for the example set forth would be as shown in 
Fig. 11. Much of this gain can be achieved, however, by pro- 
viding for one splitting arrangement such as half for free jet 
spilled air and half for engine exhaust gases. The gain for 


this arrangement is also shown in Fig. 11. a om 


Facility Nomograph 


Rapid estimation of testing capabilities can be facilitated 
by the use of a nomograph such as that shown in Fig. 12. This 
may be best illustrated by considering a typical example of 
given conditions and solving for the limiting testable alti- 
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Fig. 12 Nomograph for rapid estimation of Mach number- 
altitude test limits 


tudes. The problem can be resolved into two parts, one 
which will be determined by the facility exhauster limits and 
the other by the facility inlet pressure limits. 

For a typical determination of high altitude boundaries 
on a Mach number vs. altitude plot (exhauster limits) enter 
chart I of Fig. 12 with starting points 1, 2, 3, 4 and 5 at Mach 
2.50. Enter charts II and III to respective compression ratio 

and free jet area requirements. The intersection of these 
requirements in chart bid will yield a locus of altitude points: 


relative to the available facility operational curve. The inter- 
section of a line through these points with the facility oper- 
ational curve determines the maximum testable altitude at 
the chosen Mach number. 

Similarly, low altitude boundaries (inlet pressure limits) 
are determined by following the same procedure, except, 
when entering chart II, the P,, line is used instead of com- 
pression ratio. 

For other Mach numbers, the same procedure may be em- 
ployed provided requirements for charts II and III are 
known. 


Conclusions 


Free jet testing is perhaps the most efficient engine devel- 
opment testing method available. The test apparatus, how- 
ever, should be optimized to minimize facility cost or, con- 
versely, to permit testing over a much broader range of flight 
conditions. This optimization consists of defining minimum 
free jet nozzle sizes and providing efficient free jet diffusers. 
Further improvements of large magnitude can be obtained 
by utilizing facility exhausters to evacuate free jet spilled 


air and engine exhaust gases separately. 
2 
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Ground testing of ramjet engines must be accomplished 
under simulated flight conditions. Because the ramjet is 
an air-breathing engine that must operate at supersonic 
air speeds at high altitudes, the requirements for simula- 
tion demand extensive facilities. The Wright Aeronautical 
Division of Curtiss-Wright Corp. operates a blowdown-type 
facility capable of a wide range of simulation for testing 
large full-scale ramjet engines. 


Presented at the ARS Semi-Annual Meeting, San Francisco, 
Calif., June 10-13, 1957. 
1 Manager, Ramjet Test Division. 
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of Flight Simulation for Ramjet Engines 


RAYMOND GREENBERG: | 


Curtiss-Wright Corp., Wright Aeronautical Division, Wood-Ridge, N. J. rah hire 


HE development of ramjet engines as air-breathing pro- 
pulsion devices in the supersonic range fostered the re- 
quirement for testing in ground facilities under flight condi- 
tions. By virtue of its inherent design, the ramjet requires 
speed and altitude before it can function. In flight, this is ac- 
complished with the aid of an auxiliary booster. Conse- 
quently, in the interests of economy, time and reliability, the 
answer to ramjet testing lies in simulating flight conditions in a 
ground facility. 
Creating the environment of the flight regime is in itself a 
technique that requires extensive facilities, equipment and in- 
genuity. It is first necessary to understand the requirements 
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of high altitude supersonic flight so that we may more fully 
appreciate the techniques involved. 


r . . . q 
The ramjet engine is air-breathing by design, and for op- | 
timum performance is sized to operate in the supersonic range | 


at altitudes greater than 40,000 ft (as shown in Fig. 1). Under 
conditions, high stagnation temperatures are en- 
Although operation is virtually all in the isother- 
For ex- 


these 
countered. 
mal range, the temperature variation is considerable. 
ample, at Mach 3 the temperature is approximately 650 F, 
and at Mach 4 it becomes 1300 F. The ramjet is basically a 
simple power plant. It requires air and fuel as do conventional 
engines but is devoid of any major rotating parts. 


precisely designed supersonic inlet which, for practical pur- 
poses, can be considered comparable with a turbojet com- 
pressor. The inlet may or may not be a part of the engine, 
and, in the latter case, the inlet may be designed integrally 
vith the air frame of the flight vehicle. Thrust for propulsion 
s obtained by virtue of the mass flow and the change of ve- 
locity of the gases going through the engine from the inlet to 
the exhaust. The component sections comprising a ramjet 
ean be considered as follows: an inlet diffuser, combustion 
‘hamber, exhaust nozzle and a fuel control system. Fig. 2 
describes a typical ramjet with respect to station locations. 

As previously discussed, the inlet diffuser can be considered 
analogous to a compressor and, in addition, serves to reduce 
the air velocity from supersonic at the inlet entrance to sub- 
sonic at the face of the combustion chamber. With regard to 
ground testing, the inlet simulation poses one of the most ex- 
acting challenges because of the supersonic flow requirements 
and inlet velocity distribution. 

Combustion takes place in the combustion chamber where 
the energy level is raised. This results in the pressure drop 
necessary to expel the gases through the exhaust nozzle where 
they are accelerated to the terminal velocity required to pro- 
duce thrust. 

The fuel control system consists of a control device, fuel 
pump and distribution system. Performance of the engine is 
controlled by the fuel system which requires considerable de- 
velopment under simulated flight conditions. The fuel sys- 
tem in supplying fuel to the engine senses the effect of engine 
combustion on the inlet diffuser and automatically compen- 
sates to maintain the desired fuel-air ratio. Hence, flight 
simulation of these conditions is necessary to study the inlet 
shock-positioning and buzz characteristics of a particular ram- 
jet installation. 

Performance of a ramjet at present is considered optimum 
in a region of approximately Mach 2 to Mach 3.5, and alti- 
tudes of 40,000 to 100,000 ft. The accompanying stagnation 
temperatures at the ram inlet, which are a function of speed 
and altitude, are on the order of 250 to 800 F. The weight 
flow of air, of course, is dependent upon the size of the engine, 
which at the present level of development has become quite 
large. Consequently, the ground facility necessary to provide 
the quantities and temperatures of air and fuel and, further, to 
depress the atmosphere to the required pressure altitudes, has 
become an extremely elaborate and intricate structure. The 
Phase IV facility at the Wood-Ridge plant of the Wright 
Aeronautical Division, Curtiss-Wright Corp. (illustrated in 
Fig. 3), is in altitude facility for accomplishing these aims. 

The facility is basically a blowdown type that supplies air 
to the engine by utilizing high pressure air from storage bot- 
tles, coupled with low pressure air from continuously operated 
compressors. Large capacity heaters raise the air tempera- 
ture to that required to simulate the desired flight condition. 
Pressure altitudes are achieved by exhausting the air through 
a series of steam-fed ejectors. The air supply system (shown 
schematically in Fig. 4) can be operated in either an inter- 
mittent fashion for blowdown, which is necessary for relatively 
low altitude testing in the neighborhood of 45,000 ft, or contin- 
uously from compressors which permit higher altitude opera- 
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Fig. 2 Typical ramjet engine showing component locations 


tion. Air can be supplied to either of two test chambers 
through insulated ductwork and a surge chamber. 

One chamber is 12 ft in diam and 60 ft long with a removable 
surge chamber to permit engine installation. The other is 14 
ft in diam and 90 ft long with 36-ft side hatch for ease of en- 
gine handling. The surge chamber on the latter is also re- 
movable and, in addition, is provided with a nozzle chamber 
that may be included between the test chamber and the surge 
chamber. Both chambers are pressurized by incorporating 
hydraulic locks and water-inflated seals. Three modes of 
testing can be accomplished, namely, direct connect, free jet or 
ducted nozzle. Both chambers are operated from a common 
control room that features a schematic process control panel 
(shown in Fig. 5), an engine control panel, data recording 
equipment and periscopes for observing and recording com- 
bustion chamber flame conditions. 

The primary or blowdown system consists of 240 storage 
bottles that contain 63,000 lb of air at 3200 psi. This air is 
expanded through a control valve which automatically pro- 
grams for a selected surge chamber pressure. Air flows up to 
700 lb per see at 90 psia are obtained in this manner. The 
control system is operated either manually or automatically 
so as to give a 9-sec response with 3 per cent overshoot, main- 
taining +} per cent stability. During this operation, sec- 
ondary air is introduced and mixed with the high pressure air 
through an air-to-air ejector. Consequently, since the ram- 
jet engine provides no motive power of its own, it must be 
boosted to a minimum altitude at a supersonic velocity before 
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Fig. 3 Phase IV ramjet test facility located at the 


ignition can be realized. Utilizing the blowdown air system, 
then, permits high mass flow to simulate the low altitude con- 
dition, and, once ignition is effected, the high pressure system 
is turned off and transition to a higher cruise altitude is main- 
tained by a virtue of the secondary air from the compressor 
which is permitted to continue through the air-to-air ejector. 
This air flow is then controlled by regulating the compressor 
output. 

The secondary air is supplied by turbojet aircraft com- 
pressors that are driven by a 15,000-hp prime mover. Various 
combinations of compressors are available that provide flows 
up to 215 lb per see at 73.5’’ HgA or 130 lb per see at 200’ 
HgA. The compressor output is controlled by regulating the 
prime mover speed and through bleed valves in the com- 
pressor discharge ductwork. A comprehensive chart of 
available air is shown in Fig. 6. 

Air heating is accomplished, in the case of intermittent 
supply, by passing the air through either a horizontal tube and 
shell heat exchanger which is 6 ft in diam by 70 ft long, or a 
metal fill type vertical heat exchanger that is 8 ft in diam by 
80 ft long. The air within the shell is heated by forced air 
from direct-fired air heater that utilizes steam-atomized #2 
fuel oil for combustion. The heat realized from the tube and 
shell heater is 4,400,000 Btu per hour and can raise the air 
temperature to 800 F, while the metal fill heater releases 
7,000,000 Btu per hour and raises the air temperature to 1260 
F. Continuous air is caused to flow through a 94-ft-high 
radiant, convection-type, natural-draft heater. Steam- 
atomized #2 fuel oil is likewise used for fueling the burners. A 
heat release of 299,000,000 Btu per hour is realized in raising 
the air temperature to 800 F. Air temperature is controlled 
by regulating the heater temperature and preheating the duct- 
work to the test chambers to assure stabilization. Air require- 
ments for flight conditions up to Mach 3.25 and altitude range 
of 40,000 to 83,000 ftare achieved inthisway. 
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Fig 7 is a schematic diagram of the exhaust system which 
serves to depress the test cells to the desired pressure altitude 
and utilizes the pumping action of eleven steam ejectors. The 
steam is used as the primary source and is controlled through 

remote-operated valving to each ejector. The secondary 

ports of the ejectors are similarly valved to-the cell exhaust 
manifold. In this way the ejectors, which are sized from 24 to 
48 in. in diam, can be selectively operated to produce test 
chamber pressures in accordance with the curves in Fig. 8. 
A total capacity of 3,400,000 Ib per hour of steam is available 
through these ejectors to handle the 700 lb per sec air flow at 
45,000 ft. Steam for the ejectors is supplied from six package 
boilers which produce 410,000 Ib per hour of steam. The 
boilers charge seven partially water-filled accumulators with 
steam at 500 psi. The total storage capacity of these ac- 
cumulators is 245,000 lb. This steam is released through con- 
trol valves that regulate the pressure to approximately 156 
psi. 

The steam and exhaust products are exhausted to atmos- 
phere through a 57-ft-high silencer fitted with “sound stream” 
sound treatment for noise attenuation. 

Because of the high temperature (approximately 3500 F) 
gases emitted from the engine, an elaborate water cooling 
system is provided. Fig. 9 is a schematic diagram of the water 
system. The temperature at the exhaust gas cooler enters at 
800 F after preliminary cooling in the test cell exhaust diffuser. 
Water sprays and exhaust gas cooler drop the temperature 
to 150 F, after which the water is separated from the exhaust 
gases through a separator section that removes 98 per cent of 
the moisture. This water is regenerated through a 33,000- 
gpm circulating pump that is driven by a 1000-hp synchronous 
motor. The circulating pump returns the water to a 67-ft 
storage tank that has a capacity of 2,100,000 gal of water. 
The warm and contaminated water is returned to the top of 
the tank from where it is drawn off by a 5600-gpm pump that 
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7) Fig. 4 Schematic diagram of the air supply system — 


pumps the water to an air-flotation, oil-separating system. 
The oil and contaminating matter are floated to the surface of 
the separating basin as a result of air entrainment that creates 
bubbles which adhere to the foreign matter. Skimming pad- 
dles that are operated on an endless conveyor belt skim the 
impurities over a weir and the remaining warm water is 
pumped to the top of either of two cooling water towers. The 
towers are equipped with four controllable fans that create 
air columns against which the water rains down. The tem- 
perature is dropped in the cooling towers from 140 to 80 F. 
This cleaned and cooled water is then returned to the bottom 
of the storage tank from which it is gravity-fed to the exhaust 
cooler system. The water regenerative cycle is applied con- 
tinuously during the test operation. 

The fuel requirements for a ramjet are very demanding in so 
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far as mass flow and temperature are concerned. Therefore, 
a fuel system (shown schematically in Fig. 10), is provided 
that has a capacity of pumping 90,000 Ib per hour through 
three centrifugal pumps that raise the pressure to 600 psi. 
The pumps can be operated selectively and can be caused to 
flow through any of the three systems to either test chamber. 
Each of the three systems is sized for a particular range and 
‘an be individually controlled by remote control valves and 
pressure regulators 

High fuel temperatures are necessary to simulate the stag- 
nation temperature effect in flight. This is accomplished by 
passing the fuel through heat exchangers that are heated with 
500-psi steam. Three-way control valves are provided for 
by-passing and mixing hot and cold fuel and thereby regulat- 
ing the temperature to a maximum of 350 F. Where addi- 
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tional temperature is required, the 350 F fuel is passed through 
a chemical heat exchanger that raises the temperature to 500 F. 

Testing basically consists of investigating performance 
characteristics and durability of each system of the engine as 
well as determining acceptance of the complete engine. To 
accomplish these aims, conditions must be provided to ac- 
curately simulate the in-flight environment. Perhaps the 
most exacting are the requirements for inlet simulation. Be- 
cause of unpredictable flight paths, studies of engine operation 
must be made both at design and off-design conditions. 
Hence the problem of inlet simulation. 

The effect of off-design operation is a variation in the inlet 
air distribution at the engine face or combustion chamber en- 
trance. Some characteristic profiles plotted for Mach num- 
ber distribution in the subcritical region under varying condi- 
tions of speed and angles of attack are shown in Fig. 11. Of 
significance is the location of the high Mach number region 
observed in the center contours. The shifting of this high 
Mach number area, as noted in the diagrams, is reflected in 
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flame propagation in the combustion chamber. This may 
have a detrimental effect upon engine performance or dura- 
bility. It is also interesting to note the contour rings which 
represent the low Mach number region and the positions they 
assume when the high Mach number area changes. Under ad- 
verse conditions, flow separation takes place in the peri- 
pheral contours and can seriously affect the boundary layer 
used for cooling air flow. The penalty for starving the 
cooling air obviously is lack of durability in the combustion 
chamber and exhaust nozzle. 

For the particular cases illustrated, a side inlet was used but 
the problem is similar regardless of the inlet configuration. 
This study is of vital importance to the development engineer 
and, therefore, ground simulation of the in-flight inlet condi- 
tions must be conducted. A ducted nozzle is employed which 
provides the means to vary the inlet conditions necessary to 
affect the simulation of flight Mach numbers and angles of 
attack. 

Inlet buzz phenomena studies are similarly important in the 
development of a ramjet—particularly as they affect the fuel 
control system. Buzz can be described as an objectionable 
pulsation of the air flow through the engine. It is controlled 
by positioning the oblique and normal shock waves in the inlet 
section indicated by the dotted lines on Fig. 2. The oblique 
shock is induced by the inlet geometry and the normal wave is 
a function of the combustion chamber pressure. The oblique 
wave emanates at the tip of the inlet center body and, in a 
properly designed and operated engine, is captured by the 
outer lip of the inlet diffuser cowl. The normal wave, which 
is normal to the direction of flow, is located in the inlet diffuser 
and, since its position is a function of the thermal blockage 
from the combustion chamber, it may move fore and aft as a 
result of fuel flow changes. Consequently, instability of the 
normal shock position that causes alternate spill and capture 
cycles of the oblique shock results in buzz. This can be 
likened to compressor surge, as it imposes a limitation on the 
fuel control system operating range. 

Free jet testing is conducted to study the effects of the buzz 
cycle and determine the control system margin of perform- 
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ance. The testing technique involves using a supersonic 
nozzle in front of the inlet in conjunction with an exit diffuser 
to establish supersonic air flow simulating the in-flight condi- 
tion. The positioning of the inlet cone with regard to the 
nozzle exit or test section is of extreme importance. This is 
significant in view of the fact that the region of the desired 
Mach number falls within an envelope commonly referred to 
as the test rhombus. The rhombus is formed by the result- 
ant supersonic shock waves that emanate downstream of the 
nozzle throat as the flow expands to supersonic velocities. 
The shocks are confined within this diamond arrangement and 
it is, therefore, referred to as the test rhombus. In perform- 
ing the test, oblique shock wave of the inlet must be located 
within this area to achieve accurate Mach number simulation. 
The positioning of the inlet within the test rhombus and the 
buzz cycle effects of the oblique shock wave are viewed with a 
shadowgraph or schlieren optical viewing system. The re- 
sults are apparent on;a screen surface that may be photo- 
graphed for further analysis and record purposes. 

For small-diameter engines having symmetrical, round in- 
lets, relatively small nozzles are provided. However, for large 
engines, particularly where side inlets are concerned, much 
larger nozzles are required. A variable Mach number super- 
sonic nozzle of the latter type permitting approximately 50 
per cent secondary air flow is shown in Fig. 12. This nozzle 
incorporates two contoured blocks, the lower of which is 
stationary.- The upper, movable block is remotely controlled 
and is operated by a hydraulic ram to produce variations that 
range from Mach 1.5 to Mach 3.25. The nozzle may be op- 
erated with fixed-geometry for on-design conditions and, fur- 
ther, has provisions for rapidly moving the sliding block to 
simulate the booster phase. 


Fig. 12 Nozzle chamber for cell No. 2 containing the variable 


Mach number supersonic nozzle 


The three modes of testing, namely, direct connect, free jet 
and ducted nozzle, provide a broad variety of test simulation 
that covers all phases of testing necessary for ramjet develop- 
ment. Each has a particular advantage that is moderated by 
its limitations. Direct-connect testing is adequate for com- 
ponent evaluation and utilizes relatively simple equipment. 
However, although inlet velocity profiles may be. varied with 
flow simulators, the actual recovery conditions in‘ front of the 
engine inlet are not faithful reproductions because of the sub- 
sonic velocities in the straight inlet duct. 

Free jet testing provides good inlet simulation; however, 
it imposes external drag forces on the engine that are not 
encountered in flight, as well as severe facility limitations 
because of the high pressure drops necessary to create super- 
sonic flow. The high rate of flow necessary for free jet test- 
ing also requires prohibitive exhausting capacity. 

The ducted nozzle provides a happy medium in that the in- 
let simulation for both on-design and off-design testing is 
feasible without taxing facility capacity. However, the 
equipment involved is rather complicated and requires a good 
deal more maintenance and service. Flight and wind tunnel 
tests provide the necessary reference profile data from which 
the ducted nozzle is oriented. Full-scale tests utilizing the 
ducted nozzle have demonstrated extremely close simulation 
of these supplied in-flight data. 

It is apparent from this discussion that we have onl¥:begun 
in so far as flight simulation is concerned. Development of 
flight vehicles and powerplants is rapidly accelerating, and 
the requirements for people and technology to keep pace with 
the growth of the field of aviation offers a rigorous challenge 
for the future. 
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\ method of planning adequate test facilities to meet 
future ramjet test requirements is discussed. Data are 
presented to show that adequate planning with sufficient 
lead time can be accomplished most readily through care- 
ful review and analysis of ramjet component research re- 
sults. The analysis shows that specification of ramjet 
test facility equipment requires careful consideration of 
the over-all operating envelope of ramjet engines since 
multiple use of air supply and exhauster equipment is 
possible. Potential effects of research in the field of test 
simulation methods upon facility performance require- 


ments are discussed. 


Introduction 


| dpe ramjet test facilities are needed because the 

ramjet offers the most desirable aircraft propulsion per- 
formance over a broad range of flight conditions, and because 
methods of extrapolating performance and operational data 
procured from small-scale ramjet engines and components 
have not been evolved. Limited success has been achieved in 
scaling combustion chambers when fuel systems are properly 
scaled but the methods for scaling fuel systems are empirical 
and not completely reliable. 

A method of planning adequate test facilities to meet future 
ramjet test needs is discussed in this paper. Data are pre- 
sented to show that adequate planning with sufficient lead 
time can be accomplished most readily through careful re- 
view and analysis of ramjet component research results. 
The analysis also shows that specification of equipment for a 
ramjet test facility requires careful consideration of the 
over-all operating envelope of ramjets tested in that facility 


since multiple use of air supply and exhs auster F equipment is 


possible. 
Discussion 

A ramjet facility in which air is not recirculated is con- 
sidered in this analysis. Inlet air is taken from atmosphere, 
raised to the proper pressure level in air supply compressors, 
and temperature conditioned by the use of heaters or coolers. 
Air then passes through the test installation, is cooled, and 
then exhausted to atmosphere through exhaust gas com- 
pressors. Effects of water spray cooling on exhaust gas 
volume have been neglected. 

Operating envelopes to be simulated for ramjet engines were 
generated so that test facility performance required at each 
engine operating condition could be determined. Determi- 
nation of test facility perform: ance at each engine operating 
condition was necessary in order to show maximum perform- 
ance limits for air heaters, air supply compressors and ex- 
haust gas compressors. In addition, study of test facility 
equipment performance requirements over the complete 
range of operation for the engines provided information on 
the possibility of multiple use of air supply machines as ex- 
hausters at operating conditions at which exhaust require- 
ments are at a maximum and air supply requirements are 
relatively low. 

Presented at the ARS Semi-Annual Meeting, San Francisco, 
Calif., June 10-13, 1957. 
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Ramjet Operating Envelope 


A maximum Mach number (C-D) was chosen for the 
operating envelope shown in Fig. 1 on the basis of relative 
performance of ramjets, turbojets and rockets. A com- 
parison of powerplants was made by the use of specific im- 
pulse, thrust per unit of frontal area and thrust per unit of 
weight. Fig. 2 (1)’ shows specific impulse as a function of 
flight Mach number for the rocket, ramjet, turbojet and 
turbojet with afterburner. Specific impulse of the ramjet 
decreases above Mach 3 as shown in Fig. 6, but it appears that 
the ramjet offers attractive values of specific impulse at 
Mach 4. Thrust per unit of frontal area of a propulsion 
system is indicative of the drag of the powerplant installation 
and also indicates the variation of weight of a specific power- 
plant. Fig. 3 shows that on a thrust per unit frontal area 
basis the rocket will be superior to the ramjet at Mach 4. 
The rocket is lighter than the ramjet at the maximum Mach 
number studied. 

A maximum sea level Mach number (B) was determined on 
the basis of the fact that it is possible and beneficial to ac- 
celerate ramjet-powered missiles to Mach numbers in the 
range of 2.5 at very low altitudes with the use of rocket en- 
gines. Maximum Mach number operation was limited to 
altitudes above 50,000 ft because dynamic pressure at that 
flight condition is approximately equal to dynamic pressure 
at Mach 2.5 at sea level. 

Maximum Mach number for which the test facility should 
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Fig. 1 Typical operating envelope for use in generating basic 
machinery requirements for ramjet test facility 
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Fig. 2 Typical specific impulse of various types of engines as a 
function of Mach number 


be designed may increase due to advances in research on 
ramjet components and fuels or may decrease if strides are 
made in advancing rocket engine performance. The ramjet 
test facility planner must therefore carefully analyze re- 
search advances in ramjets and rockets as well as the relative 
complexity of development problems in the two types of 
engines in order to estimate the maximum Mach number for 
which ramjet test facilities should be designed. 

A Mach number of 1.5 at altitudes between sea level and 
about 55,000 ft (A-Z) comprises one limit of the operating 
envelope shown in Fig. 1. It is doubtful that ramjet engine 
free jet testing can be accomplished satisfactorily below 
Mach numbers of about 1.5 unless free jet testing philosophy 
is greatly changed. Below M = 1.5, the allowable over-all 
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Fig. 4 Ramjet engine diffuser pressure recovery as a function 
of Mach number 
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Fig. 3 Typical values of thrust/frontal area for various types of 
engines as a function of Mach number 
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cross-section area at the engine inlet cowl is below 50 per 
cent of the Laval nozzle exit area. Choice of minimum Mac) 
number was based on the practice in some ramjet engines of 
igniting the engine at low altitude, low Mach number condi- 
tions to assure ignition and to avoid the necessity of speci:! 
geometry changes in the engine to make ignition possible at 
high flight Mach numbers. 

Improvements in free jet test technique or the use of smaller 
engine inlet to free jet nozzle outlet area ratios to accommodate 
subcritical free jet testing would lower the minimum test 
Mach number. The engine test facility planner should, there- 
fore, carefully follow the research results on such test tech- 
niques. 

Maximum altitude to be simulated at each flight Mach 
number (D-£) was based on maintaining reasonable com- 
bustion chamber pressures. Experimental investigations of 
effects of combustor pressure on combustion efficiency indi- 
cate that minimum ramjet combustor pressures should not 
go below a critical value. One-third of an atmosphere was 
used as the minimum combustion pressure in this analysis. 
Total pressure recovery for a series of inlets, each designed for 
optimum recovery at the Mach number presented, is indicated 
in Fig. 4. These data were used along with the minimum al- 
lowable combustor pressure to calculate the maximum alti- 
tude operating limits in the following manner: 

At a chosen Mach number, diffuser pressure recovery was 
read from Fig. 4. Minimum allowable combustor pressure 
was divided by the diffuser recovery to determine engine inlet 
total pressure. 


= ambient total pressure, psfa 


H» = combustor inlet total pressure, psfa 
= diffuser pressure recovery 


Pressure ratio for the chosen Mach number was read from 
Mach number function tables (Reference 2) and ambient 
altitude pressure was calculated by use of the equation 


| 


where 
po = ambient static pressure, psfa 


Altitude was determined by interpolation using a table of 
NACA standard properties of the atmosphere (2). 

Maximum altitudes which should be simulated in the alti- 
tude facility at various Mach numbers were thus determined 
on this basis for the range of Mach numbers considered. 
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where 


Ramjet research which results in increased diffuser pressure 
recovery or lower minimum combustor pressures would raise’ 
the maximum operating altitude of the ramjet. 


Test Methods 


Fig. 5 shows three methods in use for engine development 
testing. Direct-connect burner testing is used for prelimi- 
nary phases of engine development and produces the least: 
severe test facility requirements. Ducted-nozzle testing with 
choked inlet is a new and not yet universally accepted test 
method. Research is yet to be accomplished before it can 
produce completely acceptable test results for supercritical 
engine operating conditions. Free jet testing produces satis- 
factory inlet simulation and allows engine control system 
testing and diffuser-combustor matching studies. Each engine 
inlet simulation method may be used with sufficient pressure 
ratio across the engine exhaust nozzle to establish critical 
flow at the minimum area section of the nozzle, or ambient. 
altitude pressure may be maintained at the engine outlet. 
\mbient-altitude exhaust pressure permits tests with vari- 
uble geometry and fully expanded exhaust nozzles with 
ejectors. Ducted-nozzle and free jet test methods with 
choked engine exhaust nozzle and ambient altitude exhaust 
pressures were analyzed. 


Facility Performance 


Air flows presented in this paper are for the purpose of 
comparing facility requirements for different types of testing 
and are for an engine having a combustion chamber cross- 
sectional area of 1 sq ft. Prediction of changes in ramjet 
engine size, and therefore test facility capacity, requires study 
of ramjet component research results. At a fixed flight — 
Mach number and L/D ratio for a missile, engine size would | 
increase with increased altitude. As combustion and dif- — 
fuser research produce the capability for high altitude oper- 
ation of ramjets, engine size will increase. 

Static pressure and temperature for operating conditions 
shown in the typical operating envelope in Fig. 1 were taken 
from NACA Standard Altitude Tables (2). Ambient total 
pressures and temperatures at each operating condition were 
calculated from the static values and the Mach number- 
pressure ratio-temperature ratio relations in (2). Engine 
air flow at each operating condition was determined in the 
following manner: 

1 A burner inlet Mach number of 0.2 was assumed. This 
combustor inlet Mach number was assumed based on the 
fact that momentum pressure losses in a ramjet burner in- 
crease as burner inlet Mach number increases, and because 
it becomes very difficult to attain reasonable combustion 
efficiencies at high altitudes and high burner inlet Mach 
numbers. 

2 Burner inlet total pressure was determined by multi- 
plying ambient total pressure for each test condition by the 
diffuser pressure recovery shown in Fig. 4. 

3 Burner inlet total temperature was assumed to be 
equal to ambient total temperature. 

4 Burner inlet static pressure and temperatures were 
determined from the total pressure and temperature and 
Mach number-pressure ratio-temperature ratio relations 
presented in (2). 

5 Air flow was calculated by use of the equation 


Po static pressure at combustor inlet, psfa 
Az = combustor area, ft?, assumed to be 1 in this study — 
M, = combustor inlet Mach number, 0.2 
Y = ratio of specific heats, 1.4 

g gravitational constant, 32.2 fps? 


May 1958 


\ \ 


(SUBSOMIC INLET FLOW WITHOUT ATION) 


. BION 


DUCTED NOZZLE TEST 


(SONIC INLET FLOW WITH PARTIALLY CORRECT ATION) 


“SECOND 


Fig. 5 Free jet test—supersonic inlet flow with correct simula- 
tion 
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Fig. 6 Air flows and supply pressures required to simulate test 
conditions outlined by typical operating envelope ramjet, ducted 


nozzle test 
R= gas constant for air, 53.3 ft lb/°R << 
t, = static temperature at combustor inlet, °R 


Test chamber inlet total pressure as a function of air flow 
for ramjet ducted nozzle testing is shown in Fig. 6. Ambient 
total pressure was presented because testing with a ducted 
nozzle results in losses in the ducted nozzle which are equiva- 
lent to those which would occur in the ramjet diffuser. Let- 
ters shown on the curve in Fig. 6 correspond to those shown 
in Fig. 1 and permit correlation of the air flow and pressure 
requirements with Mach number and altitude. Maximum 
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Fig. 7 Air flows and air supply temperatures required to simu- 
late test conditions outlined by typical operating envelope ramjet, 
ducted nozzle test 


Pe, PSFA 


| 


TEST CHAMBER OUTLET PRESSURE, 


AIRFLOW, Wy, LB/SEC./FT 


Fig. 8 Air flows and exhaust pressures required to simulate test 
conditions outlined by typical operating envelope—ramjet with 
ambient altitude exhaust pressure, ducted nozzle test 


air flow required is about 142 lb/sec per sq ft of engine cross 
section area at a pressure of 36,200 psfa. 

Test chamber inlet temperature is shown as a function of 
air flow for ducted-nozzle testing of the ramjet in Fig. 7. 
Maximum temperature required is about 1650 R at an air 
flow of 46 lb/sec. Temperature required at the maximum 
air flow condition of 142 lb/sec is 1165 R. 

Test chamber exhaust pressure as a function of engine air 
flow for ducted-nozzle testing is presented in Fig. 8. This 
figure presents ambient altitude exhaust pressure. Ambient 
altitude pressure at the highest altitude considered is about 
15 psfa, and engine air flow for operation at that altitude is 
2.8 lb/see per sq ft of engine cross-section area. 

In order to determine exhaust pressure required to main- 
tain critical pressure ratio across the ramjet engine exhaust 
nozzle, the total pressure losses through the ramjet down- 
stream of the inlet diffuser were determined as follows: 

1 A representative value of pressure drop for gutter flame- 
holder of 2¢ was assumed. 
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_ 2 A burner outlet temperature of 4000 R was assumed for 
all cases except those in which the engine outlet Mach num- 
ber reached a value of 1 with no contraction in the exhaust 
nozzle. For these cases, temperature ratio required to pro- 
duce choking at the burner outlet with no contraction was used. 

3 A total pressure loss of 2 per cent due to losses in the 
engine exhaust nozzle was assumed. 

4 Above values of temperature ratio were used with the 
assumed combustor inlet Mach number to calculate momen- 
tum pressure loss across the ramjet combustor. 

Total pressure ratio H;/H» across the ramjet burner in- 
cluding all of the losses discussed above are shown as a func- 
tion of flight Mach number in Fig. 9. These burner pressure 
ratios were used with values of engine diffuser exit pressure 
to determine engine outlet total pressure at each flight con- 
dition on the envelope shown in Fig. 1. Test chamber outlet 
pressure required to maintain critical pressure ratio across 
the engine exhaust nozzle was determined by dividing engine 
outlet total pressure by two. Test chamber outlet pressure 
required to maintain critical pressure ratio across the engine 
outlet is shown as a function of air flow in Fig. 10. Maxi- 
mum air flow at which exhausters are required is about 40 
lb/see per sq ft of engine cross section. Minimum exhaust 
pressure required is about 250 psfa at an air flow of about 
5.2 lb/see per sq ft of engine cross section. 

Air supply requirements are the same for ducted-nozzle 
testing with ambient altitude exhaust and with choked en- 
gine outlet exhaust pressures. At an air flow of 2.8 lb/sec 
the exhaust pressure is 15 psfa for ambient altitude exhaust 
and 295 psfa for choked engine outlet testing. These two 
pressure levels indicate the difference in exhauster pressure 
ratio required for the two types of testing. Maximum vol- 
ume flow through the exhausters for the choked engine out- 
let testing occurs at the minimum Mach number high alti- 
tude operating condition (#) rather than the maximum 
Mach number high altitude condition and is about !/, as large 
as for ambient altitude exhaust. 
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Fig. 9 Internal pressure losses in a ramjet exclusive of inlet 
diffuser losses 
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conditions outlined by the typical operating envelope—ramjet ra 
with choked exhaust nozzle, ducted nozzle test ’ 4 4 
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Free jet test facility requirements were based on the enve 3 


lope shown in Fig. 1. Engine pressure losses and air flow 
requirements discussed for the ducted-nozzle test in this re- 
port were used in the free jet study. 

Inlet total pressure and temperature required for free 
jet testing a ramjet engine throughout the operating envelope 
are the same as those required for ducted-nozzle tests. Air 
flow requirements for free jet testing the ramjet were deter- 
mined in the following manner: 

1 Engine air flow was taken from the ducted-nozzle 
calculations for each of the engine operating points. Free 
stream area through which this air flow would pass was ob- 
tained from 


Wie 
Ap = 
where 
Ay = free stream or cowl lip area, ft? 
W., = engine air flow, lb/sec 
p) = ambient density at the altitude at the operating condi- 
tion chosen, lb/ft? 
Vo = velocity corresponding to the Mach number and alti- 


tude at the operating condition chosen, fps 


2 Free stream or cowl inlet area as a function of Mach 
number is shown in Fig. 11. Cowl inlet area of the ramjet 
engine inlet is larger than the combustion chamber at Mach 
numbers above about 3.8 for the combustor Mach number 
and diffuser pressure recovery assumed in the calculations. 
An inlet area equal to the combustor area was used for deter- 
mination of free jet testing requirements at Mach numbers 
above 3.8. A lower combustor inlet Mach number would, 
because of this assumption, exist in the ramjet at flight Mach 
numbers above 3.8. 

3 Free jet nozzle exit dimension was then determined by 
using a value of about 50 per cent spill air flow for the high 
Mach number free jet test installation. Free jet tests have 
been run at high Mach numbers with spillage air flow as low 
as 35 per cent but these were tests in which the free jet nozzle 
exit was shaped to match the engine inlet, thereby precluding 
use of a variable Mach number free jet nozzle. Also, use of 
the figure of about 50 per cent spillage is justified by the need 
of tests at angle of attack. Justification of a much larger 
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Fig. 11 Ramjet cowl inlet lip area as a function of flight Mach 
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Fig. 12. Air flows and air supply pressures required to simulate 
test conditions outlined by typical operating envelope with free 
jet nozzle 


spillage could have been made on the basis of the requirement 
for subcritical engine and control system testing. The maxi- 
mum engine inlet area was doubled to arrive at a free jet 
nozzle exit area and further calculations are based on oper- 
ating a constant outlet area-variable Mach number free jet 
nozzle over the range of conditions shown in Fig. 1. It was 
concluded that in the event it is necessary to test engines 
having asymmetric inlets in an altitude facility, a free jet 
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nozzle having the correct outlet area would be built with the lb/sec, inlet pressure as low as 710 psfa is required. These 
proper exit dimensions for purposes of economy. data show that use of a given size free jet nozzle which would 
4 Air flow required for operation of the free jet nozzle permit tests of high Mach number engines over the complet: 


range in air supply capacity. For instance, point B on the 
- - operating envelope represents Mach number 2.5 at sea level. 


was determined by use of the equation range of Mach numbers results in a requirement for a broader 
4 
The curve in Fig. 11 showed that at this lower Mach number 


Wa = AnpoMy «| 


Rt 
h the design engine inlet area would be about 65 per cent of the 
i 18 ” le exi ft? Maximum value. Stated in other words, if the engine in- 
“aa = Seem, & let, and the free jet nozzle, were to be tailored for operation 


Hi - a ambient static pressure, psfa at each Mach number, the maximum facility air flow would 
reltien.--eandnsgraad a. be 268 lb/sec for each square foot of engine cross section in- 
& = ambient static temperature, °R stead of 412 lb/sec. Comparison of maximum air flow for fre« 
jet testing on this figure with maximum air flows for ducted- 
Altitude facility air supply pressure is shown as a function nozzle testing previously presented shows that the free jet 
of air flow for free jet testing in Fig. 12. A maximum air air flow is about three times the ducted-nozzle air flow for » 
flow of about 412 lb/sec per sq ft of engine cross section at a given size engine. 
pressure of 36,200 psfa is required. At an air flow of 26.5 Altitude facility air supply temperature is presented as a 
function of air flow for free jet testing in Fig. 13. Maximum 
temperature is required at an air flow of 89 lb/sec per sq ft of 


and other symbols are as previously presented. 


= engine area. Temperature required at maximum air flow oi 
on 1600 . 412 lb/sec is about 1165 R. The fact that maximum temper- 
2 \ ature is required only at air flows less than maximum shows 
” the possibility of economizing on air heater capacity. 
ui 1400 N Test chamber exhaust pressures required for altitude fa- 
- \ N\ cility free jet testing with ambient-altitude static pressure 
« downstream of the free jet nozzle exit and the engine outlet 
_ oe \ are shown as a function of facility air flow in Fig. 14. Mini- 
= | s mum test chamber exit pressure required for ambient altitude 
+ free jet testing is the same as the minimum pressure required 
2 | for ambient altitude ducted-nozzle testing. Maximum air 
= flow for free jet testing is about three times that required for 
ra ducted-nozzle testing. 
2 \ ¥ Severity of the facility equipment requirements to main- 
om \ rs tain ambient static pressures in the test cell will probably 
E 
AIRFLOW, Wa ,LB/SEC/FT* 
Fig. 13. Air flow and air supply temperatures required to simu- 
late test conditions outlined by the typical operating envelope 
with free jet 
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cause economic considerations to dictate use of test installa- 
tion recovery devices. 

Test chamber exit pressures for free jet test installations 
using second throat apparatus were calculated by use of the 
data presented in Fig. 15, and the free jet nozzle exit total 
pressures. Data presented in Fig. 15 were obtained experi- 
mentally with a variable Mach number free jet test installa- 
tion. Exhaust pressure required for satisfactory free jet 
installation operation as a function of air flow is presented in 
lig. 16. Minimum exhaust pressure required is 90 psfa com- 
pared with 15 psfa for ambient altitude free jet testing and 
250 psfa with choked engine outlet ducted-nozzle testing. 
Air flow at the minimum exhaust pressure condition for free 
jet testing with both ambient exhaust pressure and second 
throat recoveries is 5.2 lb/see compared to 3.04 lb/sec for 
ambient altitude ducted-nozzle testing, and 5.1 lb/see for 
choked engine outlet ducted-nozzle testing. Taking choked 
engine outlet ducted-nozzle exhauster volume requirements 
us one, then volume requirements for ducted-nozzle ambient 
altitude testing would be about 9, for free jet testing with 

‘cond throats about 4.5, and free jet ambient altitude test- 
ig about 17. 

It is clear that the cost of facility equipment must be bal- 
need against the relative advantages of ambient exhaust 
ressure testing. Factors such as greater facility utilization 
nd decreased complexity and cost of test installations may 
ompensate for the very large exhauster systems involved. 


Summary 


From this analysis, it was concluded that the successful 
ramjet facility planner should continuously monitor ramjet 
research and development progress and try to anticipate the 
direction and magnitude of that progress so as to satisfactorily 
anticipate requirements for ramjet test facilities. Improve- 
ments in performance of ramjet inlets, combustors, fuels, 
tailpipes, as well as structural design, can produce major 
changes in the range of operation over which ramjets offer 
the best over-all performance potential. Research on ducted- 


104 
A 
LZ) 

< 44 
103 
a 
x 
a” 
A 
a 
2 
© 
102 al 
> 
< 
x 
| 

if 

100 1000 


AIRFLOW, Wy, LB/SEC/FT? 


Fig. 16 Air flows and exhaust pressures required to simulate 
test conditions outlined by the typical operating envelope with 
free jet nozzle and second throat 


nozzle testing and modified free jet testing may show that 
ramjet test facilities would require air supply and exhauster 
equipment greatly different from that currently considered | 
necessary for such facilities. aa 
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This paper contains a formulation of the method of 
characteristics for computing the flow field in a supersonic 


jet. The method is programmed on an ERA-1103 com- 
puter. Nine of the cases calculated are reported in this 
paper. These cases represent jets issuing from two types 


of nozzles at various Mach numbers into a still or moving 
medium at a wide range of pressures. Different values of 
the ratio of specific heats of the gas have been assumed. 
Results include the effects of various parameters on the 
jet-spreading characteristics. 
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Spreading of Supersonic Jets From hil 
Symmetric Nozzles 


and J. B. PETERSON’ _« 


p= 
Nomenclature 
b = defined by Equation [3c], boundary point 
r = radial distance 7 
A = defined by Equation [8a] 
= defined by Equation [3b] 
C = defined by Equation [3d] 7 
M = Mach number 7 
P = pressure 
R = radial distance 
S = entropy divided by specific heat at constant pressure 
U_ = defined by Equation [4a] 
V = defined by Equation [4b] 
W = velocity divided by maximum velocity; subscript 1 is 
before shock, subscript 2 is after shock 
X = axial distance 
a = initial condition points at nozzle exit 
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8 = points on the singularity of the Prandtl-Meyer expansion 

y = ratio of specific heats 

e = nozzle area ratio 

¢ = defined by Equation [7d] ; 

= defined by Equation [7c] 

6 = flow angle measured from axis 

62 = defined by Equation [4d] 4 ag 

x = defined by Equation [7b] il 

= defined by Equation [7a] i 4 

= Mach angle 

vy = Prandtl-Meyer angle 

@ = angle of shock with respect to chu a flow velocity 

Subscripts 

e = exit 

b = boundary q 

c = chamber + 

t = stagnation 

= free stream + 

0 = refers to nozzle exit radius 

RM = point where the boundary has a maximum value of 
radius 

6/2 = point on boundary where @ is equal to one half its initial 


boundary point value 
Introduction 


NOWLEDGE of the flow field of an axisymmetric super- 
sonic gas jet has recently become important in the con- 
sideration of some jet propulsion problems. These problems 
include cooling the structures adjacent to or within the jet, 
and transmitting communication signals through the jet. 

While it is difficult to obtain experimental data of jet 
spreading under high altitude pressures, the method of charac- 
teristics presents itself as a convenient tool for the theoretical 
calculation of the jet properties at any pressure. Further- 
more, the effect of a supersonic external stream on the jet 
characteristics can be included in the calculation by properly 
determining the pressure variation on the jet boundary under 
interaction with the external stream. 

For a gas jet the flow conditions at the nozzle exit are de- 
termined by the nozzle geometry and can be considered as 
initial conditions which are already known and from which the 
calculation of the jet can be started. The necessary boundary 
conditions are provided by the assumption of axisymmetry 
and the fact that the boundary is under pressure equilibrium 
with the surrounding medium. The only gas property needed 
in the calculation is the ratio of specific heats, gamma. In the 
usual temperature range of most supersonic jets, gamma can 
be assumed constant. With the initial conditions, the bound- 
ary conditions and gamma all given, the method of character- 
istics can be readily applied to calculate the jet flow field until 
a shock is encountered. The shock is to be handled by using 
a separate set of equations. 

The finite-difference characteristic equations, as well as the 
machine programming and the equations for the shock wave, 
are given in this paper. The method of characteritics is pro- 
grammed on an ERA-1103 computer. The computation can 
be carried out rapidly. The cases calculated include both the 
conical and the bell-shaped nozzle types; four exit Mach 
numbers, namely 2, 3, 4 and 5; three gamma values, namely 
1.15, 1.25 and 1.35; a range of exit pressure to ambient pres- 
sure ratios from 2.8 to 6000; and three external stream Mach 
numbers, namely 4, 5 and 6. In all the cases the calculation 
was stopped before the shock caused difficulties in the isen- 
tropic computation. This was done because from considera- 
tion of a number of engineering problems it was concluded 
that the jet field would be of practical interest only in the 
vicinity of the nozzle exit. The calculation is carried out in 
most cases to at least 10 diam downstream. 

It should be pointed out that the method of characteristics 
has its limitations. First, it does not consider the mixing 
between the jet and its surrounding medium. Second, the 
method of characteristics is based on a continuum theory. 


~ 8 Numbers in parentheses indicate References at end of paper. 
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The book by. Pai (1) contains extensive references on 
supersonic jets. 


Mathematical Formulation 


Flow Variation Along Characteristic Lines 


The supersonic flow equation, being of the hyperbolic type, 
has two families of real characteristics. Characteristics of the 
first family (left running) form an angle (6 + yw) with th 
positive z-axis, where 6 and yu are the local velocity angle anc 
Mach angle, respectively. Characteristics of the second 
family (right running), on the other hand, form an angk 
(6 — yw) with the positive z-axis. The characteristics have th 
property that the variations of flow properties along thes 
lines are described by differential equations of the first orde 
(2, Equations 512 and 513). These equations are given below 
in the difference form. 

Frrst FaMILy 


AQ = AAM — BAX + CAS.............[1 
Seconp FamMILy 

A@ = —AAM + bAX CAS............ [2 
The various coefficients are defined as 


/M?—1 


— 1 cot — 1)r 
1 
(+/M? — 1 cot 0+ 1)r 
VM -1 [3d] 


The notation AS is dimensionless entropy change, and is 
defined by 
entropy change per lb m 


S = 
C, per lb m 


When r = 0, it can be shown (3, p. 680) that Equation [1] 
takes the form 
@ = }(AAM + CAS).............. [la] 


Flow Change Across Shock 


Let subscript 1 denote the flow immediately upstream of 
the shock and subscript 2 denote the flow immediately down- 
stream. If ¢ is the angle of shock with respect to the upstream 
flow velocity, the following equations can be used to calculate 
the downstream flow 


V, = W, sing 1) eset +1]... 
y+ 1 
[4e] 
6. = 6, + — are tan [4d] 
/2 


The entropy increase across shock is given by 


In [2 + (y — 1)((M sin $)?] wee [2y sin — 


re 
- 
: 
? 
; 
; 
= 
| 
> 
* 
. 


3b) 


5] 


Flow Around Nozzle Edge 


The gas expands around the nozzle edge when leaving an 
underexpanded nozzle. For a sharp edge, the expansion at 
the edge follows Prandtl-Meyer solution exactly. Away from 
the edge, the flow is governed by the three-dimensional flow 
equations. For simplicity, a sharp nozzle edge is assumed for 
all the cases investigated. 


Treatment of the External Stream 


When a jet issues into atmosphere which is still with respect 
to the rocket, it may be assumed that the pressure on the jet 
sundary is constant and equals the atmospheric pressure. 
When the rocket is in motion, however, the jet boundary 
s no longer under a constant pressure and the interaction be- 
tween the jet and the external stream must be determined in 
rder to know the conditions at the jet boundary. In general, 
he determination of the interaction is not simple. However, 

f the Mach number range of the external stream is suf- 
ciently high, the solution may be simplified by using some 
ypersonic approximation. 

In practical considerations, it turns out that when the jet 
spreading becomes significant, which is normally at high alti- 
ude, the rocket speed is usually in or near the hypersonic 
‘ange. Therefore, it appears justified to assume a hypersonic 
xternal stream in the present investigation. Under this 
issumption, the corpuscular theory of Newton, for example, 
will then be applicable and will greatly simplify the analysis 
of the interaction problem. 

According to the Newtonian approximation, the local 
boundary slope of a body in a hypersonic stream is sufficient 
to determine the pressure on the boundary at the point and its 
neighborhood. This pressure, by isentropic relationship, is 
sufficient to further determine the local boundary Mach num- 
ber of the jet flow. Knowing the Mach number, the charac- 
teristic equations will give a corresponding flow direction. 
This direction can be taken as the direction of the flow at the 
next boundary point. Thus the above procedure can be re- 


peated. 
The corpuscular theory (4) is represented by 
P. 
[6] 


where P, and P,, are the boundary pressure and the free 
stream pressure, respectively, M., the free stream Mach 
number, y,, the specific heats ratio of the external stream, 
and @ the local angle of the jet boundary with respect to the 


and y,,. The jet gas properties to be included in the data are 
the ratio of specific heats y, the boundary total pressure P,, 
and the boundary Mach number M w to which the gas expands 
around the exit corner to reach a pressure equilibrium with the 
external stream. 

In case the jet issues into still atmosphere, a constant Mach 
number M, is stored in place of M w and it applies everywhere 
on the boundary. 


Calculation of Corner Flow 


The expansion of the flow at a» around the exit corner to the 
given boundary Mach number Mw or M, can be considered 
as consisting of many small steps. During each step the 
stream line turns through a small angle and gains an increase 
in Mach number according to the Prandtl-Meyer flow. Thus 
the points @;, 82... are obtained as new data points. The 
point Bo = a is given. From each 8 a characteristic can be 
started. 


Simplifying Notations 


In order to simplify the formulas to come up in the later 
sections, some new notations are introduced here. 


In addition to these, the notations defined in Equations [8a, b, 
c, d] are also to be frequently used. 

In the following sections, difference equations are given 
which are essentially based on Equations [1, 2 and la]. 


Calculation of Interior Point from Two Noncentral Points 


The unknown point 3 can be obtained from the two known 
points 2 and 1 by the following formulas (point 1 on left- 
running characteristic, 2 on right-running characteristic) 

X 
= 
= re + A(X3 — 2 
M Ine? + 
1 


S3 = S; + 


— + AiM, + + BY Xs — Xi) + — X2) — Ci(Ss — Si) — — 82) 
= 


free stream. In the present investigation the free stream is_ 


assumed parallel to the jet axis. 


Program of Numerical Computation 


The numerical computation program employed in this in- 
vestigation is briefly described here. Only essential features 
are given, since details of any program vary, depending on the 
taste of its programmer. 


Input Data 


The nozzle exit flow is represented by several discrete data 
points @ in Fig. 1. Five pieces of information are given at 
ach data point, namely X,r, M, @and S. The a@ points may 
be on a spherical surface in the case of conical exit flow, 
may be on a plane in the case of a uniform parallel exit 
flow, and may be points after a shock. 

The data also include the external flow properties M ,, P., 
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A, + Az 


6; = 6, + A\(M; M,) B,(X3 X1) + Ci(S3 


Calculation of Interior Point From One Central and One 
Noncentral Point 


The unknown point 3 can be obtained from the two known 


WALL 


a, 


Fig. 1 Notations 
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points 1 and 2 by the following formulas, point 1 being on the 
axis 


[13] 
[14] 
S; by (10) 
M; 


~( 6) + AiMi/2 + + — Xs) — C(S; — — CAS; — S2) 


according to the foregoing procedure. Of course, the com- 
putation of corner expansion is to be omitted here. 


Results 


Cases Computed 


In order to cover a reasonably wide range of operating 


6; = 0, + 4[A,(M; — M,) + C,(S; — S,)]...... [16] 
Calculation of Boundary Point 


The unknown point 3 can be obtained from the known 
points 1 and'2 by the following formulas, point 1 being on the 


boundary 

Xo + —r — 1) 

17 
3 Saxe [17] 

1 

fs = re + — — [18] 
Ke 

P; = sin? 0, +1)........... [19] 


M; ly — P; 1 


6; = + A.(M; M2) BAX; Xe) = = . [21] 


The entropy S; is a constant along the boundary until a 
shock is encountered. The subscript ¢ refers to the stagnation 
condition. 


Calculation of Central Point 


The unknown point 2 can be calculated from the known 
point 1 as follows 


meses 
M, = M, + A E + bi( Xe s) | [24] 


Computing Procedure 


The computing procedure is to a certain extent arbitrary. 
The procedure used in this investigation is briefly described 
below. Reference is made to Fig. 1. 

1. Compute points 8; and 8... 6, at the corner by the 
Prandtl-Meyer flow equation, the number of points depend- 
ing on the mesh size desired. 

2. Calculate interior points (12, 13, . . . )of the first left 
running characteristic aby. 

3. Calculate the boundary point 6,. 

4. Repeat steps 2 and 3 for all the left running character- 
istics which originate from the data points @ until all the a are 
exhausted. 

5. Calculate the central point 30. 

6. Calculate the interior point 31. 

7. Calculate all the remaining interior points (32, 33, . . .) 
on the characteristic 30-b3. 

8. Repeat steps 5, 6 and 7 for all the downstream charac- 
teristics until a shock is encountered, or until the computation 
is sufficiently far downstream. 

9. In case the computation is stopped by a shock, locate 
the shock by tracing through the compression-wave crossings. 
Determine the flow downstream of the shock. Continue the 
calculation of flow by the method of characteristics essentially 


A,/2 + Az 


conditions of the rocket jet, a total of nine computed cases 
is presented in this paper. The results of computation are 
plotted in equal number of figures showing equiproperty lines 
and streamlines. An index to these figures is provided in 
Table 1. Every equiproperty line is marked with a Mach 
number. The three streamlines shown in most of the figures 
are the boundary streamline, the axial streamline, and a typi- 
cal streamline (dashed line) between these two. a 


Effect of External Stream an) 


The effect of external stream is shown in Fig. 11 for various 
stream Mach numbers. It is seen that the jet boundary is 
squeezed in by the external stream. However, comparison 
of Figs. 5 and 10 and other cases (not given in this paper due 
to limited space) has led to the conclusion that the flow proper- 
ties along the jet axis are practically unaffected by the ex- 
ternal stream. 


Effect of Ambient Pressure 


The ambient pressure effect is illustrated in Fig. 12. It 
will be noted that, at a given station, X/Ry = 5 in this case, 
the Mach number is a unique function of the distance R/Ro 
irrespective of the pressure ratio P./P, until the boundary is 
approached. Near the boundary, there is a sudden drop of 
the Mach number. The boundary itself, of course, flares out 
when the ambient pressure decreases. 


Effect of Nozzle Type 


Comparison of the curves in Fig. 13 quickly yields the con- 
clusion that at a given distance from the exit, within the range 
of investigation, the jet from the conical exit always has larger 
radius than the jet from the parallel exit. This difference is 
amplified by increasing the pressure ratio P./P,, and the 
Mach number M,. However, the Mach number distributions 
for the two nozzle types are not much different along the axis 
of the jet away from the exit. 
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Fig. 2 Axial section of jet 
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Fig. 13 Effect of nozzle types on boundary (y = 1.15) 


Effect of the Ratio of Specific Heats 


The ratio of specific heats y has a complex relationship with 
the nozzle area ratio e, the Mach numbers M, and M,, and 
the pressure ratios P,/P, and P,/P.. The first question in- 


volved in the study of y effect is: What is to be kept constant — 


in this study? Since the pressure ratios are known from the 
trajectory and are usually given as design conditions, it seems 
logical to investigate the y effects on the basis of fixed pres- 
sure ratios. Results of this investigation are plotted in Figs. 
l4and 15. It appears from Fig. 14 that at a given station the 
diameter of the jet is approximately proportional to the re- 
ciprocal of y. Therefore the shape of the jet boundary given 
in this paper for y = 1.15 can be used to estimate the jet 
boundary for other y values. From Fig. 15, it can be noticed 
that for constant P./P,. and P,/P, the pressure distribution. 
along the jet axis is not much affected by the value of gamma. 

Investigation has also been carried out on the basis of fixed 
values of M,and M,. The conclusion is that for a given set of 
M.and M,, the jet spread rapidly becomes less diverged when 
the value of y becomes larger. It should be pointed out, how- 
ever, that under the condition of given Mach numbers, the 
area ratio € and pressure ratios P./P, and P./P., vary rather 
rapidly with y. 


A Simplified Method for Predicting Jet Boundary 


Curves are provided for getting an approximation to the 
jet boundary for y = 1.15 and M,, = 0. Given the pressure 
ratios P./P. and P./P», the two Mach numbers M, and M, 
and the two Prandtl-Meyer angles v, and vy», can be deter- 
mined. From Figs. 16 or 18, Xpay/rz, is determined from M, 
which in turn determines Rey/r. The two quantities 
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Fig. 15 Effect of » on pressure on the jet axis; P./P,. = 52.2, 
P./P» = 62.9 


Xerwir, and Rey, are the coordinates of the boundary point 
where the jet radius is a maximum. Therefore a horizontal 
tangent to the curve is determined at this point. A tangent 
to the curve at the nozzle exit is obtained by adding the dif- 
ference of the Prandtl-Meyer angles (v, — v-,) to the wall 
angle of the nozzleiexit. A third tangent to the curve is a line 
through point (X}»/Ro, Re/2/Ro) with a slope equal to one half 
the initial angle. This point is obtaied from Figs. 17 or 18 by 
using the plot of M, vs. No2/Ro and then the plot of Re2/Ro 
vs. Xg2/Ro. The'curve is then drawn tangent to the three 
lines. 
An example follows in which M, = 2, M, = 3, y = 1.15 
and the cone angle is 15 deg. 


Ry = 7.1 from Fig. 16 
Raw = 3.5 from Fig. 16 oa 
Ro 


4 30 10 
= 
/ 4 6 8 10 16 18 
— 
+ 
26 
o7 
321 
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= — ve + 15° = 69 33 4+15 = 51° program, it is a simple matter to calculate the flow if sufficient 
Xyo/Ro = 0, Rro/Ry = 1 data are given as to the nozzle exit flow conditions, the gas 
property gamma, the pressure and the Mach number of the 
Boe = = 25.5° medium into which the jet issues. 
<allet 2 The flow graphs reported in this paper cover a fairly 
o4 Xon/Ro = 1.7; Ron/Ro = 2.3 wide range of operating conditions. A simplified method for 
predicting the jet boundary has also been derived on the basis 


Conclusions of the calculated data. If for a particular application any 

a rh 1 doubt exists as to the applicability of these results, it is recom- 

1 The spread of a supersonic jet becomes significant at mended that a special run be made on the computer for the 
high altitudes. In general, the jet flow is rather complex and particular case in question. 


far fro: 7 3 Other things being constant, an increase in divergence of 
: the nozzle exit or in the pressure ratio P./P ,, tends to increas: 
the jet divergence; an increase in the exit Mach number or in 
the Mach number of the external stream tends to decrease th: 
jet divergence; and an increase in the value of gamma tends 
to confine the jet. 


m one-dimensional. However, with the computing 


7, oe ae | 4 When the jet boundary is under finite pressure, th 
p. ane. = highest Mach number region is usually somewhere in the cor: 
6 +r Hy + Bees of the jet rather than at the boundary. Only when th 
9 pe 4 boundary pressure is zero does the maximum Mach numbe 
“, exist at the boundary. 
| 5 Ata fixed point relative to the nozzle, the gas density i 
>». -maximum when the jet boundary passes through the point. 
2) AL. i 6 In future work more emphasis should be put on gase 
sé eo with ratios of specific heats other than 1.15, and the computa 
A, Ky tion should be carried beyond the shock. 
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Technical Notes 


Vertical Re-entry into the Earth’s 
Atmosphere for Both Light and 
Heavy Bodies 


« 

7 and the subscript 1 denotes conditions at the top of an in- 
terval of constant \. The function E cannot be evaluated 

in closed form, but numerical values such as those given in 

Table 1 are easily obtained. 


R. D. LINNELL! | 
Table 1 Values of E(\p) 
Convair Division of Corp., San Diego, Np E(Xp) E(Ap) E(Xp) 
0.01 0.00995 0.2 0.1723 3. 0.4108 
Nomenclature ‘ 0.02 0 0197 0.3 0 2400 4. 0.3236 
0.08 (0.0294 0.4 0.2974 5 0.2561 
{ = cross-sectional area of body, sq ft \ | 0.04 0.0388 0.5 0.3437 6 0.2074 
‘p = drag coefficient, drag/}pv?A meine | 0.05 0.0482 0.6 0.3853 7 0. 1682 
7 = function defined by Equation [5] and Table 1 : ‘ 0.06 0.0574 0.7 0.4179 0.1428 
= acceleration of gravity, 32.174 fps? 7 | 0.07 0.0664 0.8 0.4441 9. 0.1239 
= vertical re-entry velocity, positive downward, fps ; « 0.08 0.0753 0.9 0.4675 10. 0.1139 
Vv = weight of body, lb | 0.09 0.0841 1.0 0.4841 15. 0.0694 
/ = altitude, ft 0.10 0.0927 2.0 0.4984 20. 0.0509 
3 = constant, ft?/sec? (see Equation [2]) : 
» = atmospheric mass density, slug/ft* Equation [4] is valid for any altitude interval of constant 
Subuestpte d. By applying the equation across each step of x, several 
a : i intervals of constant 8 can be used, in order to obtain a better 
B = conditions at base of an altitude interval of constant 8 7 representation of the atmospheric density. Similarly, several 
1 = conditions at top of an altitude interval of constant A of Cs. W ce A can 
The accuracy of Equation [4] is quite good, as shown in 
A SOLUTION for the velocity of a body during vertical Fig. 1 which presents the re-entry velocity for a body of 
& re-entry into the earth’s atmosphere is presented. A W/CpA = 0.1 falling from an altitude of 0.5 X 10° ft. The 
plecewise-exponential density, a piécewise-constant drag numerical solution was obtained by step-by-step integration 
coefficient and a constant acceleration of gravity are as- of Equation [1], using the Model Atmosphere (2). The 
sumed. The solution involves use of a tabulated function four-8 atmosphere (3) shown in Table 2 was used with Equa- 
of only a single variable. tion [4]. 


Allen and Eggers (1)? give an equation for the velocity of a 
body re-entering an exponential-density atmosphere vertically 


under the influence of a constant drag coefficient and constant Table 2 Approximate 6 atmosphere 


gravity acceleration. Their solution involves a series that y, ft pp, slug/ft? yp, ft B, ft?/sec? 
converges slowly for some practical values of the parameters. 0 to 0.05 X 10° 2.377 x 10-3 0 0.857 X 108 
A more tractable form is presented here for the more general 0.05 X 108 to 
case of a piecewise-exponential density and a piecewise con- 0.15 X 3.639 10-4 0.05 10° 0.695 x 108 
stant drag coefficient. 0.15 X 10° to 
The equation of motion is 0.25 10° 3.564 10-§ 0.15 0.870 x 108 
0.25 10° to 
v(dv/dy) = —g{1 — (pv?/2)/(W/CpA)}........ (1] 0.40 X 10° 7.996 X 10-8 0.25 x 10° 0.633 108 
where W/CpA is assumed to be piecewise-constant. The 
piecewise-exponential density is given by 4 
W/G,A =.1 
p = ppexp [—(y — ys)g/B]............. [2] —— NUMERICAL SOLUTION 
where the subscript B denotes the base of an altitude interval ---- FOUR-STEP @ . 
of constant 8. Ve AND W 
Using the density as the independent variable, Equation i | \ 
[1] becomes : 
where \ = 6/(W/CpA) is constant. The solution of Equa- 677— ly \ 
tion [3] is 
| 
> 4 
v? = — 2BE(Ap,)} exp [—A(o — + / \ 
29 — y) exp [—Ap] + [4] 
where WV \ 
‘ \ 
p 
E(Ap) = [(e+* — 1)/s] ds......... [5] 
2 3 4 5 
Received Dec. 13, 1957. 
Staff Scientist, Scientific Research Laboratory. Mem. ARS. 
? Numbers in parentheses indicate References at end of paper. Fig. 1 Comparison of velocity results 
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_ Some general features of the variation of the re-entry ve- 
locity can be obtained by consideration of Equation [4]. If 
the body starts from rest at a very great altitude (Ap. < 1), 
then during the first part of re-entry, where \p < 1, Equation 
[4] reduces to the expression for the vacuum free-fall velocity 


6] 
where the subscript 0 denotes the conditions at great altitude. 
If, however, the velocity at a very great altitude is large, then 
[4] reduces to 


v2? & 


These approximations are valid for \p less than about unity, 
depending on the magnitudes of the release altitude and 
velocity. The lowest altitude for which they are valid de- 
pends on the magnitude of X = B/(W/CpA); that is, on 
whether the body is light (small W/A) or heavy (large W/A). 
For high speed vertical re-entry of heavy bodies, it is known 
(1) that Equation [7] is valid almost to the earth’s surface. 

At sufficiently low altitudes (Ap > 10), the exponential 
factors are very small and Equation [4] reduces to the ex- 
pression for the equilibrium free-fall velocity 


9,2 = [8] 


For intermediate altitudes the complete Equation [4] gives 
the transition to Equation [8]. 

If the drag coefficient changes suddenly at low altitudes, 
then the new equilibrium velocity will be attained very rap- 
idly, within a few hundred feet. This is due to the exponent 
A(p — pr) = Apil(e/p1) — 1] becoming large for only a small 
increase of p/p; because the factor dp, is large. 
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_ Air Density Determination by 
Observation of a Satellite 


«ROBERT E. ROBERSON! 

Autonetics Division of North American Aviation, Inc., 
Downey, Calif. 


Previous results on the relationship between orbital _ 


eccentricity and radius for a satellite affected by air drag 
are applied to the problem of inferring air density as a 
function of altitude. 


N APPLICATION envisioned for even the most rudi- 
mentary satellites is the determination of atmospheric 
density from observations of the descent. I am not aware of 
a published description of the analytical method by which the 
air density is to be inferred, but the fact that interest usually 
is expressed primarily in the decay of altitude with time or 
revolution number implies that one of these relations is in- 
tended to be used. 
Such an approach is possible, but it does involve some 
difficulties. There are two major ones. First, there is the 
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problem discussed by Petersen (1)? of the value of the drag 
coefficient to be used for flight under satellite conditions. 
Although its approximate value can be inferred rationally, 
the coefficient inevitably enters any relationship between 
altitude and time for revolution number and must be ac- 
cepted as an additional unknown to be determined along with 
the density function. Second, there is the observation ot 
DeNike (2), based on computer solutions of the equations 
of motion over an oblate earth in the presence of air drag 
“We have an extremely complex motion . . . we have to be 
careful in using short portions of trajectory to extract such 
data as density.” This implies that an extremely precise 
determination of the density-altitude function is likely to be 
a rather sophisticated operation involving a very complex 
and time-consuming program of numerical computation. 
An alternative approach is outlined in this note. It pro- 
vides a density function without reference to the drag co- 
efficient by a very simple numerical process. These ad- 
vantages are not without cost, of course, and the proposed 
method suffers in that it may provide a cruder approxima- 
tion to the density function than can be obtained by more 
elaborate data reduction. Also, there is an implicit limi- 
tation to the case where the drag coefficient is essentially 
constant over the range of altitudes considered. However, 
tentative tests of the degree of approximation which may be 
expected indicate that the density can be estimated quite 
well this way, so that the results should be intrinsically ade- 
quate for most engineering applications. Furthermore, they 
should be useful as a first approximation in an iterative pro- 
cedure for subsequent improvements in the density function. 
The method is based on an approximate result obtained 
(3) by applying the method of variation of parameters to the 
same equations of motion as those given by Petersen. This 
is a relation between the orbital eccentricity € and the radial 
distance R of the satellite from the center of the earth, in- 
volving the air density p(R) but not the drag coefficient. 
(Strictly, R is a “characteristic” radius related to the perigee 
radius Rp by R = Rp(1 + e€).) For small eccentricities (of 
the order of « = 0.001 or less) and a = rics al atmospheric 


shell it 
_ eR 


p(R:) 


where 7 designates initial values. Given observations of 
e(R), one finds p(R)/p(R;:) directly. The constant p(R;) 
can be evaluated when drag reduces R to the point where 
p(R) can be identified with density accessible to direct meas- 
urement in the lower regions (say of the order of 100 km) of 
the atmosphere. 

If the eccentricity is initially much larger than e; = 0.001, a 
more accurate result 


p(R) _ 6k 
can be used. Here 


_2 (x) 

Io (x) 
(the functions J; and Jy) being Bessel functions of the first 
kind for imaginary argument) and 


In the right-hand side of Equation [2], one uses the first 
approximation for « and ¢ inferred from Equation [1], and by 
iteration could obtain.a more precise p(R). However, pre- 
liminary tests ofthe results against ‘“‘exact’’ solutions of the 
complete equations of motion by IBM-701 computer indicate 
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that the e(R) results for e; = 0.01 are, to within this order of 
approximation, in error by no more than about 1 per cent, 
corresponding to a density error of about 2 per cent. The 
error is not likely to become larger, except perhaps when R 
begins to change very rapidly as the satellite begins the pre- 
cipitous portion of its descent. It should be noted that 
even a considerably larger density error might be acceptable, 
inasmuch as it is now uncertain as to order of magnitude. 

It may be necessary to include the effect of the asphericity 
of the atmospheric shell arising from the oblateness of the 
carth. The analog of Equation [2] cannot be obtained unless 
the orbital eccentricity is much greater or much less than a 

ctor € related to the oblateness of the earth and defined 
= k(a?/b? — 1) | 
€ 4R sin‘ y cos 28p...... 
vhere a and 6 are respectively equatorial and polar radii of 
he earth, y is the orbital inclination (0 < y < x), and Bpis 
he angle from the ascending node to the perigee (measured 
positively in the direction of satellite motion). If xe >> €, the 
result is the same as Equation [2]. However, if xe << € 


f R dinp 


In Equation [6], the right-hand side is evaluated using the 
approximate density from Equation [1], as described follow- 
ing Equation [4]. In intermediate regions of €, one has re- 
course to an approximate first-order differential equation for 
e(R), but its solution has not been explored with the present 
application in mind. 


[p(R) eR 
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Lifetimes of Artificial Satellites’ 


ROBERT R. NEWTON? 


Applied Physics Laboratory, The Johns Hopkins University, 
Silver Spring, Md. 


Introduction 


N A recent paper (1), Henry calculated the lifetimes of 

artificial satellites. If the orbit is not originally circular, 

he showed that the life of the satellite has two separate stages: 

In the first, the orbit changes from elliptical to circular; in the 

second, the satellite spirals in from a circular orbit. Unless 

the eccentricity is originally quite small, most of the lifetime 
is spent in the first stage. 

Henry based the theories of these stages upon two particu- 
lar assumptions. In the decay from a circular orbit, he as- 
sumed that the orbit parameters remain those appropriate 
to a circular orbit, with the radius changing slowly. In the 
decay from an elliptic orbit, he assumed that the altitude at 
perigee remains constant. While the errors in these assump- 
tions are small, the lifetime is determined by only small changes 
in the orbit parameters. It is, therefore, necessary to ex- 
amine these assumptions carefully to see if they are adequate 


Received Jan. 6, 1958. 

1 This work was supported by the Department of the Navy, 
Bureau of Ordnance, under Contract Nord 7386. 

? Physicist, Research Center Division. 

3’ Numbers in parentheses indicate References at end of paper. 


May 1958 


for calculating lifetimes. We conclude that the first assump- 
tion is quite accurate, but that the second leads to an error 
which, for Sputnik I, amounts to about 40 per cent in the 
lifetime. 

We retain the assumption that atmospheric density varies 
with altitude as 


p= K exp(=—@ — {1] 


where 7 is distance from the center of the earth, R is the 
sarth’s radius, and K and c are constants. Henry uses the 
value c = 23.3km. Data on the orbital changes of Sputnik I 
are now available (2); from the rate of change of its period, 
we find that the average value of ¢ applying up to the perigee 
altitude is about 12 km. 


Lifetime in a Circular Orbit 
The equations of motion of a particle in a gravitational 
field of force, subject to a drag force pv?A, are 
= (h?/r3) — (yp/r®) — + (h?/r2)) [2] 
h = —(pA/M)h[r? + [3] 
where h = r6._ If we neglect the vertical component of ve- 
locity 7, and let h = 1/u, these become 


Using a subscript zero to denote initial values, the condition 

that the orbit be initially circular is that yurouo? = 1. 
Unless the satellite is falling rapidly, p and r, and hence wu, 

change but slightly during one period. Hence, we have 


= + Bi)...... 


where @ is small and slowly changing. Let _ 


Assuming that y is small and @ is constant, Equation [4] be- 
comes 

= —[(3 — — [8] 

From Equation [8], we see that y executes forced oscilla- 

tions; using the initial conditions, we see that the period of 


these oscillations is the same as the period of revolution. The 
vertical component of velocity, or 7, is found to be oe 


y = —[2Bro/(3 — — cos wt) 


We are interested only in the “center line”’ of the oscillations 
in *; that is, we ignore cos wt. If it is legitimate to assume 
that r and u are always connected by the relation appropriate 
to a circular orbit, we can apply this analysis at any time, and 
drop the subscript zero, obtaining 


Since 8 = w/u, this equation is found by substitution to have 
the solution 
.. [10] 


This is the relation satisfied by r and u on a circular orbit; 
therefore Henry’s assumption is justified. Since y = 0 at the 
beginning and end of each period, the satellite starts and 
ends each revolution tangent to a circle, with the speed and 
radius connected by Equation [10], which is Kepler’s law for 
the period ina circular orbit. Each revolution, the orbit 
“steps” downward; Henry’s results show that the steps in- 
crease exponentially, with most of the loss in altitude occur- 
ring during the last revolution. 


Lifetime in an Elliptic Orbit 


We begin by ‘calculating the changes in the total energy E 
and in the angular momentum Mh during one period. We 
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shall pilattaiies it use pv?A for the drag, and shall assume 
that the speed v is adequately given by r6 or h/r. For the 
change in energy, we have, obviously 


AE = ¢ = —Ah? [11] 


Neglecting r in Equation [3], we have - 
=(pA/M)(ht/r) = 

whence 


_ From AE and Af, we wish to calculate the changes in p, 
the height (from the center of the earth) at perigee, and 
q, the height at apogee. It is not accurate enough to obtain 
approximate values of AE and As from the integrals in 
Equations [11, 12], because Ap is a small difference in two 
large quantities. Instead, we must manipulate the integrals 
directly. 

E and h are connected (3) with the semi-major axis a and 
the eccentricity « by 


E = —(yuM/2a) = [yua(1 — [13] 
From these, we find by routine manipulation 
Aa = —20°(Ah*/yuM)  (p/r)a8 rm. 
Ae = —(Ah*/yuMe — €)(a/r) — (r/a)}d0 J 
p and q are related to a and by 
p =a(l q=a(l+e)........ [15] 
whence, after further manipulation, we find 
’ ) 
—(Ah*/yuM) pire) — a*(1 — €)?}d6 | 
.. [16] 


(1 + €)? — 
It is now permissible to develop approximations for Ap and 
Aq. 

Unless the eccentricity is quite small (of the order of 10~), 
most of the drag occurs near perigee. Choosing @ = 0 at 
perigee, most of the values of the integrals come from small 8. 
Keeping only the first terms in 6 which appear, but not as- 
suming that ¢€ is negligible compared with unity, Ap and Aq 
reduce to 


Ap = —(Ah?/ypM)a[(1 — €)/(1 + ©))]K exp [—(p — R)/c] X 


exp [—peO?/2c(1 + ]02d0 


Ah?/yuM)a(1 — €)—'K exp [— 


(p — R)/c] X 


exp [—peO2/2c(1 + ©) ]d0. 


These must ‘ve solved p and q. 
While we cannot solve the equations accurately unless p is 

constant, we can readily estimate the total change in p which 

occurs during the lifetime in an elliptic orbit. Eliminating n 


dq/dp = (q? — p*)/pe............... [18 | 


This can be solved exactly in terms of Bessel functions for pure 


- imaginary argument; however, if the original eccentricity €9 is 
not too large, say a few per cent, there is a simpler procedure. 
Factoring g? — p? as (¢ + p)(q — p), we can replace q + p by 

_ its average value, which is very nearly (2 + €)p. Equation 
[18] then integrates to 


= [2€opo — c(2 + X 
exp [(2 + €)(p — po)/c] + p + c(2 + &)=.. 


where jo is the initial value of p, and we have approximated 
the initial value of q by po(1 + 2¢€). The lifetime is nearly 
over when g = po; setting g = po, Equation [19] becomes, 
with x = (2 + €)(po — p)/c 

[2€0(po/c)(2 + &) — lle —xr+1=0....... [20] 


Using €) = 0.048, po = 6600 km, c = 12 km, which are 
approximate values for Sputnik I, this gives x = 3.69, and 
Po — p = 21.6 km. Over this change in altitude, the air 
density changes by a factor of 6, so that it is not legitimate to 
assume that the perigee altitude is constant. 

We can now estimate the lifetime, using Henry’s equation 


which follows his numbered Equation [16]: es 
2a) + a)“'da = 


Va(p 
exp [-(p — dn. . [21] 


The notation is the same as ours, except that a is not the semi- 
major axis; instead, a = (q — p)/2. The only place in this 
equation where changes in p are important is in the exponen- 
tial; replacing p + 2a by its approximate average pp + do = 
po(1 + €) and p + aby po(1 + €/2) 


dn = —[VW/aM/2po%(1 + €)(1 + X 
exp [(p — R)/c]da 


[19} 


From Equation [19], we get, neglecting c compared with 4a, 
and €) compared’ with 2 


= V a/€po exp (po/c) 
Substituting this into dn and integrating from a = dy = €opo 
to a = 0, we get as the lifetime in an elliptic orbit 
n = + + €0/2)V cpol 

exp [(po — R)/c]. . [22] 
This is less than Henry’s result by a factor (4/3) (1 + 6€)- 
(1 + &/2). With €) = 0.048, this factor is 1.43. Some of this 


factor arises from differences in estimating terms of order €, 
but most of it arises from the change in perigee. 


exp (p/c) 


- Applications to Sputnik I 


Analysis of Sputnik I is not the main purpose of this paper. 


grals are 


[2(1 + ©)c/pe]*”? 


respectively. 
Finally, we eliminate a, € and / in favor of pand q. Letting 
n be the number of revolutions, following Henry, Ap = dp/dn, 


and [2(1 + €)c/pe]'/? 


However, we have made some slight use of deductions from 
observations of its orbit, which should be explained. The 
most easily measured change in the orbit is the change in the 
period P. If we calculate the derivation dP/dn, which is 
directly related to dE /dn, by the methods of the preceding 
section, we get 


Aq = dq/dn, and we obtain as differential equations satisfied 
by p and q dP/dn = 
—3PV — + X 
leq/p(q — p)]'* exp — R)/c] 17) 
— 
dq/dn MAK M)a(q + p)lea/p(q — X 3 The relative error produced by neglecting here is only €/8; 


exp [—(p — R)/c] ) 


considerable algebra is needed to > show this. 
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Using the values reported in (2), we can infer the factor 
(AK/M) exp [—(p R)/c] unambiguously; this factor is 
24 xX 10-4em-. Using M = 84 kg and a diameter of 53 
em, we get K exp[—(p — R)/c], the density at the original 
perigee altitude, to be 9.0 X 10~'* gm/cm’, about 7 X 107" 
times the standard sea level value. From this, we infer that c 
is about 12 km. 

The estimate of density is uncertain by the error in the ob- 
served value of dP/dn, but this error should have a negligible 
efiect upon c. 

Having gone this far, it is necessary to go on and estimate 
the lifetime. From Equation [22], we find that the lifetime 
in an elliptic orbit is 1430 turns. From Henry’s results, we 
find the lifetime after the orbit becomes circular to be 9 turns. 
This gives a total lifetime of 1439 turns, or 96 days. 
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Vertical Ballistic Trajectories Over an 
Oblate Earth 


ROBERT E. ROBERSON! 


Autonetics Division of North American Aviation, Inc., 
Downey, Calif. 


Previous results on impact points of ballistic rockets on 
an oblate earth are modified to include the case of vertical 


trajectories. 


N A previous paper (1)? I developed expressions for the 

correction to ballistic missile impact points which result 
from the oblateness of the earth. One case was explicitly 
omitted, namely the case of vertical launch, and it was ex- 
pected that this case would have to be developed separately 
ab initio. However, further consideration of the results has 
permitted this exceptional case to be included in the previous 
theory as described below. 

Expressions for lateral deviation a,, range deviation .Agd™ 
and time deviation ur are given by Equations [35, 36, 
42], respectively, of that paper, in the general case. For short- 

ange trajectories, the results specialize to Equations [45, 46, 
47] which are in the form of the first terms of a power series 
development in the range angle Ag. (All notation is that of 
(1).) 

For a near-vertical launch case it is evident that the short 
range form suffices. As 6 — 7/2, the range approaches zero. 
All of the higher order terms in the power series in Ad can be 
neglected, for examination of the general equations shows that 
none of these can be multiplied by a coefficient that becomes 
infinite as 0 — 2/2. However, the coefficient of the series 
does have the latter property, and one is left with a dominant 
indeterminate term (tan 6) A@ in each of the short range equa- 
tions. 

To evaluate the indeterminate term, we introduce h, the 
maximum altitude of the rocket above the reference sphere 
(of radius R,) as it follows its elliptical reference path. It is 
easy to show that tan @ tan (Ag/4) = h/(R, + h), so that 


li Ah 
(tan | =; ——..... {1] 
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Using this in Equations [45-47] of (1) 


2 

R, 2 Rit+h 
a? 4h : 

a’? th 

= 1 — 3 

MT RwiRit+h ( 


= 
+ 


equation is in view of the fact that the launching 


[4] 

The last « 

velocity and maximum altitude are related by 
vy? = 2Kh ‘Rr (Rr + h) 


Next, turn to Equations [26, 32] of (1) which give latitude 
and longitude corrections to the impact point in terms of the 


variables of Equations [2-4]. Introduce AN = A," and 
AE = L‘ cos \1 as northerly and easterly corrections. 
They are given by 
12ua? h 
AN = — — [| ——— sin 2A,...............-. 5 
R, in 2A, [5] 
12ua? 
= (1 — 3sin?A,) cosA,.... [6] 


R, Wor 


where w, = WVK/R,3 is orbital angular velocity at the 
reference sphere. 

For an illustrative calculation suppose R; = a, whence 
12ua?/R, = 41.8 km and Q/a,; = 0.0588. The behavior of 
AN and AE with nominal maximum height of trajectory and 
with launching latitude is shown in Fig. 1 
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Evaluation of Flame Stability at High - 
Reynolds Numbers' 


‘a 
= diameter of injector port 


g = velocity gradient evaluated at injector wall; superscript 
t = turbulent flow 

Re = Reynolds number = (a,-d)/v 

ty = average linear gas velocity 

e = thickness of laminar sublayer; subscripts: 1 = using 


bas. ll Equation [3] for \ and one-seventh power velocity 
7 , distribution law; 2 = using Equation [3] for \ and 


one-tenth power velocity distribution law 
= friction factor 
v = kinematic viscosity of gas mixture 


EWIS and von Elbe (1) developed a theory by which the 
condition for flashback of laminar, premixed flames can 

be correlated to the velocity gradient of gas flow evaluated 
at the wall of a cylindrical burner tube. Edse (2, 3) extended 
this treatment to turbulent flames. Evaluation of the flame 
stability is straightforward for laminar flames; for turbulent 
flow, however, a solution is not readily obtainable because of 
the lack of an explicit analytical expression for the friction 
factor. The equation for the velocity gradient evaluated at 
the burner wall is obtained (4) by equating the pressure force 
that maintains the gas in a steady-state motion and the fric- 
tional force which opposes the gas flow. A general equation 
for the velocity gradient under steady-state flow conditions is: 


8d 


This equation is applicable to either laminar or turbulent 
flow fields if appropriate values of the friction factor \ are 
employed. When the Reynolds number is less than 100,000, 
the Blasius expression for the friction factor fits the experi- 
mental data quite well for turbulent flow. Then the velocity 
gradient expression becomes 


= 0.03955(Re)*/* [2] 


It has been found experimentally (2, 5) that the velocity 
gradients at flashback are dependent upon the burner diameter 
when the burners are not cooled.. However, Edse (2) showed 
that if the burner tip is cooled adequately, the diameter of the 
tube has practically no effect on the flame stability for the 
range of burner diameters employed. 

When the Reynolds number exceeds 100,000, the error in 
calculating friction factors from the Blasius expression be- 
comes quite large. When the constants of the empirical 
formula of Prandtl and Lees (4, 6) are modified as follows 


the calculated values of the friction factor fit the experimental 
data for smooth tubes to within +3 per cent for the range 3 X 
10° < Re < 10’. The deviation of values of the various expres- 
sions from Moody’s experimental data (7) is given in Fig. 1. 
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Fig. 1 Per cent error of various friction-factor equations as a 
function of Reynolds number for turbulent flow 


4 Using Equation [8] to express A, the velocity gradient equi- 
tion becomes 


= 8.38 X 10-*[Re + [4] 


Calculations of the depth of penetration of quenching for 
burners utilizing premixed gases under atmospheric and hig! 
pressure conditions (2, 3, 5) show that this distance is alway- 
less than the thickness of the laminar sublayer. From th: 
theory of turbulent flow it can be seen that the distribution o! 
linear gas velocity across a cylindrical tube varies according 
to the one-seventh power of the distance and that a lamina: 
sublayer exists in a region very close to the wall itself. The sub- 
layer thickness, obtained by considering the balance of pres- 
sure and frictional forces, the pressure drop along the tube, 
and Newton’s law of viscosity, is 


16.25d 


Values for the fricton factor can be calculated from the 
Blasius expression or from Equation [3] for very high Reyn- 
olds numbers. In the latter case the resulting expression is 
5586d 
(Re)®-78[195 + (Re)®-41] 


At a Reynolds number of 10°, the calculated sublayer thick- 
ness is 10 per cent larger when the Blasius expression for ) is 
used instead of Equation [3]; at Re = 10’, the difference in- 
creases to 63 per cent. 

According to Prandtl (4) the one-seventh power law of 
velocity distribution changes gradually to a one-tenth power 
law variation at a Reynolds number of 4 X 10°. Under these 
conditions the laminar sublayer thickness is 


2505d 
(Re)®-72[152 + 


€2 


Laminar sublayer thicknesses as calculated by Equations 
[6, 7] show that at a Reynolds number 10’ the thickness cal- 
culated by the one-tenth power law distribution is only 16 
per cent larger. Thus it appears reasonable to conclude that 
the thickness of the sublayer is relatively insensitive to the 
law governing the velocity distribution for the range of ex- 
ponents considered. 
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Effect of High Altitude Conditions ee 


Atomization Phenomena 


Cc. C. MIESSE! 


{per demands of modern defense tactics require that air- 
craft and/or rockets maintain combustion at extreme al- _ 
titudes. Since the combustion of liquid propellants is pre-_ 
dicated upon adequate atomization, evaporation and penetra- 
tion, the effect of high altitude conditions on these contribu- — 
tory phenomena is of paramount importance. 
The behavior of sprays under high altitude conditions was — 
investigated experimentally by Garner and Henny (1)? who 
used an air blast atomizer to determine that the SMD in- 
creased slightly as the pressure decreased. The data which 
they obtained with five different hydrocarbons are presented 
in Fig. 1, which reveals that the sprays of toluene and benzene 
deteriorate markedly as pressure is reduced, with lesser effects 
for the other liquids. Replotting of these data on logarithmic 
coordinates in Fig. 2 suggests that the pressure dependence — 
can be represented by 7 


Armour Research Foundation, Chicago, II. 


where the appropriate evaluation of n is noted on each curve. 

This effect of lower ambient pressures, which is accompanied 
by lower densities, was predicted theoretically by Weber (2), 
whose analysis of stream disintegration was used by Miesse 
(3) to correlate his experimental data obtained under sea- 
level conditions. In the same report, the finer atomization 
effected by higher pressures was demonstrated by comparing 
Miesse’s data with those of Kuehn (4), Lee (5) and Woltjen 
(6). The data of these investigators are presented in Fig. 3, 
which gives further indication of the decrease of drop size with 
ambient pressure. More closely akin to the air blast atomiza- 
tion is the phenomenon of secondary atomization, a simple 
criterion for which was given by Littaye (7) 


where p’ is the ambient density, U is the droplet velocity, D is 
diameter and ¢ is the surface tension of the liquid. Equation 
[2] thus predicts a value of unity for exponent n of Equa- 


tion [1], under the assumption that the density of the gas | 


Varies as its pressure. 
Received Jan. 20, 1958. 
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Fig. 2 Logarithmic variation of Sauter mean diameter with 
ambient pressure: n = negative slope 


The effect of ambient pressure on the vaporization rate for 
a single droplet is given by Langmuir (8) wy 


dm —2rDY come 


dt P 
which indicates an inverse variation with ambient pressure, 
other factors being the same. Since the SMD of a spray is de- 
fined as the single diameter which represents the vaporization 
characteristics of the distributed spray, the vaporization rate 
of the collective spray of N droplets can be expressed by as- 
suming a monodisperse spray 


dM dm 


where N can be defined in terms of the total volume V of 
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Fig. 3 Variation of relative droplet diameter with Reynolds 
number for various ambient pressures (Miesse) 


liquid in the spray 


The over-all effect of pressure on the initial vaporization of 
a spray, which provides the gaseous fuel for ignition, can be 
determined by simultaneous consideration of Equations [1, 3, 


. 5]. For kerosene which approxims ates current jet fuels, 


= 0.2, and 
6V 


dM 
dt P 


Hence it is apparent that a decrease in ambient pressure by a 
factor of ten will result in an increase in initial vaporization 
of a kerosene-type fuel by a factor of four. 

The lifetime of an evaporating droplet can be determined 
from consideration of Equation [8] which can be written in 
the form 


where X varies inversely as P, as indicated by Equation [8]. 
From Equation [6], the droplet lifetime 7 can be expressed 
as 


Simultaneous consideration of Equations [1, 3, 8] then re- 
veals that the lifetime of a kerosene droplet produced by a 
given injection system varies as P®-*, 

The penetration distance of a spray can be determined from 
an equation developed by Miesse (9) for the distance traveled 
by an evaporating droplet with an injection velocity Uo 


Since Equation [3] and the previous assumptions regarding 
the proportionality of p’ and P indicate that the product p’A 
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is independent of pressure, Equation [9] predicts that the 
penetration distance L varies as the droplet lifetime. Hence, 
from Equations [1, 3, 8] 


[10] 


again indicating that the penetration distance is decreased 
under high altitude conditions for n < 0.5. 

Actually, the evaporation rate of a decelerating droplet is 
dependent upon both the Schmidt and Reynolds numbers, as 
described by Fréssling’s equation (10) 


d(D?) 


—ro(1 + 0.3 W/Se VR) 


= —do(1 + 0.3 Wu/pD VUp'D/p) 


Since the diffusion coefficient D varies inversely as the pres- 
sure, it is apparent that the Schmidt number Sc is unaffecte:! 
by pressure variations. The Reynolds number however varie. 
as P!™ by virtue of Equation [1] and the assumed pressur:: 
dependence of the density p’. Simultaneous consideration 0: 
Equation [11] and Stokes’ drag equation by Miesse (11) lec 
to the following equations (with suitable approximation) for 
the droplet lifetime and penetration distance, respectively 


and 
L - 1 a @ 
2+ 4b 
4bA (1 — 
(1 + b)(3 + 5b) 
where 
= 0.3 WSe [14] 


and 6 < 1 is independent of pressure. Because of the small 
relative magnitude of 6, it is apparent from inspection of 
Equations [1, 12, 13] that both 7 and L vary essentially as 
P12; hence they can be expected to decrease at high alti- 
tudes, " provided n <0.5, as noted previously regarding Equa- 
tions [8, 9]. Since the effect of the Reynolds number is de- 
termined primarily by the (1 + 4A) term in the denominator 
of Equations [12, 13], it is apparent that a decrease in pressure 
leads to a decrease in the value of A, and hence to a slight 
increase in both droplet lifetime and penetration distance. 
This increasing effect, however, is of secondary importance, as 
shown by the following equations, derived from Equations 


[1, 12] 
or af KPi-™ 
| 1+ (1 + | 
K 


+ 


[15] 


Inspection of Equation [15] reveals immediately that the 
droplet lifetime, and hence the penetration distance, de- 
creases as the pressure decreases, since the second term in 
the brackets is of second order because of the smallness of both 
b and P at high altitudes. 

As a result of the foregoing analysis, it is apparent that the 
effect of high altitude, which is evidenced by lower ambient 
pressures, is threefold: (a) The initial Sauter mean diameter 
(SMD) of the spray formed by primary disintegration or 
secondary atomization is slightly larger than the correspond- 
ing SMD at sea level pressures, (b) the vaporization rate is 
greater and (c) the penetration distance is decreased as a con- 
sequence of - shorter droplet | lifetime. 


? 
| a 
i 
By 
5 | 
: 
: 
» 
* 
- 
4 


at the 
Hence, 


reased 


plet is 
eTS, 


pres- 
Tecte:| 
varies 
essur 
‘ion ¢ 


[13 


. [14] 


small 
on of 
lly as 
_alti- 
is de- 
nator 
ssure 
slight 
ance. 
pe, AS 


tions 


; the 
de- 
n in 


both 


, the 
ient 
eter 
1 or 
te is 
con- 


ION 


1 Garner, F. H., and Henny, V. E., “Behavior of Sprays 
under High Altitude Conditions,” Fuel, vol. 32, 1953, p. 151. 

2 Weber, C., “Zum Zerfall eines Fliissigkeitsstrahles,”’ 
ZAMM, vol. II, 1931, p. 136. 

3  Miesse, C. C., “Correlation of Experimental Data on the 
Disintegration of Liquid Jets,” Industrial and Engineering Chem- 
istry, vol. 47, 1955, p. 1690. 

4 Kuehn, R., “Uber die Zerstiubung fliissiger Brennstoffe,” 
Motorwagen, vols. 27 and 28, 1924 and 1925. Also NACA TM 
329, 330, 331, 1925. 

5 Lee, D. W., “The Effect of Nozzle Design and Operating 
Conditions on the Atomization and Distribution of Fuel Sprays,” 
NACA TR 245, 1932. ., 

6 Wo6ltjen, A., ‘Uber die Feinheit der Brennstoff-zerstiu- 
bung in Olmaschinen,” Technische Hochschule, Darmstadt, 
1925. 

7 Littaye, G., “Sur une Theorie de la Pulverisation des Jets 
Liquides,’’ Comptes Rendus (Paris), vol. 217, 1943, p. 99. 

8 Langmuir, I., “The Evaporation of Small Spheres,’ 
Physical Review, II, vol. 12, 1918, p. 368. 

9 Miesse, C. C., ‘Ballistics of an Evaporating Droplet,” Jer 
PROPULSION, vol. 24, 1954, p. 237. 

10 Fréssling, N., “Uber die Verdungstung |fallender Trop- 
fen,” Gerlands Beitrage zur Geophysik, vol. 52, 1938, p. 170. 

11 Miesse, C. C., “The Effect of a Variable Evaporation Rate 
on the Ballistics of Droplets,” Journal of the Franklin Institute, 
vol. 264, 1957, p. 391. 


Celestial Iconospherics, the Ultimate 
Astronomy’ 


WILLIAM J. BERGER? 


Celestial Iconospherics is defined as the study and utility 
of images formed by the refractional, diffractional and 
gravitational bending of electromagnetic radiation around 
large celestial bodies. The thesis is that, according to 
known laws of physical optics, a large spherical celestial 
body must function as a telescopic positive lens forming 
real images of almost everything in the universe; the focal 
surface of such a lens is a sphere, called an iconosphere, 
which is concentric with the celestial body. A spherical 
celestial body has iconospheres of many different sizes 
depending upon the wave length and the mode of bending 
of the electromagnetic waves. 

Physical investigation of the iconospheres requires con- 
trollable space-traveling vehicles carrying radio, image- 
transformer and television equipment. Because of the 
high optical resolving power available on the iconospheres, 
the method may suffice to detect the most remote of 
celestial objects and to detect other important objects of 
small angular diameter, such as the planets of other stars. 
Techniques of celestial iconospherics may make radio 
communication possible over interstellar distances. 


Iconospheres 

Y WELL-KNOWN laws of physical optics, there are 
three natural modes of wave-bending whereby a spherical 
body acts as a positive lens on incident electromagnetic radia- 
tion which are refraction, diffraction and gravitation. Since 
a celestial spherical body presents the same aspect to all direc- 
tions of space, the result of its positive lens action on incident 
radiation is that the body is concentrically surrounded by a 
series of spherical optically real focal surfaces, herein called 
iconospheres. The radius of any particular iconosphere de- 


pends upon the mode of bending of the radiation and upon the © 


ratio of body-diameter to radiation wave length. 


Received Dec. 23, 1957. 
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A celestial body with an atmosphere generally has icono- 
spheres of all three principal types, which, in the increasing 
order of magnitude of their radii, are refractional, diffractional 
and gravitational. A body without an atmosphere has no 
refractional iconospheres. 

A spherical celestial body is in fact a lens which looks simul- 
taneously in all directions, and its iconospheres contain 
images of almost everything in the universe. Since the re- 
solving power of a lens, as measured by the smallest angular 
diameter detectable, depends upon the ratio of radiation wave 
length to lens diameter, it follows that the optical resolving 
power on the iconospheres is extremely high, at least for the 
refractional type, and hence that the iconospheres contain in- 
formation of the greatest importance in astronomy. 


Compound Telescopes 


The fact that most moons are inside the iconospheres of 
their respective planets causes the existence of other kinds of 
optical foci. The moon may be entirely inside a conical 
bundle of rays which have been deflected by the planet and 
which have a focus on an iconosphere of the planet. Some of 
these rays, in passing the moon, are further deflected by dif- 
fraction and gravitation and thus arrive at other foci, so that 
the planet-moon system acts as a compound lens. Such a 
compound iconospherical system, however, because of its 
narrow view of the universe, is generally less useful than the 
ordinary complete iconospheres. 

Physical Investigation 

Since the minimum radius of an iconosphere is thousands of 
times greater than that of the celestial body which generates 
it, and since the iconospherical images can be oberved only by 
a detector placed on or outside the iconosphere, physical in- 
vestigation of the iconospheres evidently requires use of a 
controllable space-traveling vehicle of some sort. The vehicle 
may be manned or not. It might be a robot ship cruising for a 
long time in a gravitational orbit just outside an iconosphere, 
recording its observations by image-transformer and _tele- 
vision methods and occasionally correcting its orbit by small 
bursts of power until eventually covering all areas of interest. 
The first attempt at physical investigation may well be much 
simpler: A ship sent out from the earth to make a loop 
around some particular image (for instance, that of Mars) on 
a refractional iconosphere of the earth and thence somehow to 
return its observational data to earth. 

Image detection on visual iconospheres, i.e., those caused by 
light in the visual spectrum (including also ultraviolet and in- 
frared), is complicated by the fact that the detector may often 
be illuminated not only by the desired image but also either 
by direct sunlight or by sunlight reflected from the celestial 
body generating the iconosphere. Such difficulty can be prac- 
tically eliminated by methods of shielding, filtering, and image- 
transformation or television. Since the ratio of visual icono- 
spherical brightness to albedo of the generating body is a 
maximum for moons in eclipse by their respective planets, the 
visual iconospheres of the Galilean moons of Jupiter are very 
useful. Gravitational iconospheres have radii so great in 
comparison to diameter of their generating body that useful 
visual images probably occur on gravitational iconospheres 
generated even by the sun itself since at such distances the 
sun appears only as a moderately bright star.* 

Because of fewer shielding difficulties, image detection on 

iconospheres in the radio spectrum is usually simpler than in 
the case of the visual spectrum. 
Utility 
Celestial iconospherics has important uses in astronomy of 
the solar system, in cosmology, and in the possibility of in- 
terstellar radio communication. 

3 Zwicky’s 20-year-old theory of gravitational lenses is most re- 
cently described in: Zwicky, F., “Morphological Astronomy,” 
Springer-Verlag, 1957, pp. 215-218. 
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Before visiting other planets, we could learn much about 
them by studying their images on iconospheres of the earth. 
To illustrate in the simplest way the order of the magnitudes 
involved, let us consider the visual iconosphere of earth for 
the red light commonly seen when the sun is on the horizon. 
Since this light is refracted by nearly half a degree of arc in 
passing through our atmosphere, its focus is an iconosphere, 
having a radius of some 458,000 miles. On this iconosphere, 
the image-diameters of the sun, Jupiter, Mars and the nearer 
stars are about 4000 miles, 89 miles, 33 miles and 500 ft, re- 
spectively. With such resolving power, the planetary systems 
of other stars are definitely observable, and it may even be 
possible to determine whether or not the number of galaxies is 
finite. 

Iconospheres in the radio spectrum are useful not only in 
“seeing” across great distances through the intergalactic dust 
and gas but also in communicating over interstellar distances. 
Certain kinds of interstellar radio signals may be too weak to 
be detected on earth itself but they have strong foci on the 
radio iconospheres of the earth; thus by detecting alien radio 
signals of characteristic information-carrying type, we may 


_ prove that intelligent beings exist in other solar systems. 


By placing a radio transmitter on a suitable radio icono- 
sphere of the earth and sending a signal toward the earth, we 
cause the signal to pass as a plane wave beyond the earth. 
Thus, by the search-light principle, we may expect the signal 
to travel with virtually undiminished intensity over a vastly 
greater distance than it could achieve in its normal form of a 
diverging wave. Moreover, the cross-sectional area of the 
search-light beam produced by the earth-lens is obviously 
much greater than any which can be realized by man-made 
antennas alone, and hence by this method we may achieve 
communication over interstellar distances. Radio icono- 


_ spheres of the sun are even more effective for this purpose. 


By the same principle, experiments in interstellar radar are 
Of course, a peculiarity of two-way interstellar ex- 
periments is the long time interval (perhaps nine years) be- 
tween transmission and reception of signals. 

Conversely, since all of these things are possible to us, they 
are possible also to intelligent beings in other planetary sys- 
tems, and so we may ev ventually communicate with such 
beings. Celestial iconosphe rics is the master key to explora- 


tion of the universe. Since it uses the most powerful of all 
_ lenses, namely the large celestial bodies themse Ives, it may in- 


deed be called the “ultimate astronomy.” alee: 


This brief essay leaves many questions not explicitly an- 
swered, but the purpose has been to introduce and define 
celestial iconospherics, to forecast the proper conduct of the 
exploration of the universe, and to open the way to further 
discussion. 


Computation of Rocket Step Weights 
to Minimize Initial Gross Weight 


RICHARD P. TEN DYKE! 


The Ramo-Wooldridge Corp., Los Angeles, anes 


The technique of dynamic programming is applied to the 
problem of finding optimum step weights of multistage 
rockets. If the design velocity requirement of the rocket 
system is determined, the velocity may be allocated to the 
various stages in such a manner that the initial gross 
weight is minimized. This approach permits the study of 
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rockets with three and more stages, where the specific 
impulses for the stages are unequal, and where the jettison 
weight relations are nonanalytic. It is also possible to de. 
termine the relative flatness of optima with respect to any 
particular staging ratio with slight additional calculation. 
The approach outlined in this paper supplements but does 
not replace the use of calculus in solving the same problem 
where certain simplifying assumptions may be necessary. 


Definitions? 


Step Rocket: A rocket system employing the technique of jettison- 
ing structure weight during flight to increase performance. Also 
multistage rocket. 


Stage: The period of time starting just after the jettison of 
one structural package and lasting until just after the jettison of 
the next, or the time from take-off to just after the first jettison 
(first stage) or, if no weight is jettisoned at last-stage burnout, 
from just after the last jettison to burnout. 


Rocket weight: The initial gross weight at the start of a stage. 


Step Weight: The sum of the propellant ejected during a stage 
and the weight jettisoned at the end of the stage. Also, the change 
of weight from the start of one stage to the start of the next.? 


All weight which remains after the jettison at the 


Stage Payload: 
Rocket weight of next stage. 


end of the stage. 


Rocket Payload: Payload of last stage. 


Nomenclature 

f,{(V] = minimum weight of rocket k achieving velocity V 

g = gravitational constant 32.2 fps? 

I; = specific impulse of stage 7 (thrust + flow rate of fuel) 

N; = ratio of initial thrust to weight of rocket 7 

n = number of steps (the rocket payload may be con- 
sidered the (n + 1)-step) 

ri = burnout mass ratio W;/(Wi+1 + oiw;) 

V; = velocity ideally added by rocket 7; rocket 7 design 
velocity 

Voo = actual burnout velocity of rocket system 

V; = design velocity of rocket system; also Vi 

vi = velocity ideally added during stage 7 

6V = velocity lost by rocket to effects of atmosphere and 
gravitation 

W; = initial gross weight of rocket 7 

Wi = rocket system payload weight 

Ww; = initial gross weight of step z 

Wri = that portion of stage 7 jettison weight which is con- 
stant (independent of step size or thrust) 

aw; = portion of stage 7 jettison weight dependent on step 
weight 

B:N;W; = portion of stage 7 jettison weight dependent on thrust 

Wi = total jettison weight of stage 7 

Subscripts 


i = denotes ith step, stage or rocket — 
k = denotes kth rocket 


Introduction 


ANY authors (1-4)* have considered the problem of de- 
termining the optimum stage sizes as a straightforward 
minimization using calculus, and their results are very useful. 
To make hand computation possible certain simplifying as- 
sumptions are necessary, such as: (a) Two stages or less, (b) 
same specific impulse in all stages, or (c) structure propor- 
tional to size of step. For many problems the assumptions are 
quite reasonable, and the results are found to compare favor- 
ably with actual practice; however, problems arise where the 
simplifying assumptions must be verified or the problem 
solved by use of an alternate approach. 


Such an alternate approach is outlined here. The use of 


? Fig. 1 shows graphically the definitions of rocket and step. 
5 Numbers in parentheses indicate References at end of paper. 
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dynamic programming (5), a technique found useful in alloca- 
tion problems, permits the operator to use structure factors 
which are nonanalytic with discontinuities. It also helps the 
operator to determine flatness of optima, which is valuable if 
the rocket must be designed away from the “ideal optimum” 
for any reason. It is particularly helpful in the study of 
rockets with three or more stages. 

Hand computation with this approach is lengthy, but suit- 
able for rockets of three or fewer stages. The technique is 
ready-made for machine computation if a large number of 
eases are to be investigated. 


Suitable Optimization Parameters 


The well-known rocket equation 


yields a relationship, independent of time, between the ve- 
locity ideally added to a rocket during a stage and the stage 
burnout mass ratio. This independence does not hold for 
actual surface-to-surface rockets because two time-dependent 
efiects, atmosphere and gravitation, erode some of this 
velocity. The effects of atmosphere are: (a) Drag which in- 
creases with increasing velocity (and therefore acceleration) 
at low altitudes and (b) reduction of engine thrust due to at- 
mospherie pressure. The effect of gravitation is to penalize 
the use of low acceleration in nonhorizontal trajectories. 

However, studies of rockets flying realistic trajectories have 
yielded gratifying results. If we consider a given rocket with 
engines of constant flow rate, the velocity lost to gravitation 
and atmosphere, 6V, is principally dependent on the ratio of 
initial thrust to initial weight for each stage, N;, and is inde- 
pendent of range above 500 nautical miles. Fig. 2 presents 
this result. Thus, the simplified rocket equation, [1], may 
be used for design purposes if it is corrected by the additive 
constant 6V. The sum of the stage design velocities is 
termed the rocket design velocity V:. 


Ve = = Zig inri = Von + 8V........... [2] 


Of course range, not burnout velocity, is usually the de- 
sired result. The same studies show that for the maximum 
range trajectory, the burnout velocity essentially defines the 
range irrespective of other parameters. Although the actual 
burnout velocity required for a given range will depend on the 
burnout conditions, namely altitude, angle and down-range 
distance, deviations from nominal values will be comparatively 
small. Since burnout conditions are also greatly dependent 
upon the thrust-time program, corrections can be considered 
a part of the quantity 6V. A curve of burnout velocity as a 
function of maximum range appears in Fig. 3. 

It appears from the foregoing that the treatment of this 
optimization as an allocation (of velocity) problem is a prac- 
tical approach where Equation [2] is the constraint and V, 
is determined by the range requirement. 


The General Optimization 


The general optimzation assumes that the designer has the 
following ‘‘inputs” to the problem: (a) The specific impulse 
of fuel for each stage and (b) the jettison weight for each stage 
as a function of step weight or step propellant weight, of stage 
thrust and/or of velocity added or burnout mass ratio of each 
stage. For (a) the specific impulse is assumed constant for the 
duration of a stage but may vary in value from stage to stage. 
For (b) the function may be in analytic or tabular form except 
that if it is of certain forms computation is simplified. 

‘quation [1] yields a relationship between the design 
velocity of a stage and the stage burnout mass ratio. The in- 
verse relation also holds. 
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Fig. 3 Burnout velocity vs. maximum range 


From the definitions 


Solving for the weight of propellant 


= (" ‘) (Wit1 + [5] 


t 


From item (b) of the assumptions we say the jettison weight 
is a function A of the items listed 


Combining [5] and [6] 


1 
(: ) (Wiss + wi) + Ni, ve]. [7] 


For the nth stage of an n-stage rocket, the value W,,+1 is 
simply the payload for the last stage, Wz, and is known. 
Solving for w; we find it dependent upon three quantities. 


= 
10,000 
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per These equations yield a dynamic programming problem in two 


i variables, v and N. This can be reduced to a problem in one 

-_-variable if N is held fixed. For considerations of performance 

only, the ratio of initial thrust to initial weight of the first 

stage should fall between one and four. For succeeding stages 
values above one are preferred, but lower values are ac- 

ceptable for nearly horizontal thrust vectors. Conditions 

- other than optimality may dictate a choice of N, but no at- 

tempt will be made here to consider these constraints. 
Assume values for N; are selected, then 


Wy = 


Now define 


= minimum ferilV — ve), ve] + lV — .. [10] 
O<u<V \ f 
where — is the minimum initial gross weight of the 


+1 rocket which achieves a velocity of V — and w[fey1 
—- [V — x], x] is the weight of step k which achieves the veloc- 
ity The sum (V — + is the velocity V, and is se- 
: lected to minimize the sum of the weights. This cannot be 
accomplished unless the weight of the payload, rocket k + 1, 
is known as function of the velocity added by the payload. 
This is known for the nth stage where the payload is known 
and the velocity added by the payload is zero. 


11] 


When f,,[V] is found and plotted against V, then f,—1 can 
be found using [10] and so on to f,. Equations [10, 11] are 
_ the typical recurrence relationships found in allocation prob- 

zs ye of one variable and are solved quickly by digital 


i Sample Case 


If n is small, say three or less, and the function w[Wi+1, v;] 
_ is simple to compute, a hand solution is practical. Consider 
the case where the stage jettison weight is 


+ + Wry = h[wi, N,, Wi+i]. ° [12] 
From [7] 


= 


where a;, 8;, and w,; are constants. 


1 
) (Witi + wi) + aw + + wi) + wri 


Substituting [3] 

Tg 


If N; are selected, the weight W; is dependent only on the 
_ velocity gained, v;, and the recurrence relations are 


+ (1 — afer — vx 


= minimum [14] 
+ (1 — On)W 


Flatness of Optima 


For two-stage rockets, a single plot of f, vs. ; is sufficient 
to solve for a minimum and to get a picture of its flatness. For 
three or more stages additional plots will be required as func- 
tions of particular staging ratios. Fortunately, the informa- 
tion is available from the calculations. 
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Surface Combustion in Dissociated Air 


SINCLAIRE M. SCALA! 


General Electric Co., Philadelphia, Pa. 


Nomenclature 
A = atom 
By, = X,/m 
= u/u 
Map = combustion rate 
My = molecular weight of the nth species 
= Xn My, Mean molecular weight 
M, = first molecular species 
M, = second molecular species 
M, = 
Le = Lewis number 
Yr = Prandtl number 
= temperature 
U, v = tangential and normal components of velocity 
x, ¥, ro = body oriented coordinate system 
Xn = mole fraction of nth species 


= stagnation point velocity gradient 


= T/T, dimensionless temperature 
= viscosity 

= density 
= similarity variable 


Subscripts » 


= outer edge of boundary layer ; = 
= atomic species 

= first molecular species 
second molecular species 
= nth species 

= stagnation point 

wall 

n = denotes differentiation with respect to 7 


HE need for structural materials which are suitable for 
use as aerodynamic members at high temperature initiates 


_ the requirement for analytical techniques for studying the in- 


= 


Ordnance Systems Department. 


teraction between different classes of materials and chemically 
reacting flow fields. In order to explore the features involved 
in the theoretical determination of the combustion rate of a 
material immersed in a stream of dissociated air, two dif- 
ferent models are considered for the combustion of graphite. 

Some typical boundary layer profiles, near a surface at 
which combustion occurs, are shown in Fig. 1. It is noted 
that the physical extent of the boundary layer is not known 
a priori. Rather, in the absence of external forces, the various 
boundary layer profiles are governed by the presence of trans- 
port processes whose influence appears in the form of dimen- 
sionless groups such as the Reynolds, Prandtl and Lewis num- 
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bers. It is also noted that in the general treatment of com- 
pressible boundary layers with surface reactions, the equa- 
tions of change can be uncoupled and solved individually only 
for special circumstances, such as isothermal flows over flat 
plates. 

In addition, unlike surface catalysis (1, 2),? in which the 
macroscopic gas velocity vanishes at the surface, during sur- 
face combustion the normal component of velocity is not zero 
at the surface but is related to the rate at which the surface is 
consumed. 

Reactions at surfaces generally occur in five steps, the slow- 
est of which will determine the over-all rate of the process. 
Ii is usually considered that these successive steps include the 
transport of reactants to the surface, adsorption of the gases, 
combustion, desorption of the products or reactants, and 
transport of the liberated products back to the main stream 
by convection and diffusion. 

If the surface reactions are extremely rapid, the rate-con- 
trolling step will be the diffusion process. That is, the com- 
bustion rate will depend on the quantity of reactant that can 
diffuse to the surface, in the presence of a counter-diffusion of 
combustion products. 

For simplicity, it will be assumed here that gas phase reac- 
tions can be neglected; i.e., the boundary layer is frozen. 
Then, upon utilizing a similarity transformation, the stagna- 
tion point boundary layer equations for a three-component 
compressible chemically reacting gas, in which homogeneous 
reactions are frozen, become (3) 


in + MijX jn) + JPrBiy 
jy + MuiXiq)|y + fPrBjy = 0.....(2] 


which are the equations for the conservation of the ith and jth 
diffusing species 


1] Me 


which is the momentum equation, and 


fPr6, [Le 5X in + MijM ;X = (@... [4] 


which is the energy equation. In deriving Equations [1-4] 
which constitute a coupled ninth order system of nonlinear 
differential equations, it was assumed that the Prandtl and 
Lewis numbers were constant, and that the specific heats 
were constant and equal for all species. 

In the theoretical models considered here, it was assumed 
that by appropriate “lumping” of those species whose thermal 
and chemical behavior is sufficiently similar, dissociated air 
and the gaseous combustion products may be represented by 
a three-component system, denoted by subscripts 7, 7 and & 
(see Fig. 2). 

The two models considered were: (a) Nitrogen is an inert 
diluent (model I). (b) Oxygen and nitrogen react completely 
to form CO and CN, i.e. ; 


X iy = = 0 (model IT) 


The relevant surface boundary conditions were derived and 
solutions were obtained (3) by means of a Reeves analog com- 
puter (REAC). 

The combustion rate is given by 


du. 


and a universal solution for all combinations of free stream 
conditions and wall temperatures is shown in Fig. 3 for model 
II, a Prandtl number of 0.75 and a Lewis number of 1.2. 
Note that since Xx, is always zero (see Fig. 1), and since 
=X, = 1, the specification of Xi, determines the free stream 


n 
dissociation level. That is, when Xi, = 0, the free-stream air 


2 Numbers in parentheses indicate References at end of paper. 
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Fig. 1 Typical combustion boundary layer profiles 


Model I Model II 


Fig. 2 Models for the combustion of graphite. Model'I: no 
nitrogen reactions; model II: 


nitrogen reactions allowed _ 


.20 
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Te 
Te 
Fig. 3 Variation of combustion rate with surface temperature, 
model II 


is undissociated and when Xi, = 1.0, it is entirely dissociated. 

The solutions obtained for model I yielded values for m, 
which were of the order of 25 per cent of those obtained for 
model IT. a 
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Technical Comments 


Comment on ‘Stresses in Rocket 


Photoelastic Design Data for Stresses in Rocket Grains,”? 
we have arrived at unreasonably high calculated circum- 
ferential stress on the internal surface of a case-bonded pro- 
pellant charge. Moreover, the sign of the estimated stress 
was opposite from what one would expect. Subsequently this 
matter was called to the attention of one of the authors, M. L. 
Williams, who has explained the discrepancy as being the re- 
sult of omitting certain terms in expanding the denominator 
of Equation [8]. This omission limits the applicability of 
Equation [9] to cases in which the ratio Ec/E is of the order 
one. 

Since this ratio may be several hundred, even for low tem- 
perature conditions and high rates of strain under which the 
mechanical properties of the propellant approach those of an 
elastic material, caution should be used in employing the 
simplified Equations [8a and 9]. 


Received Feb. 7, 1958. 

1Ordahl, D. D., and Williams, M. L., ‘‘Preliminary Photoelas- 
tic Design Data for Stresses in Rocket Grains,” JET Prorui LSION, 
vol. 27, June 1957, pp. 657-662. a 


? Senior Technical Specialist. Mem. ARS. ~ iv) 


Corrections to ‘Comments on ee 
_ Powered Flight Trajectory 
of a Satellite’! 


4 
BURTON D..FRIED? 


The Ramo-Wooldridge Corp., Los Angeles, Calif. 


HE optimum steering program for the powered flight 
trajectory of a satellite launcher is given by* 


where y is the angle between the thrust vector and the hori- 
zontal at time ¢ and a and b are constants whose values depend 
on the thrust and mass programs. The conditions for the 
validity of [1] are: (a) The missile is treated as a point par- 
ticle. (b) The distance covered during powered flight is small 
compared to the earth’s radius R., so that the gravitational 
force field may be taken as uniform during powered flight. 
(c) Aerodynamic forces are neglected. (d) The quantity W 
to be optimized (i.e., made stationary) is a function only of 
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N ATTEMPTING to apply Equation [9] of “Preliminary 


‘burnout velocity and altitude, for instance, the orbit altitude 
at perigee (if W depends also on the range at burnout, tan y 
red becomes instead a quotient of two linear functions oi 
‘The of powered flight is short to® 


95min 


where g is the gravitational acceleration. 

In a recent note,! Edelbaum remarks that the analysis lead- 
ing to [1] does not include “the effect of centrifugal force. 
which is significant for satellite launching vehicles,” thereby 
inferring that “centrifugal force” effects may invalidate [1]. 
This is incorrect. So long as conditions (a) through (e) are 
satisfied, the program [1] will give a stationary solution for a 
W of the form specified in (d) independent of whether the 
velocity during powered flight becomes equal to or even 


greatly exceeds the standard satellite velocity v, = VV gR.. 

Edelbaum’s error is an understandable one, stemming from 
a common misinterpretation of “centrifugal force.”” The 
point is simply that 7f (0) ts satisfied, then we may analyze the 
powered flight trajectory in Cartesian coordinates. The total 
force acting on the missile consists only of the thrust force 
and the gravitational force; there is no “centrifugal force’’ 
(the earth’s rotation being neglected, of course), regardless of 
the magnitude of the missile’s velocity. Of course, one may 
choose to use polar coordinates (r, 6) in place of (x, y) to 
describe the motion, as Edelbaum does, and one then finds in 
the equations of motion a “centrifugal” pseudo-force, mv’r/r?. 
However, so long as (b) is satisfied the choice of coordinates 
clearly cannot affect the answer and a correct calculation in 
any coordinate system must again lead to [1]. 

To demonstrate this analytically, we use the general result*® 
that if condition (b) is absent and the equation of motion for 
the point mass missile is assumed to be 


VR./29 = 


a 


where 7 is the acceleration due to thrust, N is a unit vector 
along the thrust, and @ is the gravitational acceleration 


then the optimum direction of thrust is along the vector }(¢) 
defined by the differential equation 


[4] 


and boundary conditions which depend on the particular op- 
timization problem at hand. For instance, if we wish to make 
stationary a function W(r, v,) of the burnout position and 
velocity for given initial conditions fo, Vo and given n(t), then 
we have 


As a result of the vector methods, it is clear that when condi- 


3 Fried, B. D., “On the Powered Flight Trajectory of an Earth 
Satellite,’ Jer PRopuLsION, vol. 27, June 1957, pp. 641-643. 

4 Lawden, D. F., “Optimal Flight Trajectories,” Jet PRopu.- 
SION, vol. 27, Dec. 1957, p. 1263. 

5 Fried, B. D., “Space Technology,’’ Wiley and Sons, Chap. 2b 
(to be published). 
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For records you can read atonce... 
new Kodak Linagraph Direct Print Paper | . 


This paper is designed for use when you must know 
what’s happening while it’s happening . . . in such ap- 
plications as rocketry, control, computing, product 
design, component testing and nuclear research. It re- 
quires no processing, comes out of the recording in- 
strument ready to read. 

To get such instant monitoring of high-speed varia- 
bles you eliminate the powerful energy amplification 
factor of photographic development. Traces must be 
put down by a high-intensity u-v light source. 

Kodak Linagraph Direct Print Paper gives you a 


EASTMAN KODAK COMPANY 
Graphic Reproduction Division 
Rochester 4, N. Y. 


May 1958 


wide range of recording speeds—from 0 to 3,000 cps. 
Records have adequate permanence for extended stor- 
age. If used under sunlight or fluorescent lamps they 
should be stabilized. The surface of the paper takes 
notations in pencil or ink. Records are suitable for 
use in recording system readers. 

Several major manufacturers now offer recording 
instruments based on the characteristics of Kodak 
Linagraph Direct Print Paper. We’ll gladly send you 
their names on request. If you want more detailed in- 
formation on the paper, we'll send that too. 

ye 
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tion (b) is satisfied so that we can set Va = 0, then }, = 


and for the boundary conditions [5] 


x Vo,W + (4 
When 0W/dz, = 0, this is just equivalent to [1]. 


neglecting the distance from launch to burnout after comput- 


ing Va rather than before, gives 


[6] Va 


One may, 


as stated above, resolve [2, 4] into polar coordinates and then 


solve for A, and Ag, but in the resultant algebraic complexity Az 


it is possible to lose sight of the simple result [6]. 


Of course, when Va ¥ 0, [1, 6] are no longer valid and we i, 
must solve the coupled system of Equations [2, 4]; but this is 


simply a consequence of the failure of condition (b) and not 


of any “centrifugal force” effects. 


It is interesting to note 
that the first correction to the Va = 0 result, obtained by 


If (4; — tb) <1/ 
reduces to [6]. 


and hence, from [4 and 5] in Cartesian coordinates 


(w? = g/R.).......-. [7] 


—(0W/d7) cos w(t — th) + 


[(OW/dz2;)-sin w(t — 


—(W/dny) cos h w(t — + 
[(OW/dy)-sin h /2 w(t — 


2 w, as required by condition (e), then [8 | 


LAUNCHING 
RAMP 


MOTHER 


TRANSMITTER 


2,826,378 

Apparatus for radio control of guided mis- 

siles (2,826,378). J. N. Childs Jr., Wynd- 

moor, Pa. 

Mother equipment is a directional trans- 
mitter for providing a very narrow rapidly 
sweeping radio beam. When the time 
intervals between successive signals at the 
missile are equal, the missile is traveling 
down the center line of the swept sector 
directly toward the target. 

Control system for guided missiles (2,824,- 
710). A. C. Hale, Bloomfield Township, 
Mich., assignor to U. S. Air Force. 

Homing missile of the glider type in 
which aerodynamic control is effected by 
elevons. Homing information is provided 
in the form of right-left and up-down error 
signals and means for measuring differ- 
ences between roll and pitch angles, 


Vehicle for testing control systems at 
supersonic speeds (2,824,711). H. H. 
Porter, Washington, D. C., assignor to 
U.S. Navy. 

Rocket engine powered guided missile 
having a forward compartment for con- 
taining guidance equipment to be tested. 
Device for indicating the relative move- 
ment of the earth and a heavenly body 
R. H. Farguhar, Philadel- 
phia, P 

Time. ring within a transparent globe 


_ MISSILE 


~ 


~ 


SWEPT SECTOR Se, 


and fixed to a shaft ine, through its axis 
of rotation. A second shaft carries a head 
representing a satellite rotatable in a great 
circle path adjacent to the outer surface 
of the globe. 

Aerodynamic valves (2,825,203). J. H. 
Bertin and A. R. Mihail, Asnieres, France, 
assignors to SNECMA Co. 

Pulsejet unit with resonant firing com- 
bustion chamber. Valveless unilaterally 
conductive air inlet device offers greater 
resistance to backflow of hot gases from 
the chamber than the flow of fresh air in- 
to the chamber. 
= propulsion units (2,825,204). M. 

adosch, F. G. Paris, J. Le Foll and J. 
Bertin, Neuilly-sur-Seine, France, assign- 
ors to SNECMA Co. 

Combustion apparatus for generating a 
continuous flow of motive gas. Two di- 
vergent jet pipes have means for urging the 
gas to flow successively and intermittently 
toward one inlet and then the other. 
Combustion devices especially suitable for 
gas turbine engines and propulsion units 
for aircraft (2, 825,205). é& co Paris, 
France, assignor to SNECMA Co. 

Jet propulsion engine with a flame- 
holding device comprising an auxiliary 
compressor which discharges fluid into the 
combustion chamber at an inclination from 
the axial flow, forming transverse fluid jets 
creating a wake zone in which a flame can 
be held. 


Polycyclic aldehydes as hypergolic fuels 
(2,826,037). C. R. Scott and J. V. Smith, 
Bartlesville, Okla., assignors to Phillips 
Petroleum Co. 

Separate and simultaneous injection of 
an oxidant and fuel into contact with each 
other in the combustion chamber of a 
reaction motor. Fuel is composed of an 
aldehyde of hydrogen atom, a methyl! rad- 


George F. McLaughlin, Contributor 


ical, a hydrogen atom and an alkyl radical! 
having carbon atoms. 

Gas inlet structure for combustion cham- 
bers (2,826,039). P. F. Ashwood, Farn- 
ham, England, assignor to Power Jets (Re- 
search and Development) Ltd. 

Channel extending from the combustion 
chamber wall inwardly across the combus- 
tion reaction flow path. The open side of 
the channel faces upstream and a part of 
the gas is directed into it from the chamber. 
Vortex tube free air thermometer (2,826,- 
070). H. C. Box and L. S&S. Packer, 
Eggertsville, N. Y. 

Cylindrical inlet conduit with a portion 
terminating in a right angle defining an 
elliptical-shaped opening for exhausting 
the air. A fairing maintains a static pres- 
sure region in the air stream adjacent to 
the exhaust opening. 

High speed sampler (2,826,075). R. J. 
Vernier, Baltimore, Md., assignor to U. S. 
Army. 

Atmospheric sampler comprising a filter 
disk locked in a hollow holder ring to which 
a suction tube is attached. Suction ap- 
plied to the ring causes ejection of the disk. 
Torpedo suspension band (2,826,119). 
A. H. oa, Attleboro, Mass., as- 
signor to U. 8. Navy. 

Apparatus for suspending a missile from 

an aircraft having a bombrack and mov- 
able supporting hook. 
Mechanism for ejecting rockets from an 
aircraft (2,826,120). J. M. Land and R. P. 
Graham, Stratford, Conn., assignors to 
Chance Vought Aircraft, Inc. 

Pair of drums mounted within the body 
of an airplane. Plungers eject rockets 
from opposite sides of the aircraft through 
openings in the drums. 
Signal device (2,826,144). 
Tokyo, Japan. 

Tubular waterproof container having a 
storage chamber with inflammable pow- 
der. Release of a spring-tensioned trigger 
causes a rod composed of ignitable material 
to rub against an abrasive surface on the 
container wall. 

Alloy steel for use at high temperatures 
(2,826,380). W. R. Kegerise, — 
Pa., assignor to Carpenter Steel C 

Alloy capable of being hardened ad show 
a martensitic structure, good corrosion 
resistance, ductility and life at tempera- 
tures up to 1200 F, high impact strength, 
and machineable at room temperature. 


Yukichi Inata, 


Eprtors Note: Patents listed above were selected from the Official Gazette of the U.S. Patent Office. 


Printed copies of patents may 


be obtained from the mnenasd al Patents, sieaaseee 25, D. C., at a cost of 25 cents each; design patents, 10 cents. 
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OF A SERIES 


New Olds-developed machine 
makes wheel balancing three 
times more accurate! 


Out-of-balance wheels and tires are 
not only a source of annoyance and 
tire wear, but also in extreme cases, 
a detriment to safety by causing ex- 
cessive shimmy at higher speeds. 


To virtually eliminate this problem, 
Oldsmobile engineers, in conjunction 
with the General Motors Research 
Section, have developed a machine 
that automatically balances every 
wheel and tire with a degree of pre- 


cision not previously possible on a 
production basis. With this equip- 
ment. balancing is now accurate to 
2 inch-ounces, or approximately three 
times more precise than before. 


The heart of such accuracy is an auto- 
matic electronic’ computing device. 
After the tire and wheel are located 
on a delicate sensing table, supported 
on an air bearing, four differential 
transformers signal the out-of-bal- 
ance to an electronic computor. This 
computor then resolves the vector 
forces and a signal of the proper mag- 
nitude and direction is transmitted 
to the stamping head which automati- 


TIPPING THE BALANCE IN YOUR FAVOR 


cally revolves to the correct location 
on the wheel. The stamping head then 
prints the correct weight, accurate to 
.25 ounce. The entire assembly is 
then moved to a station where the 
weights are attached. 


It has often been said that “Olds 
really knows how to put a car to- 
gether.” This reputation grew from 
a sincere concern for just such little- 
noticed details. A warm welcome 
awaits you at your Olds dealer’s. He in- 
vites you to try a 58 Olds on the road. 


OLDSMOBILE 
GENERAL 


DIVISION 
MOTORS CORP. 


.-.Famous for Quality Manufacturing 
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Pioneer in Progressive Engineerin 


we have 


We have developed new tech- 
niques, new methods, new 
processes that effect production 
economy so necessary to a success- 
ful missile program. 


Here at Newbrook you will find 
men with experience gained 
from doing ...a modern plant 
with up-to-date equipment... 


pecializing in 
@ Motor Cases 
Solid and Liquid 


propellants = precision inspection to meet your 
© Jato most exacting quality control 
requirements. 
@ Nozzles 


And most important, Newbrook 
specialization results in strict 
reliability ! Let us help you with 
your Missile Hardware problems. 


@ Plenum Chambers 


@ Blast Tubes 

@ Fuel 


Phone or } 


NEWBROOK MACHINE CORPORATION 
20 Mechanic Street Phone 45 
SILVER CREEK, NEW YORK 


‘vite 


. 


the “Know-How!” 


Contains carbon, silicon, manganese 
chromium, molybdenum, copper, nitrogen 
and iron. 

Bearing cooling system (2,826,895). W. | 
English, Northport, N. Y., assignor io 
Fairchild Engine and Airplane Corp. 

Gas turbine engine housing with a dis- 
charge opening and an exhaust nozzle 
forming an air chamber in which the e 
haust gases create an aspirating effect for 
cooling, the tail cone interior. 


Rocket dispensers (2,826,961). A 
White, Washington, D. C., assignor to 
ACF Industries, Ine. 

For installation on high speed aircraft, 
an outer shell with a conical nose spaced 
from a tubular body. A dispensing tube 
in the shell has a conical sheath with its 
rear end attached to a bulkhead. The 
sheath expands radially to permit passage 
of a rocket between the nose and body. 
Container for high pressure fluids (2,827,- 
195). ‘I’. F. Kearns, Arlington, Va 

Spherical liner over which a wire mesh 
is wound in layers. Tension in the wires is 
increased in each layer, moving outward] 
from the liner. 

Gas turbine aircraft powerplant (2,827,- 
759). 1B. W. Bruckmann, Dayton, Ohio, 
assignor to the USAF 

Engine with contraflow air-fuel combus- 

tion system. An air-fuel escape nozzle is 
directed forwardly in the combustion 
chamber to deliver the mixture directly 
against the oncoming compressed air. A 
gas exhaust jet leads from the outlet side 
of the turbine to free air. 
Combined anti-icing and generator cool- 
ing arrangement for a gas turbine engine 
(2,827,760). F.C. 1. Marchant and J. E. 
Bell, Mangotsfield, England, assignor to 
Bristol Aero ingines, Ltd. 

Air outlets from hollow spider arms dis- 
charge compressed air to a double-walled 
entry nose bullet housing the electrical 
generator. 

Constant speed reaction motor (2,827,- 
762). M.C. Towns Jr., Morristown, N. J., 
assignor to Reaction Motors, Inc. 

Means for admitting a propellant to the 
reaction chamber in regulated amounts to 
automatically maintain a given speed, the 
supply being unaffected by centrifugal 
forces. 

Ion tracer airspeed indicator (2,827,786). 
B. Boyd, R. G. Dorsch and G. H. Brodie, 
Bedford, Ohio, assignors to U. S. Navy. 

Apparatus for indicating the speed of an 
airstream flowing along a surface by means 
of an ion bundle created in the airstream 
by corona discharge. 

Gyroscope (2,827,788). A. M. Campbell, 
Weston, Mass., assignor to Minneapolis- 
Honeywell Regulator Co 

Internally threaded hollow stator 
mounted within a hollow rotor. Static 
balance is achieved by an_ externally 
threaded mass in the stator moved along 
the rotor axis upon rotation of the mass. 


Control arrangements for craft operable 
in space (2,827,789). K. A. Oplinger and 
I. M. Holliday, Venice, Calif., assignors to 
Westinghouse Electric Corp. 

Gyroscope with an inertia member 
mounted for rotation about a pivot axis. 
The pointed end of an electromagnet core 
remains within the boundary of an open- 
ing in an armature over the range of rela- 
tive movement of the frame members. 
Planetarium projector (2,827,827). A. N. 
Spitz, Lansdowne, Pa. 

Spherical mirror centered in a hemi- 
spherical perforated globe. High intensity 
light source in the zenith position in the 
globe directs rays to the mirror, passing 
through perforations | and reflected to the 
globe, 
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AiResearch units power the contro/s of 
America’s most potent defense weapons 


(2,827,- 
1, Ohio, 


ombus- 
ozzle is 
bustion 
directly 
air. <A 
let side 


cool- 
engine 
d J. E. 
mor to 
Key defense and population centers are now 


being ringed with batteries of Army Nike 
Hercules missiles to deter or destroy aggressors. 
Supplying power for flight controls is the 
AiResearch auxiliary power unit pictured above, 
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ectrical 


2,827,- 
»N.J., now in production. 
ee As a member of the Army-industry team pro- 
ints to ducing the Nike Hercules (Army Ordnance, 
ed, the Western Electric-Bell Telephone Laboratories 
rifugal and Douglas Aircraft). AiResearch was chosen 
7 786 to design, develop and manufacture this vital 
Benalla. accessory power source for the missile because 
avy. of nearly two decades of experience in light- 
1 of an weight turbomachinery. 
eee This experience includes applications utiliz- 
a ing solid propellants, liquid mono-propellants, 
npbell bi-propellants, atomic power, cryogenic gases as 
apolis- well as gasoline and air. AiResearch’s ability for 
high capacity production as well as in research 
stator and development. made it the logical choice. 
rahe Garrett’s AiResearch divisions have also 
along designed systems and components for 18 other 
nass, missiles and rockets in the U.S. defense arsenal. 
erable We invite your inquiries. 
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Jet Rocket 
Engines 


Propellant Vaporization as a Criterion 
for Rocket-Engine Design; Calculations 
Using Various Log-Probability Distribu- 
tions of Heptane Drops, by Richard J. 
Priem, NACA TN 4098, Oct. 1957, 29 pp. 

Design and Evaluation of a Turbojet 
Exhaust Simulator, Utilizing a Solid- 
Propellant Rocket Motor, for Use in Free- 


_ Flight Aerodynamic Research Modeis, 


Aerodynamics of Jet 


by Carlos A. deMoraes, William K. 
Hagginbothom Jr. and Ralph A. Falanga, 
NACA RM L54115, Dec. 1954, 25 pp. 
(Declassified from Confidential by au- 
thority of NACA Res. Abstracts 120, p. 11, 
10/3/57. 

Investigation of Blade Failures in a J34, 
Eleven-Stage Axial-Flow Compressor, by 
Howard F. Calvert, Andre J. Meyer Jr. 
and C. Robert Morse, NACA RM ¥52116, 
May 1953, 24 pp. (Declassified from 
Confidential by authority of NACA Res. 
Abstracts 120, p. 8, 10/3/57.) 

Jet Engine Thrust Reversers, by Rob- 
ert C. Kohl and Joseph 8. Algranti, SAE 
J., vol. 65, no. 11, Oct. 1957, pp. 38-40. 

Seporinent Investigation of a Rocket 
Source for a High Temperature Wind 
Tunnel, by F. K. Hill and H. A. Wallskog, 


Johns Hopkins Univ., Applied Phys. 


 Lab., CM-907, July 1957, 46 pp. 


Dynamic Response Characteristics of 
Gas Turbines, by S. L. Soo and William 
W. S. Charters, ASME Prep. 57-F-22, 
Sept. 1957, 17 pp. 

General Design Considerations for 
Smaller Gas Turbines, by W. T. von der 
Nuell, ASME Paper 57-F-13, Sept. 1957, 
17 pp. 

British Aero Engines, Aeronautics, vol. 
37, Sept. 1957, pp. 102-127. 

Designing the World’s Most Powerful 
Turbojet, by Randolph Hawthorne, Avia- 
tion Age, vol. 28, Oct. 1957, pp. 26-31. 

Scale Effect in a Gas Turbine Combus- 
tion Chamber Fitted with a Swirl Atom- 
izer, Part 1, by M. V. Herbert and J. A. 
Bamford, Combustion and Flame, vol. 1 
no. 3, Sept. 1957, pp. 360-378. 

Mach 4 Turbojets, by J. 8. Butz Jr., 
Aviation Week, vol. 67, no. 17, Oct. 28, 


1957, pp. 86-89. 


Propelled Vehicles 


Theoretical Aerodynamic Force and 


Moments for a Slender Body Moving 


Through the Flow Field of a Slender Body, 
by Henry W. Woolard, Cornell Aeron. 
Lab. Inc., Rep. QC-910-C-15, Oct. 1956, 
pp 

Laminar Boundary Layer on a Spinning 
Circular Cone in Supersonic Flow at a 
Small Angle of Attack, Part II, by Martin 
Fiebig, Cornell University, Grad. School of 
Aeron. Engng. (AFOSR-TN-57-68; AS- 
bia AD — Feb. 1957, 14 pp. 


Eprtor’s Note: Contributions from E. R. 


oe Eckert, J. P. Hartnett, T. F. Irvine Jr. 


and P. J. Schneider of the Heat Transfer 
Laboratory, University of Minnesota, are 
gratefully acknowledged. 


Inviscid Hypersonic Flow over Blunt 
Nosed Slender Bodies, by Toshi Kubota, 
Stanford Univ., Stanford, Calif., Preprints 
of papers, Heat Transfer and Fluid Me- 
chanics Inst., 1957, 439 pp., pp. 193-210. 

Measurement of the Moment of In- 
ertia of Missile-Type Bodies, by A. B. 
Schwartz, S. Malick and J. R. Friesen, 
Aircraft Engng., vol. 29, Sept. 1957, pp. 
271-274. 

Formulas for the Velocity Potential of 
Linear Subsonic and Supersonic Flow for 
Not-so-slender Bodies at Zero Lift, by 
Friedrich Keune, Z. Flugwiss., vol. 5, Aug. 
1957, pp. 243-247 (in German). 


Heat 
Flow 


Condensation Interfaces in Two Phase 
Flows, by P. Chiarulli and R. F. Dressler, 
Nat. Bur. Standards Rep. 5341, April 1957, 
14 pp. 

On Menard Inserts in Supersonic Noz- 
zles, by J. L. Amick, 8S. A. Harrington and 
H. P. Leipman, Univ. Michigan, Engng. 
Res. Inst., Tech. Rep. CF 2656, May 1957, 
8 pp. 

Mass Transfer Cooling of a Laminar 
Boundary Layer by Injection of a Light 
Weight Foreign Gas, by E. R. G. Eckert, 
P. J. Schneider, A. A. Hayday and R. M. 
Larson, Univ. Minnesota, Inst. Tech., 
Mech. Engng. Dept., TN 14 (AFOSR- 
TN-57-323; ASTIA AD 1382395), June 
1957, 45 pp. 

Growth of Drops by Condensation, by 
Ignace I. Kolodner, N. Y. Univ., Inst. 
Math. Sci., Rep. IMM-NYU 215 (2 ASTIA 
AD 121943 ), March 1955, 28 pp. 

Evaporation of Fuel Sprays, II: Ex- 
perimental Work, by G. W. Bensor and 
R. J. Brisebois, Canada, Nat. Aeron. 
Estab., Lab. Rep., LR-182 (ASTIA AD 
121802), Nov. 1956, 15 pp. 

The Thermal Properties of Carbon Di- 


Transfer and Fluid 


oxide, Nitrogen, Air, Hydrogen and 
Helium, by A. W. Nicklin, Gt. Brit., 
Atomic Energy Authority, Industrial 


Group, TN 36, Nov. 1956, 9 pp. 

The Sound Generated by Interaction of 
a Single Vortex with a Shock Wave, by 
G. I. Ram and H. 8S. Ribner, Preprints of 
papers, Heat Transfer and Fluid Me- 
chanics Inst., 1957, pp. 1-22. 

The Chemical Kinetics of Air at High 
Temperatures; a Problem in Hypersonic 
Aerodynamics, by Saul Feldman, Pre- 
prints of papers, Heat Transfer and Fluid 
Mechanics Inst., 1957, pp. 173-192. 

An Investigation of Stagnation Point 
Heat Transfer in Dissociated Air, by 
Peter H. Rose and F. R. Riddell, Pre- 
prints of papers, Heat Transfer and Fluid 
Mechanics Inst., 1957, pp. 231-248. 

An Experimental Investigation of the 
Secondary Flow Occurring in a Compres- 
sor Cascade, by W. D. Armstrong, Aeron. 
Quart., vol. 8, Aug. 1957, pp. 240-256. 

The Smoothing of Temperature Fluc- 
tuations by Radiative Transfer, by Ed- 
ward A. Spiegel, Astrophysical J., vol. 126, 
July 1957, pp. 202-207. 

Heat Transfer to Boiling Methanol, by 
A. J. Lowery Jr. and J. W. Westwater, 
Indust. Engng. Chem., vol. 49, Sept. 1957, 
Pt. 1, pp. 1445-1448. 


M. H. Smith, Associate Editor, and M. H. Fisher, Contributor 
The James Forrestal Research Center, Princeton University 


Combustion, Fuels and 
Propellants 


What’s All This Talk about Exotic Pro- 
pellants? by John F. Tormey, Missiles 
and Rockets, vol. 2, Sept. 1957, pp. 93-94 

Fluorine: Tamed for Rockets, by H. R. 
Neumark and F. L. Holloway, Missiles 
and Rockets, vol. 2, Sept. 1957, pp. 97, 
99-100. 

Liquid Propellants for Manned Rocket 
Aircraft, by William Mitchell, Missiles 
and Rockets, vol. 2, Sept. 1957, pp. 109- 
110, 112, 114. 

Dimazine Comes of Age as Rocket Fuel, 
by William G. Strunk, Missiles «nd 
Rockets, vol. 2, Sept. 1957, pp. 116-118, 
120. 

LOX—Mainstay Oxidizer for the Mis- 
sile Age, by William H. Thomas, Missiles 
and Rockets, vol. 2, Sept. 1957, pp. 131 
132. 

An Investigation of the Behavior and 
Reaction Mechanisms of Nitric Acid- 
Hydrocarbon Flames, by M. H. Boyer 
and P. E. Friebertshauser, Combustion 
and Flame, vol. 1, no. 3, Sept. 1957, pp. 
264-280. 

Extraction of Ions from a Flame, by J. 
Deckers and A. van Tiggelen, Combustion 
and Flame, vol. 1, no. 3, Sept. 1957, pp. 
281-286. 

Predicting the Laminar Flame Speed in 
Gases with Temperature-Explicit Reac- 
tion Rates, by D. B. Spalding, Combustion 
and Flame, vol. 1, no. 3, Sept. 1957, pp. 
287-295. 

One-Dimensional Laminar Flame 
Theory for Temperature-Explicit Reac- 
tion Rates, by D. B. Spalding, Combustion 
and Flame, vol. 1, no. 3, Sept. 1957, pp. 
296-307. 

Flame Decomposition of the Propyl 
Nitrates, by J. Powling and W. A. W. 
Smith, Combustion and Flame, vol. 1, no. 3, 
Sept. 1957, pp. 308-320. 


Notes on Problems in Combustion In- 
stability of Rocket Motors, by T. Paul 
Torda, Aeron. Engng. Rev., vol. 16, no. 11, 
Nov. 1957, pp. 34-37, 47. 

Shift to Solid Propellants, Aviation 
Week, vol. 67, no. 17, Oct. 28, 1957, p. 34. 

The Premixed Ozone-Cyanogen Flame, 
by A. G. Streng and A. V. Grosse, J. 
Amer. Chem. Soc., vol. 79, no. 20, Oct. 20, 
1957, p. 5583. 

Studies in Bomb Calorimetry, VII: 
Principles for the Nagy of Calorimeter 
Systems, by R. Mott, Fuel, vol. 36, 
Oct. 1957, pp. 147-268. 


Rocket Thrust and the Energy Power of 
Ergols, by P. Glansdorff, A. Jaumotte and 
J. Passelecg, Fusées Recherche Aeron., 
vol. 2, no. 3, July 1957, pp. 205-212 
(in French). 


Stabilization of Concentrated Solutions 
of Hydrogen Peroxide, by W. C. Schumb, 
Indust. Engng. Chem., vol. 49, no. 10, 
Oct. 1957, pp. 1759-1762. 

Borane High Energy Fuels, /nteravia, 
vol. 12, no. 10, Oct. 1957, p. 1023. 


Density of Liquid Ozone, by Robert I. 
Brabets and John M. McDonough, /. 
Chem. Phys., vol. 27, Oct. 1957, pp. 880- 
882. 
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New extreme-high-temperature lubricants for missiles 


7 and supersonic aircraft SHELL ETR GREASES 
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One of the serious lubricating problems _ tures. Co-operation with representatives _ special thickener—an organic vat dye— _ 
faced by designers of missiles and super- of bearing manufacturers and military _ which has exceptional heat stability and — 
sonic aircraft has been solved by scientists personnel resulted in a completely new _jelling efficiency. 
at Shell Research Laboratories. class of greases—SHELL ETR GREASES. If you are presently in the market for — 

The problem: to find a grease which These greases can easily withstand tem- _—an_ultra-high-temperature-range grease, — 
would permit components to operate with _ peratures up to 600°F. They give superior _ we will be glad to provide more informa- 
certainty under extreme high tempera- lubricating performance because of a tion on Shell ETR Greases. a 


SHELL OIL COMPANY 


50 West SOth Street, New York 20, N.Y. 
100 Bush Street, San Francisco 6, Calif. 


May 1958 
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Aluminum mandrels for 
forming solid fuel propellant are 
now being cast in production 
by the unusual foundry methods of 
Morris Bean & Company. While 
we assume there is no present 
need for a mandrel as large as the 
one on the left, it can be cast. 
Currently we are working 
on solid and hollow mandrels 
up to 8 feet long. Their smooth 
surfaces and accurate contours 
eliminate much difficult machining; 
cost is drastically reduced. 
In addition to large size, 
we would be happy to 
explore with you ways to 
produce intricate star-lobes. 
Telephone or write. 
Morris Bean & Company, 
Yellow Springs 6, Ohio. 


Materials of Construction 


Effect of Overheating on Creep-rupture 
Properties of S-816 Alloy at 1,500°F, 
by John P. Rowe and J. W. Freeman, 
University of Mich., NACA TN 4081, 
Dec. 1957, 75 pp., diagrs., photos., tabs. 

Effect of Overheating on Creep-rupture 
Properties of HS-31 Alloy at 1,500°F, 
by John P. Rowe and J. W. Freeman, 
University of Mich. NACA TN 4083, 
Dec. 1957, 78 pp., diagrs., photos., tabs. 

The Useful Heat Capacity of Several 
Materials for Ballistic Nosecone Construc- 
tion, by Jackson R. Stalder, NACA TV 
4141, Nov. 1957, 19 pp., diagrs. 

Preliminary Investigation of the 
Strength and Endurance of Plansic- 
impregnated Fiberglas Compressor 
Blades, by Donald F. Johnson and Andie 
J. Meyer Jr.. NACA RM E54127:, 
Jan. 1955, 21 pp., diagrs., photos. (D»- 
classified from Confidential by authority 
of NACA Res. Abs. 122, p. 9, 12/3/57.) 

Buckling at High Temperature, by N. JJ. 
Hoff, J. Roy. Aeron. Soc., vol. 61, Nov. 
1957, pp. 756-774. 

Study of Chemical Factors Affecting 
Corrosion of Carbon Steel A1S1 1020 by 
Liquid Phase Fuming Nitric Acid, by John 
B. Rittenhouse and David M. Mason, 
Calif. Inst. Techn., Jet Propulsion Lai. 
Mem. 20-147, Oct. 1957, 31 pp. 

Strength of Cylinders, by W. R. D. 
Manning, Ind. & Engng. Chem., vol. 49, 
Dec. 1957, pp. 1969-1978. 

Stress Concentration in Heavy-walled 
Cylindrical Pressure Vessels, by J. G. 
Faupel and D. B. Harris, Ind. & Engng. 
Chem., vol. 49, Dec. 1957, pp. 1979-1986. 


Instrumentation, Data 
Recording, Telemetering 


French Telemetering Equipment, by J. 
Idrac, North Atlantic Treaty Organization, 
Advisory Group for Aeron. Res. and Dev., 
Rep 33, Feb. 1956, 18 pp. (in French). 

Weapon System Requirements for Air 
Data Computers, by A. K. Hammell, 
Aeron. Engng. Rev., vol. 17, Jan. 1958, pp. 
32-35. 

Transistor Inverters Supply 400-Cycle 
Power to Aircraft and Missiles, Aviation 
Age, vol. 28, Dec. 1957, pp. 112-118. 

Investigation of the Conditions in the 
Plasma in Mercury-vapor Arc Discharges 
and Oscillating Electron Ion Sources, 
by Max Hoyaux and Paul Gans, Aleliers 
Construc. Elec., Charleroi, Rep. EOARDC 
TR 54 (AFOSR TR 56-23), (AD 88353), 
1956, 51 pp. 


Guidance Systems and 
Components 


Bomarc Guidance Spawns New Testing 
Techniques, by Henry P. Steier, American 
Aviation, vol. 21, Dec. 30, 1957, pp. 32-33. 

Theoretical Investigation of the Per- 
formance of Proportional Navigation Guid- 
ance Systems—Effect of Missile Con- 
figuration on the Speed of Response, by 
Marvin Abramovitz, NACA RM, A52- 
J22, Jan. 1953, 20 pp. (Declassified from 
Confidential by authority of NACA Res. 
Abs. 121, p. 14, 11/4/57.) 

Application of Statistical Theory to 
Beam-rider Guidance in the Presence of 
Noise, I: Wiener Filter Theory, by EI- 
wood C. Stewart, NACA RM A55E11, 
Aug. 1955, 40 pp. (Declassified from Con- 
fidential by authority of NACA Res. 
Abs. 121, p. 15, 11/4/57.) 
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EADY AND 


ELIABLE 


. . - Requires no refrigeration ... Easier to store than liquid oxygen 


Performance of Nitrogen Tetroxide as an oxidizer for many fuels is comparable 
to that of liquid oxygen. In convenience it is far superior to liquid oxygen. 

N20. can be shipped easily and stored indefinitely without refrigeration in 
ordinary carbon steel containers. It is a dense, mobile liquid that is noncorrosive 
if kept dry. The quantity of oxygen contained per unit volume of N2Qx is 1.01 
Kg/liter at 20°C. 

N20sz is available in tonnage quantities from Allied’s Hopewell, Virginia plant. 
For experimental purposes, N20 is 
available in 125 lb. cylinders and 2000 lb. 
llied containers. Write for technical data, and 

information on prices and delivery. 
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N,0, Oxidant for liquid 


rocket propellants _ 


Molecular weight 

Boiling Point 

Freezing Point 

Latent Heat of 
Vaporization 

Critical Temp. 

Critical Pressure 

Specific Heat of Liquid 


Density of Liquid 
Density of Gas 


Vapor Pressure 


Ethanolamines « Ethylene Oxides Ethylene Giycols Urea Formaldehyde « U. F. Concentrate— 
85 «Anhydrous Ammonia*e Ammonia Liquors Ammonium Sulfate Ammonium Nitrate « Sodium 
Nitrate Methanol « Nitrogen Solutions « Nitrogen Tetroxide Fertilizers & Feed Supplements 


Nitrogen Division ¢ Department NT4-12-1 + 40 Rector Street, New York 6, New York 


92.02 
21°C 
—11.3°C 


99 cal/gm @ 21°C | 
158°C 
99 atm 
0.36 cal/gm 
—10 to 20°C 
1.45 gm/ml at 20°C | 
3.3 gm/liter 
21°C, at 1 atm 
2 atm at 35°C 
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Pressure Vessel 


Technology 


and now Titanium... ~ 


ad 


Spheres of titanium alloy represent 
the optimum configuration and 
material of maximum strength/weight 
for the storage of gases and liquids 
under extremes of pressure and 
temperature. 

Rheem Aircraft, presently supplying 
pressure units to current ICBM pro- 
grams, is directing its years of experi- 
ence to the development and pro- 
duction of all types of vessels and 
propellant chambers. 


RHEEM MANUFACTURING CO. 
AIRCRAFT DIVISION 


11711 woodruff avenue, downey, california 


Engineers: Join Rheem in challenging technical work. 


Servomechanisms and 
Controls 
Servo Analysis Charts, by Erwin Bieser 


and Samuel Adler, Electronics, vol. 30, 
Dee. 1, 1957, pp. 173-176. 


Flight Vehicle Design and 
Testing 


Dynamic Response of Beams and Plates 
to Rapid Heating, by B. A. Boley and A. D. 
Barber, J. Appl. Mech., vol. 4, Sept. 1957, 
p. 413. 


Some NACA Research on Effect of 
Transient Heating on Aircraft Structures, 
by J. E. Duberg, Trans. ASME, vol. 79, 
July 1957, pp. 1014-1018. 

Achievement of High Heat Fluxes in a 
Wind Tunnel, by Antonio Ferri and Mar- 
tin H. Bloom, J. Aeron. Sci., vol. 24, 
Oct. 1957, pp. 772-773. 

Some Structural Penalties Associated 
with Thermal Flight, by J. W. Mar and 
L. A. Schmit, Trans. ASME, July 1957, 
pp. 990-1004. 


Space Flight 


Man in Space Needed to Recapture 
Lead, by Evert Clark, Aviation Week, 
vol. 67, no. 20, Nov. 18, 1957, pp. 26-27. 

What the Future Holds for the Earth 
Satellite, by Robert P. Haviland, J. 
Astronautics, vol. 4, no. 3, Autumn, 1957, 
p. 41 

Space Vehicles as Tools for Research in 
Relativity, by Fred Singer, J. Astronautics, 
vol. 4, no. 3, Autumn, 1957, p. 49. 

Long Term Human Confinement in 
Space Equivalent Vehicles, by Charles 
A. Dempsey, Franklin D. VanWart, 


| John H. Duddy and Jack K. Hockenberry, 


J. Astronautics, vol. 4, no. 3, Autumn, 
1957, p. 52. 

After the Sputniks, by Robert Saundby, 
Aeroplane, vol. 93, Nov. 22, 1957, pp. 
783-784. 


_ Astrophysics, Aerophysics 


eet, 


Atmospheric Temperature Observa- 
tions to 100,000 Feet for Several Climato- 
logical Regions of the Northern Hemi- 
sphere, by H. B. Tolefson, NACA 7N 
4169, Nov. 1957, 26 pp. 

Experimental Aspects of General Rela- 
tivity Theory, by Louis Witten, J. Astro- 
nautics, vol. 4, no. 3, Autumn, 1957, p. 
46. 

The Ratio of Helium and Hydrogen 
Abundances in Planetary Nebulae, by 
John 8. Mathis, Astrophysical J., vol. 125, 
Nov. 1957, pp. 493-502. 

High-dispersion Spectra of Mars, by 
C. C. Kiess, C Corliss, Harriet k. 
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GIANT CAPACITOR SYSTEM POWERS ADVANCED 


Lockheed experiments on electrical power 
sources for space craft are using electrical 

7 energy stored in a giant capacitor system, 

and Scientists from Lockheed 7 : one of the largest in the aircraft industry. 

This propulsion research—based on 
Missile Systems—where expanding electrohydrodynamics— is an example of the 
advanced work being conducted in Division 
missile programs insure laboratories near Stanford University in 


lectrical, chemi 
more promising careers the electrical, chemical and nuclear fields 
The giant capacitor bank is capable of 


a 5 million kilowatt jolt of electricity, 

used to fire Lockheed’s new “hotshot” wind 
tunnel to 18,000 degrees Fahrenheit 

with speeds up to 15,000 miles an hour. 


Research in such sophisticated areas by 
March = Division scientists has earned Lockheed 
leadership in missiles. Nose cone reentry 


‘ali 

a “ studies, for example, have helped speed 
delivery dates on Polaris, Lockheed’s 

E submarine fleet ballistic missile. 

ents. With expansion on our IRBM, Earth 4 

Rev., Satellite, and other advanced programs, 


_ new positions have opened for 
Mag- engineers and scientists... positions that 


a call for the utmost in creative ability. | 
These areas include: Propulsion, 
Aerodynamics-Thermodynamics, 
pp. Guidance, Electronics, Flight é 
Controls, Structures, Design. 


If you are qualified, you are invited to write 
Research and Development Staff, Palo Alto 
ilies, 20, California for further information. _ 
Inter- 
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Discussing characteristics of the capacitor 
system are research scientist Dr.Edward Fisher 
of the Propulsion Staff, right, and Dr. Terry E. 
Turner, Hot Shot Tunnel Group Leader. 
The two groups are engaged in a joint effort 
on advanced propulsion experiments. 


= 


ics 
serva- 
imato- 
Hemi- 
TN 
A TN 
EXPERIMENTS ON ELECTRICAL PROPULSION i 
e, by 
y 
1. 126 
¢ 
phys:- 
84. 
rer to 
J 
y Jr, 
r 


NOW IN 


VOLUME PRODUCTION 
...the AT C.A.E. 


(USAF TYPE MA-1A). 


PORTABLE 
STARTING UNIT 
for 
LARGE JET 
AIRCRAFT 


MODEL 141 
TURBO-COMPRESSOR 
ENGINE 


Typical of the fine results of Continental development is 
the TC-106 portable starting unit for large jet aircraft. This 
aa advanced new model, with a high performance turbine 

compressor as its heart, weighs one-third less than its 
_ predecessor, yet has 17 per cent higher output, and in 
_ addition, other important qualities: greater mobility, less 
_ noise, and a completely automatic control system. .. . It 


_is now in volume production at the Continental Aviation 
Engineering Toledo plant. 


rey 


C.A.E. gas turbine models—the J69-T-9, the J69-T-2, and the J69- 
T-19A are being built for Cessna’s T-37A twin jet trainer, Temco's 


TT-1 Navy jet trainer, the Beech jet Mentor trainer, and the Ryan 
Q-2A Fire Bee target drone. 


CONTINENTAL AVIATION & ENGINEERING CORPORATION 


12700 KERCHEVAL AVENUE, DETROIT 15, MICHIGAN 


SUBSIDIARY OF CONTINENTAL MOTORS CORPORATION 
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tion in relation to environment including 
temperature extremes, surface erosion an | 
surface penetration. 
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1956. (Monograph 2.) Discusses desig! 
orbit, propagation of electromagnetic wav» 
signals, and nuclear power source for ‘‘Astro 
—an artificial satellite time-and-radio orbit 

Navigation by Satellites. Missiles an: 
Rockets 1: 48-52, Oct. 1956. Indicates tha 
an artificial satellite time-and-radio orbi 
celestial navigation system that is simple, 
yet reliable, may be feasible. 

Ley, Willy 

Rocket Fuel Research Paves Way for Spac 
Flight. Indus. Labs. 7: 100-101, July 195¢ 
Discussion of liquid and solid rocket fuels 
advantages and disadvantages of each, and 
their applications to Project Vanguard. 


Ludwig, G. H., and Van Allen, J. A. 
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ment for the Minimal Earth Satellite. J 
Astronautics 3: 59, 60, 90, Autumn/Winte 
1956. Discussion of a proposed experiment 
in which a single Geiger-Mueller counter 
would be used to study the geographical and 
temporal variations of the primary cosmic- 
ray intensity. The counter would incor- 
porate a system for storing its output during 
the orbit of the satellite, and for a short-time 
read-out of this information. 
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What Can We Expect of the Man-Mad 
Moon? Reader’s Dig. 69: 162-166, July 
1956. Truths and untruths behind one of 
science’s most ambitious ventures. 


The Martin Co., Baltimore, Md. 

Effect of the Earth’s Oblateness and Atmos- 
phere on a Satellite Orbit, by George Fosdick 
and Marvin Hewitt. 41 pp., Aug. 1956. 
(Rpt. ER-8262.) Presents results of inves- 
tigations to determine the effects of the 
earth’s oblateness and atmosphere on a satel- 
lite orbit. 

Mazur, D. G. 

Telemetering and Propagation Problems of 
Placing the Earth Satellite in Its Orbit. Inst. 
Radio Engrs. Proc. 44: 752-754, June 1956. 
Reviews the problems and general planning in 
connection with the earth satellite launching 
vehicles. 


Mengel, J. T. 

Tracking the Earth Satellite and Data 
Transmission by Radio. Inst. Radio Engrs. 
Proc. 44: 755-760, June 1956. Description 
of ‘‘Minitrack’’ which is to be used to track 
and send data from Vanguard. 

Merrill, H. J. 

Satellite Tracking by Electronic Optical 
Instrumentation. In Van Allen, J. A., ed. 
Scientific Uses of Earth Satellites, pp. 29-38, 
Ann Arbor, University of Michigan Press, 
1956. Essential differences between satellite 
tracking and astronomical observations or 
missile tracking warrant improvements in 


| optical instrumentation. These are associ- 


ated with superior time constants, back- 
ground elimination, utilization of the opti- 
mum spectral region, and the superior 
characteristics of the electronics in the servo 
tracking loop. 


Missiles and Rockets 
International Geophysical Year Earth Satel- 
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Dr. H. H. Kurzweg, pioneer German rocket scientist who is now Associate Director for Aeroballistics at NOL, 


Bia Na interprets nose-cone data obtained in supersonic wind tunnel test. 


FROM UNDERSEAS TO OUTER SPACE—WEAPONS OF THE 


U. S. Navy Photo 


NAVAL ORDNANCE LABORATORY GUARD OUR FREEDOM 


Basic and applied research—and the development of 
advanced weapons for underwater, surface and air 
warfare—are the vital missions of the U. S. Naval Ord- 
nance Laboratory near Silver Spring, Md. 

Established in 1918 as a small experimental group 
to design sea mines and fuzes, NoL has become one of 
the nation’s most respected institutions for military re- 
search and development. Three thousand persons, in- 
cluding 1,000 scientists and engineers, work at its 
modern $50 million facility. 

In its forty years, the Laboratory has made many 
notable contributions to science and technology, and 
to our national defense. In World War II, Nou-devel- 
oped sea mines helped to strangle Japan’s supply sys- 
tem. More recently, Not’s work on TARTAR, TERRIER 
and TALOs greatly speeded the development of these 
needed missiles. The atomic depth charge weapon 
BETTY, conceived and perfected by Not in less than 
five years, now provides the U. S. Fleet a massive de- 
fense against modern enemy submarines. 

Today, Nox’s work embraces virtually the entire 
field of advanced weapons design. One highlight is the 


This is 
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one of 
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Laboratory’s present key role in the development of 
POLARIS, the nuclear-capability missile that our sub- 
marines will soon be able to launch while submerged 
off enemy shores. This is one of many new weapons 
that Nox is developing or assisting with at present — 
weapons that are needed if the tree world is to survive. 

But in the long run, some of Not’s other work in 
basic and applied research may be of even more sig- 
nificance to mankind. For example, the advanced aero- 
ballistics studies carried on at NoL, which already have 
solved some of the rBM and 1cBM problems, may help 
unlock the remaining doors to interplanetary space 
travel. Or NoL’s work in advanced chemistry may pro- 
vide entirely new and better types of metals, plastics 
and other substances. Or Now’s basic research in phys- 
ics and mathematics may furnish a key to the internal 
structure of subatomic particles, and open up whole 
new worlds for science and mankind. 

Thus the Naval Ordnance Laboratory aids our 
nation in two ways: by developing weapons that will 
insure the survival of freedom, and by finding new sci- 
entific knowledge that will enrich our lives. — 


or Wetense. 


FORD INSTRUMENT CoO. 
DIVISION OF SPERRY RAND CORPORATION 
31-10 Thomson Avenue, Long Island City 1, New York 
Field Sales Offices: Beverly Hills, Calif.; Dayton, Ohio 


Engineer at Ford Instrument Co. checks out 


fire control computer developed and built by 


‘é Ford for the U.S. Navy's Bureau of Ordnance. 
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CY SHOWER 
EYEWASH FOUNTAIN 


This safety station is part of the compre- 
hensive safety program maintained by 
the General Electric Company at its Ap- 
pliance Park, Louisville, Kentucky. Gen- 
eral Electric considers its safety record 
a prime achievement,and provides Haws 
Emergency Facilities at critical locations 
throughout the production areas. 


HAWS 
EMERGENCY 
FACILITIES 


e Drench Showers ¢ Eye-Face Wash Fountains 


Not just a shower ...a HAWS DRENCH SHOWER! A quick pull of the 
chain releases a solid sheet of water to rid body and clothing 

of dangerous chemicals and caustics. Safety in seconds—saving 
" yital moments until medical aid arrives, perhaps avoiding 
serious injury. HAWS Eye-Wash Fountain, too, is ready to 
flood eyes from specially designed fountain heads. These are 
only two items from HAWS complete line of emergency facil- 
ities to meet every industrial need. Safety authorities stress 
the necessity of instant irrigation for eye or body contamina- 
tion. Provide for it with HAWS equipment. Get the full facts 


now, by writing... 


Leader in drinking water facilities since 1909 
1443 FOURTH STREET « BERKELEY 10, CALIFORNIA 


lite Vehicle Program. 1: 120-121, Nov. 1956. 
Complete organizational and _ functional 
charts for the United States IGY satellite 
program. 

Newell, H. E., Jr. 


The International Geophysical Year Earth 
Satellite Program. Franklin Institute. 
Earth Satellites as Research Vehicles. Pro- 
ceedings of a Symposium, Philadelphia, 
Franklin Institute, Apr. 18, 1956, pp. 3-10, 
Lancaster, Pa., Journal of the Franklin Insti- 
tute, June 1956. (Monograph 2.) 

Article with same title in Sci. Mon. 83: 
13-21, July 1956. 

Odishaw, Hugh 

The IGY Earth Satellite Program. Inst: 
Soc. Am. J. 3: 298-301, Sept. 1956. Obje 
tives of the IGY and methods used to imple- 
ment them. Includes description of pro- 
posed rocket and satellite operation anc! 
instrumentation. 

Ovenden, M. W. 

The Astronomical Uses of Small Artificia 
Satellites. Sci. Prog. 44: 472-479, July 1956 
Summarizes possible astronomical uses of th 
IGY satellites. 

Perkins, F. M. 

Flight Mechanics of Ascending Satellit 
Vehicles. Jet PRopuLsIon 26: 352-358, May 
1956. Solution of the optimum vacuun 
stage thrust axis pitch program for a satellit« 
using continuous thrust all the way fron 
ground to orbit, with an explanation oi! 
trajectory optimization to minimize propel- 
lant consumptions for a vehicle using « 
power-off coast up to the orbit, and of analog 
and digital computor techniques for trajec- 
tory calculations. Includes discussion of 
types of orbits and trajectories. 

Petersen, N. V. 

Factors Affecting the Lifetime of Earth 
Satellites. Aero. Dig. 73: 74-76, 78, 80, 82, 
July 1956. Influence of gravitational field 
and atmospheric envelope of the earth are 
considered. 

Lifetimes of Satellites in Near-Circular and 
Elliptic Orbits. Jet 26: 341 
351, diags., May 1956. The problem of 
determining the approximate lifetime or 
duration for simple satellite configurations as 
influenced by the gravitational field and 
atmospheric envelope of the earth is con- 
sidered. 

Pfister, Wolfgang 

Study of Fine Structure and Irregularities of 
the Ionosphere and Rockets and Satellites. In 
Van Allen, J. A.,ed. Scientific Uses of Earth 
Satellites, pp. 283-291, Ann Arbor, Univer- 
sity of Michigan Press, 1956. Real advan- 
tage of the experiment will be that the chang- 
ing aspect allows scientists to gain a three- 
dimensional picture of the configuration of 
the blob, while so far they are used to looking 
only at the two-dimensional diffraction pat- 
tern at the ground. 


Porter, R. W. 

Recovery of Data in Physical Form. In 
Franklin Institute. Earth Satellites as Re- 
search Vehicles. Proceedings of a Sympo- 
sium, Philadelphia, Franklin Institute, Apr. 
18, 1956, pp. 103-112, Lancaster, Pa., Jour- 
nal of the Franklin Institute, June 1956. 
(Monograph 2.) The chairman of the Panel 
on Earth Satellites outlines the difficulties of 
recovering data, possible methods of attack, 
proposed design, trajectory computations and 
thermal considerations. 


Rand Corp., Santa Monica, Calif. 

Skin Temperature Variation During Re- 
entry of Scientific Satellite, by D. J. Masson. 
4 pp., Mar. 30, 1956. (Res. Memo 1693.) 
From the point of view of destruction by 
heat, the recovery of a satellite is considered 
possible. 

Rosen, M. W. 

Placing the Satellite in Its Orbit. Inst. 

Radio Engrs. Proc. 44: 748-751, June 1956. 
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Here at Douglas we’re involved in a greatly 
accelerated missile and space program. This 
requires one of the most intensive engineering 
and research efforts in our history. 

The problems are great ones as we move into 
the new dimension of unmanned and manned 
space vehicles. They require specialists in almost 
every engineering field. But their solution will 


May 1958 


STRAIGHT TALK TO ENGINEERS 
from Donald W. Douglas, 


ji ai President, Douglas Aircraft Co., Inc. 


Jn 


result in great benefits not only to our own nation 
but to all mankind. 

If you’re interested in tackling these problems 
with us...in giving your best in an all-out drive 
to solve them... we’re interested in you! 

Please write to Mr. C. C. La Vene 
Douglas Aircraft Company, Box J-620 
Santa Monica, California 


ce 
— 


AEROJET-GENERAL 
CORPORATION 


TURBO-MACHINERY 
DIVISION 


Offers Key Positions 


for specialists, design engineers 
and test engineers with creative 
ability and broad. experience in 
cycle and performance analysis, 
design, and development in fol- 
lowing fields: 


Propellant-powered, high-speed 
machines from one to several thou- 
sand horsepower, operating at 
high pressure ratios and tempera- 


tures. 


Axial and radial flow pumps of 
high suction specific speed for 
mono-propellant and bi-propellant 
systems. 


Mono-propellant and bi-propellart 
systems for wide ranges of design 
and operating conditions. 


¥ ELECTRICAL 
q 


Frequency, voltage, sequencing, 
and other electrical controls for 
turbo-machines, auxiliary power 
and magnetic amplifier systems, 
and related ground support equip- 
ment. 


You are cordially invited 
to send a detailed resume to: 


AEROJET-GENERAL 
CORPORATION 


Box 296-H 


L. THOMPSON 


Azusa, California 


A Subsidiary of The 
General Tire & Rubber Company 


Orbits and Lifetimes of Minimum Satellites. 


Description of Vanguard launching vehicle. 
Republished in J. Astronautics 3: 61—64, 82, 
illus., Autumn/Winter 1956. 
Saelman, B. 

Note on an Artificial Satellite. 
Phys. 24: 630, Dec. 1956. 
Sedwick, J. L., Jr. 

Interpretations of Observed Perturbations 
on a Minimal Earth Satellite. In Van Allen, 
J. A., ed. Scientific Uses of Earth Satellites, 
pp. 44-48, Ann Arbor, University of Michigan 
Press, 1956. Gives a method for calculating 
the perturbing forces from observations of 
the orbit, considering air resistance and 
effects of the earth’s oblateness. 

Shternfel’d, A. 

Iskusstvennye 
Earth Satellites). 


Am. J. 


Sputniki Zemli (Artificial 
180 pp., Moscow, State 
Publishing House of Technicotheoretical 
Literature, 1956. In Russian. Contains 
information about propellants for outer space 
rockets, satellites and ballistic missiles. 
Sicinski, H. S., Spencer, N. W., and Boggess, 
R. L. 

Pressureand Density Measurements Through 
Partial Pressures of Atmospheric Components 
at Minimum Satellite Altitudes. In Van 
Allen, J. A., ed. Scientific Uses of Earth 
Satellites, pp. 109-118, Ann Arbor, Univer- 
sity of Michigan Press, 1956. Pressure and 
density measurements by means of a ‘“‘Syn- 
chrometer.”’ 

Singer, S. F. 

Application of an Artificial Satellite to the 
Measurement of the General Relativistic ‘‘Red 
Shift.”” Phys. Rev. 104: 11-14, Oct. 1, 1956. 
Some experimental schemes for the measure- 
ment of the clock shift are discussed. 

The Artificial Earth Satellite. Discovery 
17: 140-145, Apr. 1956. <A description and 
plan of MOUSE are given. 

The Artificial Earth Satellite—Past, Present 
and Future. Interavia 11: 956-959, illus., 
Dec. 1956. How satellite experiments may 
lead to a better scientific understanding of 
our environment which may have important 
consequences for our economy and way of 
life. 

The Effect of Meteoric Particles on a Satel- 
lite. Jet PRoputsion 26: 1071-1075, Dec. 
1956. Discusses some physical models to 
allow calculations of the erosion rate of high 
velocity meteoric dust on space vehicles. 
Important indications are derived for the 
choice of material and construction of a 
satellite skin to minimize meteor effects. 
An erosion experiment for the Vanguard 
satellite is described. 

Geophysical Research with Artificial Satel- 
lites. In Landsberg, H., ed. Advances in 
Geophysics, pp. 301-367, New York, Aca- 
demic Press, 1956. A review of experiments 
which have been suggested for a satellite, a 
geophysical tool whose scope extends further 
than any other known technique. 

Measurements of Interplanetary Dust. In 
Van Allen, J. A., ed. Scientific Uses of 
Earth Satellites, pp. 301-316, Ann Arbor, 
University of Michigan Press, 1956. Indi- 
cates that a study of interplanetary dust can 
give a great deal of new information on the 
electromagnetic conditions of interplanetary 
space. 

Measurements of the Earth's Magnetic 
Field from a Satellite Vehicle. In Van Allen, 
J. A., ed. Scientific Uses of Earth Satellites, 
pp. 215-233, Ann Arbor, University of 
Michigan Press, 1956. Discusses the applica- 
tion of satellite observations to magnetic 
measurement on three different time scales. 

Satellites for Physicists. Phys. Today 9: 
21-23, diags., Apr. 1956. Points out the 
educational and missionary aspects of satel- 
lites as well as other possibilities with which 
the geophysicist and the astrophysicist are 
concerned. 

Studies of a Minimum Orbital Unmanned 
Satellite of the Earth (MOUSE). Part II. 


MARTIN 
DENVER 


has created an Applied Research 
Section for which the following 
specialists are required. Edu- 
cation requirements in all cases 
are a Ph.D. Experience should 
be heavy enough to permit in- 
dependent research activity with 
a minimum of supervision: 


| 
| 


Solid State Physicist to work on 
electronic devices. 


Nuclear Physicist with experi- 
ence in nuclear propulsion and 
general applications of nuclear 
power. 


Applied Mathematician with 
experience in computers, aero- 
dynamics, trajectories, etc. Must 
be a senior man. 


Magneto-Hydrodynamicist to do 
developmental work ad- 
vanced propulsion systems. 


Information Theory and Noise 
Analysis Engineer experienced 
mathematical analysis and 


in 

theoretical electronics. 

. Engineer experienced in ad- 
vanced propulsion techniques 


such as ion motors, free radicals, 
etc. 


Physicist experienced in re-entry 
field, familiar with problems of 
radiation from high temperature 
gases, hypersonic flow, etc. 


Qualified persons are urged to 
write, giving complete details 
of education, experience (in- 
cluding personal technical ac- 


complishments) and salaries 
earned to: 
DAVID POTTER H-10 DIV. 
MARTIN-DENVER 
P.O. BOX 179 
DENVER 1,COLORADO 
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Astronautica Acta 2: 125-144, 1956. Pur- 
pose is, first, to describe the (elliptical) orbit 
of an artificial satellite in terms of the launch- 
ing conditions (i.e., at burnout); in particu- 
lar, to give perigee (minimum) and apogee 
naximum) altitudes of the ellipse as a func- 
tion of launching altitudes and of errors in 
launching velocity and angle. Second, to 
investigate the subsequent behavior of the 
orbit under the influence of the drag of the 
upper atmosphere. 
Simas, V. R. 
Low-Noise Preamplifier for Satellite 
racking. QST 40: 42-43, 172, 174, 176, 
illus., Dee. 1956. Describes a preamplifier 
nd converter fashioned for use in the Mini- 
rack system for tracking the earth satellite. 
Sixth Annual Congress of IAF 
Abstracts of Papers. J. Astronautics 3: 
5-16, 25, Spring 1956. Includes the follow- 
g selected abstracts: Uses of Satellite 
ehicles, by H. E. Canney and F. I. Ordway; 
‘he Visibility of an Earth Satellite, by R. 
‘ousey (see also entry under Tousey); The 
satelloid, by K. A. Ehricke; Some Problems 
n Rocket Development, by G. A. Partel; 
‘stimated Lifetimes of Satellites from Near 
Circularity, by N. V. Petersen. 


Sky and Telescope 

Layout and Organization of a Moonwatch 
Station. 15:3-4(follows p. 552), Oct. 1956. 
Spitzer, Lyman Jr. 

Astrophysical Research with an Artificial 

Satellite. In Franklin Institute. Earth Satel- 
lites as Research Vehicles. Proceedings of a 
Symposium, Philadelphia, Franklin Insti- 
ute, Apr. 18, 1956, pp. 69-74, Lancaster, 
Pa., Journal of the Franklin Institute, June 
1956. (Monograph 2.) Lists and dis- 
‘usses, briefly, problems that could be inves- 
tigated such as: Earth’s atmosphere, plan- 
etary atmospheres, matter in space, and 
tellar atmospheres. 

On the Determination of Air Density from 
t Satellite. In Van Allen, J. A., ed.  Sciten- 
tific Uses of an Earth Satellite, pp. 99-108, 
Ann Arbor, University of Michigan Press, 
1956. Discusses various possible methods of 
evaluating the atmospheric density, including 
observation of the satellite’s orbit, and 
measurement of the change of spatial orien- 
tation. 

Stakutis, V. J., and Brennan, J. X. 

Visibility from a Satellite at High Altitudes. 
In Van Allen, J. A., ed. Scientific Uses of 
Earth Satellites, pp. 137-146, Ann Arbor, 
University of Michigan Press, 1956. The 
simple design of two types of instrumentation 
for measuring luminous quantities is de- 
scribed to show their applicability to satellite 
installation. 

Stedfield, R. L. 

Engineering the Earth Satellite. Mach. 
Design 28: 82-86, Nov. 1, 1956. Descrip- 
tion of the orbit, the launching rocket and 
of the Vanguard satellite itself and its 
instrumentation. 


Steel 

Metalworking Tackles Satellite. 139: 120- 
123, Oct. 15, 1956. Steps to be employed in 
manufacture of shell—building of aluminum 
magnesium skinned launching vehicle which 
will be a three-stage rocket about 72-ft long, 
Steier, H. P. 

Tracking the IGY Satellite. Missiles and 
Rockets 1: 76-77, Oct. 1956. Details of the 
transistorized subminiature transmitter called 
the Minitrack. 

What Guides the Vanguard? Missiles and 
Rockets 1: 70-72, Nov. 1956. Minneapolis- 
Honeywell builds complex gyro reference 
system. Description and photos included. 
Stroud, W. G. and Nordberg, W. 

Meteorological Measurements from a Satellite 
Vehicle. In Van Allen, J. A., ed. Scientific 
Uses of Earth Satellites, pp. 119-132, 
Ann Arbor, University of Michigan Press, 
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A quiz for 
Rocket Propulsion Engineers 


7 (Time Limit: 10 years) 


: 1. What is the maximum theoretical specific impulse obtainable 
_ for a rocket motor with zero combustion volume? 
; = What is the basic mechanism involved in jet separation in 
nozzles, and how can the separation point be predicted as a func- 
tion of nozzle divergence angle? 
3. What is the effect on turbine blade design of a working gas 
: which is still reacting and changing its composition? 
4, What are the dynamic interrelations among the various com- 
7 - ponents of a rocket engine, and what control systems will best 
_ solve the rocket engine control problems? 
_ 5. What is the flame holding mechanism in a rocket motor and 
how can fundamental understanding of such a mechanism be 
used to increase reaction ratio? 
7 6. What types of fluid transport systems offer simplicity, effi- 
ciency, and light weight (other than the conventional 


° turbopump)? 

° 7. What are some workable methods of translating chemical 
: energy into thrust — other than by the standard combustion at 
° high pressure and then conversion to velocity? 

: If you would like to work in these areas and have an outstand- 
. ing technical background, an inquiring mind, and the creative 
: ability to visualize new solutions, possess a M.S. or Ph.D. degree 
. in Physics, Chemistry, M.E., A.E., Ch.E., E.E., or a B.S. degree 
. in the above with a thorough background in applied research, we © 
have a place for you at Rocketdyne. 
° Please write to Mr. A. E. Jamieson, Engineering Personnel, 

° 6633 Canoga Avenue, Canoga Park, California. oer 
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VARIABLE FREQUENCY 
MOTOR GENERATOR SETS 
ADJUST 360 TO 440 CPS. 


Generator mounted controls 
include reset buttons, limit 
switch. Motor and generator 
remain stationary. Variable 
drive pulley adjustment con- 
trolled by small motor. Remote 
control panels available 


PRECISE 
CURRENT 


for your 
TESTING 
NEEDS! 


Units can be equipped with 
synchronous motor starter and 
magnetic amplifier, auto- 
matic voltage regulator. 


KATO 400 CYCLE 
MOTOR GENERATOR SETS 
NOW UP TO 250 KWI 
KATO MOTOR GENERATOR 


SETS are available in frequencies, 


speeds and sizes for every special- INPUT 

ized use . . . operating high cycle 60 CYCLES 
tools, testing components an OUTPUT 
electronic equipment. 400 CYCLES 


WRITE FOR NEW FOLDERI 


Builders of Fine Electrical Machinery Sinee 1928 


1498 First Ave., Mankato, Minn. 


1956. Usefulness of observations is reviewed 
and various levels of information available 
from various levels of complexity of instru- 
mentation are outlined. A specific instrumen- 
tation suitable for the initial IGY satellites 


| is presented and basic aspects of operational 


| problems outlined. 


It is shown that the 
necessary orientation information can be ob- 
tained from the albedo data. 
Symposium on the U. S. Earth Satellite 
Program—Vanguard of Outer Space 

Inst. Radio Engrs. Proc. 44: 741-767, 
illus., June 1956. See entries under Hagen, 
Rosen, Mazor, Mengel, Whipple and Van 
Allen for papers given at this symposium. 


Tousey, R. 


The Visibility of an Earth Satellite. Astro- 
nautica Acta 2: 101-112, 1956. Conditions 


for naked eye and telescopic visual observa- 
tion of white 21-in.-diam earth satellite are 
calculated using known values of sky bright- 
ness, solar illumination and visual thresholds 
of eye. See also entry under Sixth Annual 
Congress. 
U. S. National Committee for the Inter- 
national Geophysical Year 

USNC-IGY Rocketry and Satellite Pro- 
grams. 41 pp., Washington, National Acad- 
emy of Sciences, National Research Council, 
1956. Collection of five articles: The IGY 
Rocket and Satellite Program, by Joseph 


| Kaplan; The U. S. IGY Rocketry Program, 


| by H. E. Newell Jr.; Observations with 
| Satellite-borne Apparatus, by J. A. Van 
Allen; Radio Tracking—Earth Satellite 


Program, by J. P. Hagen; The Optical Pro- 
gram, by F. L. Whipple. 
Van Allen, J. A. 

The Artificial Satellite as a Research Instru- 
ment. Sci. Am. 195: 41-47, Nov. 1956. 


Cosmic-Ray Observations in Earth Satellites 
In Van Allen, J. A., ed. Scientific Uses of 
Earth Satellites, pp. 171-187, diags., Ann 
Arbor, University of Michigan Press, 195«. 
A single Geiger tube or scintillator carried in 
a satellite will make possible the study of the 
cosmic-ray intensity above the atmosphere on 
comprehensive geographical and temporal 
bases for the first time. 

Ed.: Scientific Uses of Earth Satellites. 
316 pp., illus., Ann Arbor, University of 
Michigan Press, 1956. Includes thirty-three 
papers given at the Upper Atmosphere 
Rocket Research Panel at the University of 
Michigan on Jan. 26-27, 1956. Principa! 
purpose of the book is, according to the edi- 
tor’s preface, to bring potential value of 
artificial satellites to the attention of the 


scientific community at large. See also 
entries under individual authors: Bennet 
Berning, Chubb, Daniels, Davis, deBey, 


Dubin, Gast, Hartman, Heppner, Hintereg 
ger, Hok, Hudson, Jones, Katz, LaGow 
Merrill, Pfister, Sedwick, Sicinski, Singe: 
Spitzer, Stakutis, Stroud, Van Allen, Vestine 
Ziegler and Zirker. 

Scientific Value of the Earth Satellite Pro- 
gram. Inst. Radio Engrs. Proc. 44: 764-767, 
June 1956. Highest ‘‘flight-priority” in the 
Vanguard Project has been assigned to the 
following: (a) Monitoring of the intensity of 
the solar ultraviolet, (b) monitoring of cosmi« 
ray intensity and the measurement of its lati- 
tude, longitude, and altitude dependence, (« 
measurement of the size spectrum and the 
number density of interplanetary dust, and 
(d) measurement of the earth’s optical albedo 
over large area. 

Study of the Arrival of Auroral Radiations. 
In Van Allen, J. A., ed. Scientific Uses of 
Earth Satellites, pp. 188-193, diags., Ann 


PHYSICISTS 


MONOBALL” 


Self-Aligning Bearings Research and Development 


Outstanding theoretical and experimental physi- 
cists needed for research programs in our Aero- 


ROD END 


PLAIN TYPES vicina physics group. Challenging research positions 
al are availabe in the following fields: oi 
PATENTED USA 


All Werld Right Reserved 


_ Shock and Detonation Wave Phenomena 

Hypersonic Flow 
_Rarefied Gas Phenomena 
_ Transport Phenomena 


CHARACTERISTICS 


RECOMMENDED USE 


ANALYSIS 


Stainless Steel For types operating under high temper. 

a Ball and Race ature (800-1200 degrees F.). sy Plasma Physics 

Chrome Alloy For types operating under high radial 
Stee! Balland Race —( ultimate loads (3000-893,000 Ibs.). If you have an M.S. or Ph.D. degree and experi- 


ence or interest in these fields, we offer you an 
opportunity to use your initiative and creative 
ability. 


Bronze Race and [ Seusepee operating under normal loads 
Chrome Steel Ball with minimum friction requirements. 
Thousands in use. Backed by years of service life. Wide variety 
of Plain Types in bore sizes 3/16” to 6” Dia. Rod end types in 
similar size range with externally or internally threaded shanks. 
Our Engineers welcome an opportunity of studying individual 
requirements and prescribing a type or types which will serve 
under your demanding conditions. Southwest can design special 
types to fit individual specifications. As a result of thorough 
study of different operating conditions, various steel alloys 
have been used to meet specific needs. Write for revised Engi- 
neering Manual describing complete line. Dept. JP-58. 


SOUTHWEST PRODUCTS CO. 


OUNTAIN AVE., MONROVIA RNIA 


Excellent employee benefits including liberal 
vacation policy. Please send resume to: 


E. P. Bloch 


ARMOUR RESEARCH FOUNDATION 
of Illinois Institute of Technology 
10 West 35th Street 

Chicago 16, Illinois 
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lites Arbor, University of Michigan Press, 1956. Winterberg, F. 
ses of A satellite in a nearly pole-to-pole orbit is Relativistische Zeitdilatation Eines Kuenst- 
Ann considered a splendid vehicle for plotting out lichen Satelliten (Relativistic Time Dilatation ae a a 
1956, auroral zone by direct observation. in an Artificial Satellite). Astronautica Acta i ).> ‘fe 
ied in Vestine, E. H. 2: 25-29, 1956. In German. Applies the aa 45 Saad 
of the Baploring the Atmosphere with a Satellite- general theory of relativity to the problem of th t / t d = 
are on Borne Magnetometer. In Van Allen, J. A., time measurement in a satellite. 0 e tatente ; 
poral ed. Scientific Uses of Earth Satellites, pp. Winternitz, P. F. . . . 
198-214, diags., Ann Arbor, University of The Physical and Chemical Fundamentals engineer & scientist 
lites. Michigan Press, 1956. Summary of explora- of Satellite Flight: Part I and II. J. Astro- 
ty of tions for consideration. Various electric nautics 3: 43-47; 65, 68-70, 80, Spring, 
three current systems suggested for explaining the Autumn-Winter 1956. Physical principles of 
phere known geomagnetic variations at ground rocket operation, chemical principles, applica- A PL 
ity of level are briefly described. tion to the design of rocket power plants for OFFERS 7 
cipal Webb, J. A satellites, computation of circular and escape 
Interplanetary Communications. J. As- period of revolution of a satellite, | CREATER 
: the tronautics 3: 29-30, 36, illus., Summer 1956. _—— | 
Considers, in part, an earth satellite radio H. K. | 
— relay link. ‘omponents for Instrumentation of Satel- 
nett, Whi le. F.L lites. In Van Allen, J. A., ed. Scientific FREEDOM OF 
oe Scientific Value of Artificial Satellites, Uses of Earth Satellites, pp. 55-67, Ann Arbor. 
ereg- University of Michigan Press, 1956. En- 
Gow Frank. J. 262: vironmental conditions are reviewed and | ACTIVITY 
stine, cia "GY power sources. Present and _ predictable 
Pro- pes batteries, electron tubes, transistors and APL has responsibility for 
“460, Whipple, F. L., and Hynek, J. A. frequency control devices is discussed- the technical direction of 
the Research Program Based on the Optical 
Zirker, J. B., Whipple, F. L. and Davis, R. J. h of th ded l 
— l'racking of Artificial Earth Satellites. Inst. Time Available for the Optical Observation much OF the guided missile 
ae ltadio Engrs. Proc. 44: 760-764, June 1956. of an Earth Satellite. In Van Allen, J. A., ed. program of the Navy Bureau 
and Scientific Uses of Earth Satellites, pp. 23-28, of Ordnance. As a result 
ae sible photographic techniques for observing a A Arbor, University of Michigan Press, | . . 
lithe 1956. Light curves during a twilight observa assignments of challengin 
ee White, E. C., Foley, J., and Wilkins, R. G. tion are computed for specularly and dif- ging a 
baa Power Supplies and Telemetry for an fusely reflecting spheres of 20-in. diam. scope that range from basic ‘ 
Instrumented Artificial Satellite. Part I. The specular reflector is visible over a larger research to prototype test- 
Kore Orbital Consideration. Part IT. Instrumen- fraction of its path. The time available for ing of weapons and weapons 
be of tation and Telemetry. Part III. Avail- observation of the specular reflector is com- systems 
Ann ability of Power. Brit. Interplan. Soc. J. puted for several cases. Sky brightness, - 
i 15: 177-191, July-Aug. 1956. geometry, and the earth’s shadow are factors. ’ 
A high degree of freedom 
of action enables APL staff 
members to give free rein 
to their talents and ideas. 
Thus, professional advance- 
ment and opportunities to 2 
accept program responsi- 
bility come rapidly. Promo- - 
tion is rapid, too, because : 
of our policy of placing 
professional technical men 
ENGINEER, ME at all levels of supervision. 
FLUTTER & VIBRATION 
— ments include: the first ram- 
jet engine, the Aerobee high 
A N A LY S T altitude rocket, the super- 
7 . 
a Senior Position With General Electric sonic Terrier, Tartar, and 
Talos missiles. Presently the 
Several years experience in vibration anal- pioneering effort to achieve nuclear Laboratory is engaged in 
ysis of turbo machinery assemblies — powered flight will find themselves in the solving complex and ad- 
plus a strong background in fundamental _ forefront of their profession. vanced problems leading to 
concepts and applications in this field—are Is Nuclear Experience Necessary? No. future weapons and weap- 
essential for this advanced assignment. In-plant courses and graduate study —on ons systems vital to the 
M.S. desirable. full tuition refund basis—provide essential national security. Interested 
The work concerns the development of nuclear theory and technology. engineers and physicists are 
nuclear power systems for aircraft at Gen- Publication of technical papers is encouraged invited to address inquiries 
eral Electric. It deals specifically with high- Location: Cincinnati, Ohio, known both as to: 
temperature reactor core components. an engineering center and a fine place to 
Engineers who have contributed to the live. Excellent housing available. Professional Staff 
Appointments 
7 Please write in confidence, stating salary requirements to: % 
Mr. P. W. Christos 
Tl 
. AIRCRAFT NUCLEAR PROPULSION DEPT. e Johns Hopkins University 
| a GENERAL ELECTRIC Applied Physics Laboratory 
8617 Georgia Avenue, Silver Spring, Md. 
Div. 34-MQ P. O. Box 132 Cincinnati 15, Ohio 
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The first 
volume of an 
important 
new work.. 


AIRCRAFT 
AND MISSILE 
PROPULSION 


Volume 1 


Thermodynamics of Fluid Flow and 
Applications to Propulsion Engines 


By M. J. ZUCROW 
Purdue University 


This is the initial volume of a three- 
volume work designed to give you a bet- 
ter understanding of the fundamental 
principles governing the functioning and 
operating characteristics of engines used 
to propel winged aircraft and missiles, 
both guided and unguided, at high 
eds. Volume I discusses the thermo- 
ynamics of fluid flow and its application 
to propulsion engines. 


These original features help you follow the 
subject matter with ease: 


The nomenclature for a chapter precedes the 
chapter. 


All derivations are explained in detail. 


An unusually large number of i!lustrative 
problems are worked completely. — 


Numerous exercise problems are included 
with their answers. 


Tables or figures needed for solving problems 
are provided. 


® References at the end of each chapter guide 
you in further study. 


1958 538 pages Illus. $11.50 


Mail this coupon for your ON-APPROVAL 
copy TODAY! 


JOHN WILEY & SONS, Inc. JP-58 
440 Fourth Ave., New York 16, N. Y. 
Please send me a copy of AIRCRAFT AND MISSILE 
PROPULSION, Vol. I, to read and examine ON AP- 
PROVAL. Within 10 days I will return the book and 
owe nothing, or I will remit $11.50, plus postage. 


...Zone.. ..State.. 
SAVE POSTAGE! Check here if you pay- 
ment, in which case we pay postage. Same return 
privilege, of course. 
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